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-p  ■  1.  SUMMARY 

:  r  .-> 

“The  project  on  "Finite  Elements  and  Localized  Failure"  was  comprised 
of  research  on  three  different  levels: 

^‘(i ) ^Experimental  Work: 

Stroke  controlled  post-peak  experiments  on  plain  concrete 
specimens  which  were  subjected  to  direct  tension,  triaxial 

compression  and  direct  shear. 

J 

0  (iij ^Constituti ve  Work: 

Homogenization  of  localized  failure  due  to  tensile  cracking  and 
frictional  slip  in  the  form  of  the  "Composite  Fracture  Model" 
which  describes  the  degradation  of  strength  within  an  elementary 


volume  of  a  fracture  energy  equivalent  continuum, 
(iii ) -'Computational  Work: 
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Stabilization  of  strain-softening  computations  of  localized 
failure  within  solids  subjected  to  tensile  cracking  and/or  shear 
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2.  RESEARCH  OBJECTIVES 


The  principal  objective  of  the  research  project  on  "Finite  Elements 
and  Localized  Failure"  was  the  limited  ductility  or  rather  the 
"fragility"  of  particulate  composites  such  as  concrete  under  low 
confinement.  This  facet  of  the  material  behavior  controls  the  strength 
and  deformabi 1 ity  of  indeterminate  structures  which  exhibit  substantial 
reserves  as  long  as  sufficient  ductility  accompanies  the  localization  of 
failure. 

Within  this  endeavour  several  critical  issues  had  to  be  addressed  in 
order  to  gain  insight  into  the  fundamental  aspects  of  brittle-ductile 
fracture.  Since  failure  is  a  structural  rather  than  a  local  constitutive 
process  involving  stress  redistribution  and  non-uniform  deformation,  the 
research  efforts  enconpassed  several  closely  interrelated  tasks  on  the 
experimental,  the  constitutive  as  well  as  the  computational  level. 

2. 1  Experimental  Issues 

1.  Under  what  condition  does  concrete  failure  reduce  to  a  local 
constitutive  issue  such  that  the  degradation  of  strength  and 
stiffness  on  the  structural  level  does  directly  reflect  the 
effects  of  localized  strain-softening? 

2.  To  what  extent  is  progressive  concrete  failure  a  surface  driven 
fracture  phenomenon  of  singular  macrodefects  rather  than  the 
result  of  distributed  microdefects  within  a  volume  driven  damage 
process? 

3.  Is  there  a  transition  between  brittle  fracture  and  ductile 
flow  which  separates  softening  from  the  hardening  response 
regime  and  where  is  it  located? 
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4.  Are  tensile  cracking  and  shear  faulting  separate  failure 
mechanisms  due  to  debonding  and  decohesion  or  is  the  development 
of  shear  bands  simply  the  result  of  coalescing  tensile  cracking? 

5.  Is  dilatancy  a  frictional  mechanisms  associated  with  the  loss  of 
normality  or  is  it  simply  another  facet  of  fracture  under  mixed 
mode  conditions? 

6.  To  what  extent  is  brittle  failure  the  result  of  tensile 
straining  rather  than  tensile  stressing? 

7.  Is  the  loss  of  strength  in  the  post-peak  regime  accompanied  by 
simultaneous  degradation  of  stiffness  which  in  turn  introduces 
strong  elastic-plastic  coupling? 

2.2  Constitutive  Issues 

1.  Can  the  differences  of  discrete  and  distributed  failure  analysis 
be  reconciled  by  incorporating  fracture  energy  concepts  into  the 
strain-softening  formulation  of  "equivalent"  continua? 

2.  Can  a  "nonlocal"  constitutive  formulation  be  developed  which 
captures  the  transition  of  failure  from  tensile  cracking  to 
shear  faulting  which  gradually  diffuses  into  continuous 
hardening  under  increasing  confinement? 

3.  Is  failure  an  anisotropic  phenomenon  which  requires  storage  and 
monitoring  of  vector-  or  tensor-  valued  state  variables  instead 
of  scalars  only? 

4.  Can  the  elastic-plastic  coupling  be  neglected  as  a  first 
approximation  or  do  we  have  to  account  for  the  simultaneous 
degradation  of  strength  as  well  as  stiffness? 
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5.  Is  there  a  necessity  for  rate  sensitive  constitutive 
formulations  such  as  viscoelasticity  or  viscoplasticity  in  order 
to  assure  positive  definiteness  of  the  underlying  constitutive 
operators? 

6.  Does  the  strain-space  setting  of  strain-softening  plasticity 
resolve  the  difficulties  of  stable  peak  and  post-peak 
descriptions  or  is  the  latter  simply  an  issue  of  the  driving 
control  variables? 


2. 3  Computational  Issues 

1.  Can  the  computational  difficulties  of  discrete  fracture  analysis 
with  continuously  changing  fracture  topologies  be  resolved  by 
"equivalent"  strain-softening  formulation  of  the  constitutive 
behavior? 

2.  What  is  the  basic  difference  between  the  proposed  composite 
fracture  model  and  the  discrete  crack  analysis  in  fracture  mech¬ 
anics  with  regard  to  fracture  energy  release  rate  computations? 

3.  To  what  extent  do  progressive  failure  computations  within  the 
proposed  composite  fracture  model  circumvent  the  topological 
difficulties  of  the  fictitious  crack  model  of  Hillerborg  et  al. 
and  those  of  the  imbricate  continuum  theory  by  Bazant  et  al.? 

4.  What  are  the  restrictions  for  uniqueness  and  stability  of 
contained  strain-softening  computations  with  regard  to  the 
maximum  mesh  size? 


5.  What  are  the  numerical  requirements  for  capturing  localization 


and  propagation  of  contained  strain-softening  regions  within  a 
sol  id? 
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6.  To  what  extent  can  post-peak  response  predictions  of 

strain-softening  computations  be  stabilized  by  displacement  or 
arc-length  control  of  the  underlying  incremental/iterative 
solution  strategy? 
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3.  SIGNIFICANT  ACCOMPLISHMENTS 


The  diverse  research  efforts  on  "Finite  Elements  and  Localized 
Failure"  resulted  in  several  significant  accomplishments  which  are 
summarized  below  within  three  groups  of  activities: 

3. 1  Experimental  Results 

Three  different  sets  of  experiments  were  carried  out  to  examine  the 
post-peak  response  behavior  of  plain  concrete  specimens  subjected  to 
direct  tension,  triaxial  compression  and  direct  shear  loading  histories. 
All  experiments  were  carried  out  with  a  servo-controlled  MTS-testing 
system  in  order  to  study  the  load-displacement  response  and  fragility 
behavior  of  concrete  failure  under  low  confinement  and  the  localization 
of  failure  regions  when  the  specimen  was  loaded  into  the  post-peak 
regime. 


a.  Direct  Tension  Experiments 

The  principal  results  of  these  experiments  are  illustrated  in  Figs. 
1,  2,  3  and  4.  The  response  curves  depict  the  axial  force-displacement 
(stroke)  behavior  when  NX-size  cylindrical  concrete  specimens  of  D  =  54 
mm  diameter  and  H  =  108  mm  height  were  subjected  to  monotonic  loading 
histories  with  partial  unloading. 

The  experimental  observations  support  the  following  conclusions: 

(i)  Concrete  failure  in  direct  tension  is  an  extremely  brittle 
fracture  phenomenon  which  is  difficult  to  control  in  a  stable 
fashion  if  unnotched  specimens  are  utilized. 

(ii)  Tensile  cracking  is  clearly  a  fracture  energy  dominated  process 
in  which  the  energy  release  per  unit  of  newly  generated  crack 
surface  plays  a  fundamental  role.  In  fact  the  direct  tension 
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experiments  with  unnotched  cylindrical  concrete  specimens  of 
different  height  clearly  support  the  hypothesis  of  a  constant 
fracture  energy  release  rate  in  tension,  i.e.  =  const. 

(iii)  The  degradation  of  tensile  strength  during  debonding 

consistently  follows  an  exponential  decay  function  when  the 
nominal  tensile  stress  is  plotted  versus  the  overall  extension 
of  the  tension  specimen. 

(iv)  The  tensile  fracture  mechanism  fully  supports  the  formation  of  a 
singular  macrodefect  which  is  oriented  perpendicular  to  the 
direction  of  major  principal  stress,  while  the  contiguous 
regions  of  the  specimens  unload  elastically.  As  a  result,  the 
overall  extension  of  the  tension  specimens  furnishes  a  global 
measure  of  the  crack  opening  displacement  during  tensile 
debonding. 

(v)  The  degrading  bond  strength  is  fully  recoverable  under  cyclic 
loading-unloading  conditions.  Separation  takes  place  along  an 
extremely  tortous  crack  surface  with  strong  "aggregate 
interlock"  in  direct  tension. 

(vi)  There  is  considerable  degradation  of  tensile  stiffness  in  the 
post-critical  response  regime,  i.e.  the  unloading-reloading 
branches  indicate  a  progressive  reduction  of  the  slope  near  zero 
residual  strength. 

( v i i )  Direct  tension  experiments  on  cylindrical  test  specimens  of 

different  height  and  constant  cross-section  support  the  fracture 
energy  release  concept  of  a  single  macrodefect  rather  than  the 
damage  concept  of  distributed  microdefects.  Representative 


values  for  the  fracture  energy  release  rate  of  low  strength 
concrete  with  r  =  22.1  MPa  are  in  the  range  of  =  60  N/m 
for  full  height  and  half height  specimens 

I  3^f 

Gi  ~  -  — -  >  -  f  a.  du,  =  56  to  62  N/m  (1) 

)  ^  A  'll  Z 
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b.  Triaxial  Compression  Experiments 

The  principal  results  of  these  experiments  are  illustrated  in  Figs. 
5,  6,  and  7  which  depict  the  axial  force-displacement  (stroke)  and  the 
lateral  versus  the  axial  displacement  (stroke)  response  when  NX-size 
cylindrical  concrete  specimens  were  subjected  to  confined  triaxial 
compression.  The  experimental  observations  support  the  following 
conclusions : 

(i)  The  strength  and  dictility  increase  dramatically  with  increasing 
confining  stress. 

(ii)  There  is  a  distinct  transition  from  brittle  to  ductile  fracture, 
whereby  the  softening  failure  mechanism  gradually  transforms 
into  continuous  hardening, 

(iii)  The  partial  unloading  branches  during  axial  loading  in  the  pre- 
and  post-peak  regime  indicate  surprisingly  little  degradation  cf 
sti ff ness . 

(iv)  The  lateral  deformations  illustrate  the  strong  frictional 

dilatation  which  diminishes  rapidly  with  increasing  confinement, 

(v)  The  failure  mode  changes  gradually  from  "compressive  splitting" 
under  uniaxial  compression  to  inclined  shear  faulting  under 
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Dilatation  Response  of  Convencional  Triaxiai  Compression 
Tests  with  Partial  Unloading 


increasing  confinement.  Beyond  the  transition  point  failure 
transforms  slowly  into  diffuse  damage  within  the  entire 
specimen. 

(vi)  The  modulus  cf  toughness  provides  a  first  approximation  of  the 
strain  energy  release  during  compressive  splitting  in  unconfined 
compression  when  the  energy  dissipation  in  the  hardening  regime 
is  discarded 


c  r 

Uf  a  /  pz  duz  “  -10*91  (2) 

UP 

In  view  of  the  compressive  splitting  mode  of  failure,  the  crack 
surface  area  is  ill-defined.  However,  if  we  assume  that  the 
fracture  energy  release  rate  6^  =  60  N/m  controls  the 
splitting  mode  of  failure  the  unknown  crack  surface  area 
can  be  evaluated  from 

I  aUf 

G*  = - 1  ♦  A  a  0.18  m2  (3) 

f  3Ac  c 

This  value  is  almost  hundred  fold  the  cross-sectional  area  of 
the  cylindrical  test  specimen  which  corresponds  to  the  crack 
surface  area  in  the  direct  tension  test 

A  =  A^  =  2.29  x  10"3  m2  (4) 

Fig.  8  summarizes  these  fracture  energy  interpretations  of  the 
post-peak  experiments  which  are  consistent  with  the  8-10  fold 
increase  of  uniaxial  compressive  strength  ove'  the  tensile 
strength. 
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Direct  Tension  Test: 

\\ 

_L 

^  =  2.29  x  10- 3 

Uniaxial  Compression  Test: 


Direct  Shear  Test: 


U”  *  0.14  Nn 


v? 

=  -  10.! 


___  ♦  A  *  0. 
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Figure  8  Fracture  Energy  Release  Calculations 


c.  Direct  Shear  Experiments 

The  principal  results  of  these  experiments  are  illustrated  in  Figs. 
9,  10  and  11  which  depict  the  experimental  set  up  and  the  nominal 
shear  stress  versus  relative  lateral  displacement  response.  To  this  end 
prismatic  concrete  specimens  of  the  size  103  x  117  x  204  mm  were 
subjected  to  different  normal  confining  pressures  and  subsequently 
sheared.  The  experimental  observations  support  the  following 
conclusions: 

(i)  The  peak  strength  and  residual  strength  values  exhibit  a 
substantial  increase  with  increasing  confining  pressure  in 
accordance  with  the  concept  of  frictional  materials.  The 
residual  shear  strength  does  not  increase  proportionally  with 
the  confining  pressure  indicating  that  the  angle  of  internal 
friction  decreases  with  confinement. 

(ii)  Spontaneous  energy  release  under  low  confinement  diminishes  with 
increasing  confinement  level  in  favor  of  a  stable  post-peak 
response  throughout  the  shearing  process, 

(iii)  The  contained  failure  process  in  the  direct  shear  specimen  is 
initiated  by  distributed  tensile  cracking  well  before  reaching 
peak  strength.  This  initial  cracking  phase  is  followed  by 
compressive  strut-action  at  peak  and  subsequent  localization  of 
a  discrete  shear  band  due  to  mixed  mode  fracture  at  the  residual 
strength  level. 

(iv)  The  highly  non-homogeneous  failure  process  during  mixed  mode 
fracture  requires  inverse  identification  of  the  underlying 
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Direct  shear  apparatus.  Side  view  of  shear  device  and 
support  system. 


Nominal  shear  stress-relative  displacement  response 
curves  for  concrete  specimens  in  the  direct  shear 
apparatus. 
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material  parameters  in  order  to  interprete  the  response  behavior 
and  in  particular  the  "fragility"  of  the  specimen. 

(v)  The  traditional  measures  of  ductility  in  terms  of  the  ratio 
between  the  modulus  of  resilience  and  the  modulus  of  toughness 
is  of  little  meaning  particularly  since  continuous  shearing 
takes  place  at  the  residual  strength  level  in  form  of  relative 
rigid  body  translation  of  the  upper  versus  the  lower  portion, 
(vi)  The  fracture  energy  release  rate  for  developing  the  shear  band 
provides  a  more  appropriate  measure  of  ductility  which  remains 
fairly  insensitive  of  the  level  of  confining  normal  stress.  The 
area  under  the  descending  portion  of  the  lateral  force-lateral 
displacement  diagram  provides  a  first  approximation  of  the 
strain  energy  release  rate  for  Mode  II  fracture  when  the  energy 
dissipation  in  the  pre-peak  regime  is  discarded 

u  v 

s  r  r 

uf  =  /  pi  du!  +  /  pn  dun  3  '50  (5) 

UP  VP 

The  second  term  in  Eq.  5  represents  the  release  due  to  shear 
dilatanoy  in  the  post-peak  regime  which  contributes  less  than  20%  to  the 
energy  balance  in  the  case  of  a  normal  confining  pressure  of  on  * 

1.38  MPa.  Fig.  8  illustrates  the  interpretation  of  the  post-peak 
response  in  terms  of  fracture  energy  considerations  for  Mode  I  and  Mode 
II  failure.  Since  the  fracture  surface  at  the  residual  strength  level 


% 


corresponds  to  the  cross-sectional  area,  the  fracture  energy  release  rate 
for  mode  II  type  failure  follows  from  eq.  5 


.II 


aUj 
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2430  N/m  with  A$  =  2.06  x  lQ-^2  (6) 


which  is  almost  fifty  times  the  fracture  energy  release  rate  for  tensile 
cracking,  oj  »  60  N/m.  This  large  value  is  consistent  with  the  large 
increase  of  strain  energy  release  in  compression  and  lends  itself  to  a 
fracture  mechanics  interpretation  of  the  direct  shear  test  with  a  well 
defined  slip  surface  because  of  the  box  confinement.  Clearly,  in  view  of 
the  different  orders  of  magnitude  of  g|  and  G^l  concrete  failure 
will  always  be  preceded  by  tensile  cracking.  However,  the  post-peak 
behavior  in  direct  shear  is  controlled  by  the  separate  fracture  mechanism 
of  decohesion  which  is  responsible  for  the  formation  of  a  distinct  shear 
band  when  the  inclined  compressive  struts  are  sheared.  It  is  interesting 
that  the  "structural  aspects"  of  the  direct  shear  test  clearly  support 
two  distinct  fracture  mechanisms  in  spite  of  the  complexity  of  this 
experiment . 

(vii)  Interpretation  of  the  underlying  softening  mechanism  is 

complicated  by  the  interaction  of  debonding  during  tensile 
cracking  in  the  pre-peak  regime  and  the  subsequent  decohesion 
during  shearing  of  the  compressive  struts  in  the  post-peak 
regime.  Parametric  finite  element  studies  clearly  indicated 
that  these  mechanisms  are  distinctly  different  and  that  the 
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degradation  of  shear  strength  in  the  post-peak  regime  is 
primarily  controlled  by  decohesion  rather  than  debonding. 

(viii)  Ouring  shearing  the  confined  specimen  is  subjected  to 

substantial  dilatancy  which  results  in  considerable  vertical 
uplift  during  softening  in  the  post-peak  regime. 

(ix)  There  is  little  change  in  the  shear  strength  degradation  for 
nominal  confinement  levels  between  0.69  £  on  £  6.9  MPa.  In 
fact  mode  II  type  fracture  of  concrete  is  always  preceded  by 
tensile  cracking  and  thus  mobilizes  failure  below  the  transition 
point  for  brittle-ductile  fracture. 

3.2  Constitutive  Results 

The  principal  modelling  effort  was  directed  towards  the  development 
of  strain-softening  relations  for  tensile  cracking  and  frictional  slip 
with  the  understanding  that  strain-softening  does  not  constitute  a  local 
material  property.  The  fracture  energy  concept  of  Gf  *  const  was 
introduced  in  order  to  homogenize  the  surface  dominated  failure  mechanism 
within  an  elementary  volume  element  of  finite  dimension.  For  mode  I  type 
cracking  this  infers  that  the  fracture  energy  of  opening  a  discrete  crack 
of  surface  dAt  is  equivalent  to  the  energy  released  during  straining 
the  elementary  volume  dVe  in  tension  up  to  rupture 

Gf  A  «  /  ,  duf  dA  =  /  T  d*fdVe  (7) 

\  Ve 

This  fracture  energy  equivalence  furnishes  the  relationship  between  the 
crack  opening  displacement  duf  in  the  cracked  discontinuum  and  the 
tensile  fracture  strain  de^  of  the  homogenized  continuum 
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At/Ve  is  a  measure  of  the  crack  density  which  lends  itself  to  the 
interpretation  of  crack  spacing  ht=Ve/Aj..  Fig.  12  illustrates  the 
composite  fracture  model  for  uniaxial  tension.  In  this  case  the 
homogenization  of  a  single  tensile  crack  within  the  elementary  volume 
does  not  require  further  explanation,  it  remains  to  be  said  that  the 
generalization  to  two-  and  three  dimensions  is  fairly  straight  forward. 

In  fact  the  composite  fracture  theory  was  successfully  adopted  to  develop 
a  fracture  energy  based  strain-softening  plasticity  model  for  tensile 
cracking,  see  Appendix  B.  In  this  case  the  equivalent  tensile  fracture 
strain  Cf  was  utilized  as  control  variable  in  order  to  monitor 
strain-softening  of  the  fracture  energy  equivalent  continuum  in  terms  of 
the  following  relationship  between  the  softening  moduli  Ef-Ed 


E 


f 


(9) 


Note  that  the  strain-softening  modulus  Ef  is  directly  expressed  in 

d°t 

terms  of  the  fracture  modulus  E,  =  — -which  represents  the  post-peak 

a  dUf 

slope  in  terms  of  the  primary  test  variables,  the  nominal  strength  ot 
*  P/A  and  the  overall  crack  extension  Uf.  In  this  context  the 
characteristic  length  ht  =  Ve/At  converts  the  "crack-opening 
displacement"  into  equivalent  tensile  strain  of  the  homogenized 
continuum. 

Except  for  this  novel  definition  of  the  softening  modulus  Ef  the 
resulting  strain-softening  plasticity  model  follows  the  established 


25 


Fig.  1 2  Composite  Fracture  Hooel 
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concepts  of  plasticity  theory  and  leads  to  the  following  incremental 
stress-strain  relationship 

•  «-  •  .  .r  lE  m  m^iE  , 

d  =  ETt  where  lE-j.  =  lE -  (10) 

E^-Hn^E  m 

For  isotropic  degradation  of  the  tensile  strength  at  =  o(ef)  a  single 
softening  relationship  suffices  to  monitor  the  cracking  status  in  all 
directions  in  terms  of  the  equivalent  tensile  fracture  strain 

ef  *  /<e}>£  +  <i2>f  +  <e3>f  (ID 

In  view  of  the  inherent  directionality  due  to  tensile  cracking  this  is 
clearly  a  gross  simplification  which  is,  however,  conservative  since  the 
major  stress  component  oj  *  ot  controls  the  degradation  of  tensile 
strength.  Extension  to  anisotropic  softening  along  the  three  principal 
direction  does  not  introduce  conceptual  difficulties  but  requires 
eparate  monitoring  of  three  independent  strain-softening  processes. 
However,  the  usual  assumption  of  permanent  crack-memory  does  require 
special  previsions  for  "induced  anisotropy"  such  as  kinematic  rules, 
back-stress  meaures  etc.  in  order  to  accommodate  rotation  of  the  crack 
directions  during  rotation  of  the  principal  directions. 

The  "main  accomplishment"  of  the  constitutive  development  was  the 
generalization  of  the  composite  fracture  theory  from  Mode  I  type  cracking 
to  mixed  mode  fracture.  In  fact  compressive  splitting  as  well  as 
strain-softening  during  shear-faulting  under  triaxial  compression  could 
be  reconciled  within  the  format  of  the  fracture  energy  based 
strain-softening  plasticity  model  delineated  in  Appendix  B.  The  increase 
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of  strength  and  ductility  with  increasing  confinement  was  accomplished 
via  the  Leon  generalization  of  the  Mohr-Coulomb  criterion  plus  tension 
cut-off  and  the  adoption  of  in  Eq.  11  as  control  variable  for 
tensile  softening.  Thereby  the  tensile  strength  degrades  simultaneously 
with  the  internal  cohesion,  while  the  internal  friction  increases  in 
order  to  retain  a  residual  strength  envelope  which  passes  through  the 
transition  point  of  brittle-ductile  failure.  The  softening  modulus  Ef 
diminishes  to  zero  with  increasing  confinement  since  the  value  of  the 
characteristics  length  h  decreases  rapidly  in  compression 


(12) 


It  is  particularly  appealing  that  softening  in  compression  as  well  as 
tension  can  be  unified  within  the  concept  of  the  composite  fracture 
theory,  whereby  the  equivalent  tensile  fracture  strain  Cf  is  the  main 
control  variable  for  the  entire  failure  process. 

3.3  Computational  Results 

Recent  computational  Issues  have  focused  on  the  wel 1 -posedness  of 
quasistatic  and  dynamic  strain-softening  computations  and  the  underlying 
uniqueness  and  stability  of  the  numerical  results.  In  the  context  of  the 
research  on  "Finite  Elements  and  Localized  Failure"  these  issues  had  to 
be  addressed  because  of  the  structural  aspects  of  strain-softening 
computations.  In  other  words,  stable  overall  response  behavior  of  the 
structure  is  only  possible  if  the  strain-softening  regime  is  "contained" 
within  the  intact  elastic  portion  of  the  structure  and  if  the  solution  is 
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stabilized  by  displacement  or  general  arc-length  control.  Therefore, 
spatial  localization  of  the  material  instabilities  is  of  crucial 
importance  for  the  stable  control  of  the  overall  structural  response.  It 
is  this  strong  interaction  between  the  space  and  time  (load)  domain  which 
is  critical  for  the  survival  of  the  numerical  failure  simulation.  The 
papers  in  Appendices  C  and  0  summarize  the  fundamental  aspects  of 
quasistatic  and  dynamic  failure  analyses  with  the  aid  of  one  dimensional 
model  problems. 

Considering  the  incremental  energy  balance  between  the  intact 
elastic  storage  region  V-Vf  and  the  energy  release  region  Vf  the 
difficulties  of  progressive  failure  computations  become  apparen' . 
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The  finite  element  discretization  of  the  space  domain  demonstrates 
clearly  the  equivalence  between  the  crack  driving  force  in  fracture 
mechanics  of  discrete  crack  extension  and  the  equivalent  fracture  load 
due  to  progressive  strain-softening  in  the  homogenized  continuum. 


B  dofdV 
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and  do*  -  E  de* 


The  principal  difference  between  the  discrete  and  the  distributed 
fracture  approach  centers  around  the  role  of  the  crack  surface  area  dA, 
which  is  in  the  first  case  the  primary  control  variable  while  the 
strain-softening  region  Vf  is  only  a  secondary  result  of  distributed 
failure  analysis.  In  the  latter  case  the  geometric  extent  of  the  failure 
zone  is  only  defined  indirectly  in  terms  of  the  "loading  condition"  in 
the  displacement  (strain)  driven  conputational  algorithm 

TnSE  i  >  0  (14) 

This  condition  has  to  be  interpreted  for  load  increments  which  are  truly 
finite,  and  it  leads  to  the  concept  of  for  "iterative  reversibility"  in 
order  to  allow  for  elastic  unloading.  This  aspect  is  of  utmost 
importance,  if  elastic  snap-back  and  sharp  localization  of  the 
strain-softening  regime  are  to  be  captured. 

Appendix  D  describes  a  dynamic  model  problem  which  illustrates 
strain-softening  effects  and  localized  deformation  trapping  in  the  case 
of  wave  propagation  in  an  axial  force  member.  The  bar  was  subdivided 
into  9,  21  and  97  finite  elements  of  equal  length  and  subjected  to 
constant  end  velocities  c  =  +0.5  Cp  amounting  to  50%  of  the  acoustic 
wave  speed.  Fig.  13  describes  the  lay-out  of  the  model  problem  in  which 
the  predictions  of  the  traditional  strain-softening  approach  with  Ef  * 
const,  (independent  of  mesh  size)  were  compared  with  those  of  the 
composite  fracture  model  with  Gf  =  const,  (mesh  size  dependent  value  of 
the  strain-softening  modulus,  Ef  *  const.). 

The  resulting  stress  histories  are  depicted  in  Fig.  14  which  clearly 
demonstrates  the  strong  mesh-sensiti vity  of  the  traditional  Ef  = 
const,  approach.  In  fact  Fig.  15  illustrates  the  rapid  decrease  of  the 
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fracture  energy  which  diminishes  to  zero  with  decreasing  mesh  size,  i.e. 
Gf  +  0  as  h  ♦  0  when  Ef  *  const.  On  the  other  hand  the  composite 
fracture  model  maintains  fracture  energy  release  rates  which  remain 
constant  irrespective  of  the  mesh  size.  Fig.  16  illustrates  the  time 
history  of  strain  localization  at  the  center  of  the  tension  specimen  for 
the  three  idealizations.  Clearly,  it  is  primarily  a  question  of  spatial 
resolution  to  capture  tensile  cracking  by  strain-softening  of  an 
equivalent  continuum.  It  is  rather  intriguing  that  both  strain-softening 
models  reproduce  the  actual  separation  phenomenon  in  the  form  of 
deformation  trapping  in  the  center  element,  whereby  the  conventional 
strain-softening  approach  with  Ef  =  const,  leads  to  a  more  pronounced 
strain  localization  according  to  the  analytic  solution  of  strain 
singularity  in  a  domain  of  measure  zero,  i.e.  e  ♦  »  as  h  +0.  Obviously, 
both  continuum  models  cannot  capture  the  actual  separation  and 
concomitant  displacement  discontinuity  of  the  center  of  the  bar,  however, 
both  strain-softening  formulations  simulate  the  fracture  process  in  the 
form  of  deformation  trapping,  whereby  the  composite  fracture  model 
maintains  a  constant  fracture  energy  release  rate  irrespective  of  the 


CONSTANT  NRWWO  *U>H 


Fig.  14  Global  Response  Behaviour 


Fig.  15  Mesh  Sensitivity  of  Fracture  Energy 
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Fig. 16  Localization  of  Tensile  Failure 
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In  the  course  of  the  research  on  "Finite  Elements  and  Localized 
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published  in  technical  conference  proceedings  and  technical  Journals. 
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Computational  Aspects  of  Strain-Softening  and  Localization  in 
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Computational  Aspects,  New  Orleans,  Dec.  10-14,  1984,  ASME  Symp. 
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should  complete  nis  u : ssercat l on  witnin  another  1.5  years. 


6.  LIST  OF  ORAL  PRESENTATIONS 


In  view  of  the  strong  interest  on  concrete  failure  the  principal 
investigators  participated  very  actively  in  workshops,  conferences  and 
international  symposia  in  order  to  present  early  results  of  the  Air  Force 
project  on  "Finite  Elements  and  Localized  Failure".  In  fact,  one  of  the 
principal  investigators  organized  a  Symposium  on  "Constitutive  Equations: 
Macro  and  Computational  Aspects"  during  the  ASME  Winter  Annual  Meeting 
1984,  New  Orleans  Dec.  10-14,  1984  with  17  papers  and  a  set  of 
substantial  Conference  Proceedings  of  272  pages.  The  chronological  list 
below  summarizes  the  presentations  of  technical  meetings  and  university 
seminars  during  the  funding  period  1982-1985. 


6.1  International  Conferences 

1.  "Numerical  Modeling  of  Strength  and  Deformation  Behavior  in  the 
Direct  Shear  Test",  Inti.  Conf.  Constitutive  Laws  for  Engineering 
Materials,  Tucson,  Arizona,  Jan  10-14,  1983. 

2.  "Finite  Element  Analysis  of  Concrete  Failure  in  Shear",  Air  Force 
Symp.  Interaction  of  Non-Nuclear  Munitions  with  Structures,  U.S.  Air 
Force  Academy,  Colorado  Springs,  May  10-13,  1983. 

3.  "Experimental  Modeling  of  Strength  and  Deformation  Behavior  of 
Concrete  in  Direct  Shear",  Air  Force  Symp.  Interaction  of 
Non-Nuclear  Munitions  with  Structures,  U.S.  Air  Force  Academy, 
Colorado  Spring,  May  10-13,  1983. 

4.  "Finite  Element  Modeling  of  Concrete  Failure  in  Shear",  4th  ASCE  - 
Engineering  Mechanics  Specialty  Conference,  Purdue  University, 
Lafayette,  May  23-25,  1983. 

5.  "Computational  Stability  and  Uniqueness  of  Strain-Softening 
Materials",  W.  Prager  Symposium  on  Mechanics  of  Geomaterials, 
Northwestern  University,  Evanston,  Sept.  11-15,  1983. 

6.  "Strain-Softening,  a  Material  Property,  that  is  the  Question?"  ASCE 
Spring  Convention,  Atlanta,  Georgia,  May  14-18,  1984. 
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7.  "Composite  Damage  Model  for  Localized  Fracture",  5th  ASCE  - 
Engineering  Mechanics  Specialty  Conference,  University  of  Wyoming, 
Laramie,  Aug.  1-3,  1984. 

8.  "Continuous  Fracture  Computations  and  Strain  Localization  in 
Cementitious  Solids  Subjected  to  Tension  and  Shear",  XVIth  IUTAM 
Congress  of  Theoretical  and  Applied  Mechanics,  Lyngby,  Denmark, 

Aug.  19-25,  1984. 

9.  "Experimental  and  Computational  Aspects  of  Concrete  Fracture,  Inti. 
Conf.  Computer  Aided  Analysis  and  Design  of  Concrete  Structures", 
Split,  Yugoslavia,  Sept.  17-21,  1984. 

10.  "Constitutive  and  Computational  Aspects  of  Strain-Softening  and 
Localization  in  Solids",  ASME-WAM'84  Symp.  Constitutive  Equations: 
Micro,  Macro  and  Computational  Aspects,  New  Orleans,  Dec.  10-14, 

1984. 

11.  "A  Composite  Fracture  Model  for  Localized  Failure  in  Cementitious 
Materials",  2nd  Air  Force  Symp.  Interaction  of  Non-Nuclear  Munitions 
with  Structures,  Panama  City  Beach,  Florida,  April  15-19,  1985. 

12.  "Experimental,  Constitutive  and  Computational  Aspects  of  Concrete 
Failure",  U.S. -Japan  Symposium  on  Finite  Element  Analysis  of 
Reinforced  Concrete  Structures,  Tokyo,  May  21-24,  1985. 

13.  "Computational  Modeling  of  Material  Softening",  2nd  Joint  ASCE/ASME 
Mechanics  Conference,  Albuquerque,  June  23-26,  1985. 

14.  “Stability  and  Uniqueness  of  Strain-Softening  Computations", 
Europe-US  Symposium  Finite  Element  Methods  for  Nonlinear  Problems, 
NIT  Trondheim,  Norway,  Aug.  12-16,  1985. 

15.  "Computational  Modeling  of  Distributed  and  Localized  Fracture",  8tn 
Inti.  Conf.  Structural  Mechanics  in  Reactor  Technology,  Aug.  19-23, 

1985. 

16.  "Composite  Fracture  Model  for  Strain-Softening  Computations  of 
Concrete",  Int.  Conf.  Fracture  Mechanics  of  Concrete-Fundamentals 
and  Applications,  EPF  Lausanne,  Switzerland,  Oct.  1-3,  1985. 
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6.2  Workshops  and  University  Seminars 

1.  "Finite  Element  Modelling  of  Concrete  Failure  in  Shear",  Graduate 
Seminar  at  Northwestern  University,  Evanston,  May  20,  1983. 

2.  "Computational  Aspects  of  Concrete  Mechanics,"  US-Dutch  Cooperative 
Research  Symposium,  Delft  University  of  Technology,  Delft,  June 
22-24,  1983. 

3.  "Nlchtl ineare  Stof fmechani k  und  Finite  Elemente"  Graduate  Seminar, 
University  of  Innsbruck,  June  29,  1983. 
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4.  "Finite  Elements  and  Localized  Failure",  US-AFWL  Material  Properties 
Meeting,  Albuquerque,  Oct.  3,  1983. 

5.  "Strain-Softening,  Experimental  and  Computational  Aspects",  U.S. 

Army  Waterways  Experiment  Station,  Vicksburg,  Oan.  17,  1984. 

6.  "Experimental  Results  on  Concrete  Failure",  Workshop  on  Inelastic 
Deformations  and  Failure  Modes,  Northwestern  University,  Evanston 
Nov.  18-21,  1984. 

7.  "Experimental  and  Computational  Aspects  of  Strain-Softening",  Hardy 
Cross  Lecture,  University  of  II  line  is  at  Urbana-Champaign,  April  19, 
1985. 

8.  "Strain-Softening  Computations  r,'  Brittle  Failure  Processes", 
Graduate  Seminar  at  Lund  Institute  of  Technology,  Lund,  Sweden, 

Aug.  16,  1985. 

9.  "Experimental,  Constitutive  and  Computational  Aspects  of  Concrete 
Failure,"  Seminar  at  Applied  Mechanics  Division,  Sandia  National 
Laboratories,  Albuquerque,  Nov.  8,  1985. 
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7.  NEW  FINDINGS 


The  research  effort  on  "Finite  Elements  and  Localized  Failure"  led 
to  a  better  understanding  of  brittle-ductile  failure  phenomena  in 
particulate  composites  such  as  concrete.  Most  of  the  findings  were 
listed  in  Section  3  above  under  the  heading  of  significant 
accomplishments  on  the  experimental,  the  constitutive  and  computational 
fronts. 

The  most  important  result  of  the  investigation  on  concrete  fragility 
is  that  concrete  failure  is  a  fracture  driven  process  in  which  the 

1 

nominal  tensile  strain  plays  an  instrumental  role  to  homogenize  the 

effect  of  coalescing  microdefects  into  a  single  macrodefect.  In  fact  I 

| 

strain-softening  is  not  a  local  material  phenomenon  but  rather  the  j 

i 

consequence  of  distributing  the  degradation  effects  of  microdefects  in  an  j 

elementary  but  finite  volume.  Therefore,  the  principal  issue  of  j 

extracting  useful  fragility  information  from  servo-controlled  post-peak 
experiments  on  laboratory  specimens  centers  around  the  homogenization  of 
the  fracture  process  in  the  discontinuum  by  an  equivalent 
strain-softening  continuum. 

To  this  end  a  series  of  servo-controi led  experiments  were  performed 
on  cylindrical  and  prismatic  concrete  specimen  . 1  'rovided  unique 
data  of  concrete  strength  and  deformation  behavio  ell  as  fragility 
when  subjected  to  direct  tension,  triaxial  compression  and  direct  shear. 

The  experiments  furnished  quantitative  evidence  on  the  brittle-ductile 
transition  of  ylindrical  failure  which  diffuses  entirely  under 
sufficient  confinement.  The  different  tests,  also  substantiated  two 
different  failure  modes,  (i)  a  planar  failure  mode  under  direct  tension 

l 

I 

i 

i 
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and  direct  shear  due  to  progressive  debonding  and  decohesion,  and  (ii)  a 
cylindrical  failure  mode  under  trlaxial  compression  due  to  simultaneous 
decohesion  and  Increase  of  internal  friction.  The  two  different  failure 
modes  support  the  concept  of  tensile  strain  as  thi  control  variable  for 
the  degradation  of  strength,  whereby  planar  failui^  is  caused  by  the 
single  value  of  major  principal  strain,  while  axisymmetric  or  cylindrical 
failure  is  the  result  of  double  valued  major  principal  strains. 

For  interpretation  of  these  failure  phenomena  the  "Composite 
Fracture  Theory"  was  developed  which  resorts  to  the  fracture  energy 
release  rate  in  order  to  monitor  the  fracture  process  in  an  elementary 
volume  of  finite  size  by  “equivalent"  strain-softening.  In  this  way  the 
fracture  mechanics  of  tortuous  cracking  In  brittle  solids  is  replaced  by 
fracture  energy  equivalent  strain-softening.  In  fact  the  traditional 
format  of  strain-softening  plasticity  could  be  readily  adopted,  except 
for  the  new  definition  of  the  strain-softening  modulus  Ef  in  terms  of 
the  fracture  energy  release-rate  and  the  crack  density  per  unit 
volume 

dot  dc^  duf  Ve 

f  d cf  duf  def  d  \ 

d2°f  i 
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The  resulting  tangential  material  law 
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maintains  constant  fracture  energy  release  rate  when  the  elementary 
volume  is  subjected  to  fracture  up  to  rupture  where  the  stress  reaches 
the  residual  strength  level.  Although  the  composite  fracture  model  was 
initially  developed  for  tensile  cracking  it  was  further  extended  to 
include  the  prevalent  strength  and  failure  phenomena  in  compression,  such 
as  “splitting  compression"  and  "shear  faulting".  Beyond  the 
brittle-ductile  transition  point  the  strain-softening  formulation  reverts 
to  strain-hardening  during  which  volume  driven  damage  accumulation 
replaces  the  surface  driven  fracture  process.  As  a  result  the  composite 
fracture  model  covers  the  entire  spectrum  of  failure  phenomena  of 
cmcrete,  starting  from  brittle  tensile  cracking  under  Mode  I  type 
fracture,  to  Increasingly  ductile  fracture  up  to  continuous  hardening 
beyond  the  transition  point  under  sufficiently  high  confinement. 

The  ensuing  strain-softening  computations  of  fracture  energy 
equivalent  elastic  plastic  solids  lead  to  stringent  prerequisites  of  the 
numerical  simulation  procedure.  The  spatial  discretization  with  finite 
elements  is  limited  by  mesh  size  restrictions  for  local  stability  and  the 
numerical  integration  has  to  make  proper  provisions  for  elastic  unloading 
in  order  to  allow  for  snap-back  and  concomitant  strain  localization  in 
load  increments  which  are  truly  finite.  Contained  fracture  calculations 
require  numerical  stabilization  with  displacement  and/or  arc  length 
control  with  computational  strategies  which  are  sufficiently  robust  to 
survive  local  instabilities  and  which  are  on  the  other  hand  sufficiently 
sensitive  to  detect  failure  of  the  overall  structure. 


8.  CONCLUDING  REMARKS  AND  FUTURE  WORK 

Considerable  progress  was  made  in  improving  the  understanding  of 
concrete  fragility.  Novel  post-peak  experiments  were  performed  in 
unnotched  concrete  specimens  which  demonstrated  the  brittle-ductile 
transition  of  failure.  In  addition  the  composite  fracture  theory  was 
developed  in  order  to  monitor  the  fracture  energy  release  rate  within  the 
framework  of  strain-softening  plasticity.  The  computational  issues  of 
contained  strain-softening  computations  were  addressed  in  quasistatic  and 
dynamic  environments. 

Several  Issues  remain,  however,  unresolved  and  deserve  further 
attention: 

(i)  For  a  better  understanding  of  the  actual  failure  process  in 

particulate  composites  such  as  concrete  "mesomechanical"  studies 
are  needed  on  the  experimental  as  well  as  computational  level. 

In  this  case  concrete  could  be  resolved  into  its  constituents, 
the  aggregate,  the  hardened  cement  paste  and  the  interface  In 
order  to  examine  the  mesomechanical  aspects  of  damage  and 
fracture  in  particulate  composites.  In  view  of  the  experimental 
difficulties  with  SEM,  CAT-SCAN  X-ray,  hypersonic  and  acoustic 
emission  techniques  a  probabilistic  numerical  simulation  would 
provide  valuable  information  on  the  internal  microprocesses  and 
the  associated  macrophenomena  in  conjunction  with  a  substantial 
experimental  data  base. 

(’1)  The  "anisotropy"  or  rather  the  load  history  effects  of  failure 
should  be  studied  in  more  detail  in  view  of  the  strong 
directionality  of  damage  and  fracture  and  the  large  difference 
in  ductility  between  planar  and  cylindrical  failure.  To  this 
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end  servo-controlled  laboratory  experiments  could  provide  very 
valuable  information  for  loading  histories  involving  large 
rotations  of  the  loading  direction. 

(iii)  The  stabilization  of  contained  strain-softening  computations 
deserves  immediate  attention.  The  development  of  robust 
algorithms  for  post-critical  response  calculations  is  of  great 
practical  importance  for  reliable  fragility  analyses  of  actual 
structures. 
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ABSTRACT 

The  current  study  examines  the  degradation  ot  strength  In  cementitious  mate¬ 
rials  loaded  In  tension  and  shear.  Based  on  experimental  observations  a  straln- 
soltcnlng  model  la  developed  which  encompasses  both  continuous  decohesion  due  to 
distributed  microdefects  as  well  as  the  fracture  energy  release  rste  concept  for 
single  macrodefecti.  Distributed  and  localized  fracture  are  combined  in  the  compo¬ 
site  damage  model  of  an  equivalent  homogenized  continuum  of  finite  dimension.  In 
one  case  the  modal  reduces  to  che  ctsattionmi  strain-softening  formulation  of 
continuous  damage  mechanics.  In  the  other  esse  It  furnishes  a  non-local  softening 
law  Which  depends  on  the  characteristic  fracture  dimension  as  well  as  the  fracture 
configuration. 

The  performance  of  the  composite  damage  model  Is  Illustrated  with  esperl- 
eent.l  data  of  the  direct  tension  and  trlaxlal  compression  tests  on  NX-slze 
mortar  specimens.  Posc-crltlcsl  response  predictions  illustrate  that  the  frac- 
turt  energy  release  rate  should  reaein  Invariant  If  localized  tensile  cracking  Is 
to  be  captured.  In  contrsdlstlnctlon,  the  shear-softening  law  should  be  a  constant 
fracture  property  If  distributed  fracture  le  to  be  described  In  the  form  of 
frictions’,  slip- 

These  basic  concepts  for  distributed  and  localized  fracture  are  subsequently 
utilized  to  formulate  a  comprehensive  crlsxlal  fracture  model  for  both  tensile 
debonding  and  frictional  decohesion  In  cementlclous  solids.  These  two  mechanisms 
are  described  by  a  non-aasoc  iatlve  Mohr-Coulomb  formulation  with  tension  cut-olt 
in  unlch  the  tensile  strain-softening  is  controlled  by  the  critical  fracture  ener¬ 
gy  release  rate.  The  performance  of  the  overall  model  Is  verified  by  experimental 
results  obtained  In  the  direct  shesr  test  in  which  both  censile  cracking  as  well 
as  frictional  slip  take  place  Interactively  similar  to  contained  fracture  process¬ 
es  In  highly  Indetermlnace  structures. 

1.  INTRODUCTION 

In  the  recent  past  a  vigorous  debate  has  been  concerned  with  the  postcrltl- 
cal  or  post-peak  strength  behavior  of  engineering  materials.  In  particular  the 
limited  ductility  of  cementitious  otcerlals  has  been  of  considerable  Interests  ad¬ 
vocates  of  continuous  strain-softening  models  have  been  in  open  conflict  with  the 
(racture  mechanics  proponents  on  this  Issue.  Clearly,  che  computational  advanta¬ 
ges  of  che  smeared  approach  provide  the  predominant  argument  [n  favor  of  the  con¬ 
tinuous  damage  models  in  wnlch  progressive  decohesion  due  to  distributed  slcrode- 
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fects  Is  described  by  an  equivalent  strain-softening  behavior  of  a  homogenized 
conrlnuua.  On  the  other  hand,  serious  questions  arise  with  regard  to  the  validity 
of  computational  predictions  la  the  case  of  fracture  localization  when  the  struc¬ 
tural  Integrity  Is  determined  by  the  sev.  rlty  of  a  single  macrcdefect.  This  type 
of  failure  has  been  traditionally  the  realm  of  fracture  mechanics  In  which  the 
initiation  and  propagation  of  a  single  discontinuity  Is  considered  the  determining 
factor. 

The  main  objective  of  the  present  investigation  Is  to  compare  the  conflict¬ 
ing  strategies  and  to  reconcile  the  fundamental  differences  of  strain-softening 
and  fracture  mechanics.  To  this  end,  a  composite  damage  formulation  will  be  deve¬ 
loped  which  embodies  the  basic  concepts  of  single  and  distributed  fracture  e- 
vents.  On  one  end  of  the  spectrum  the  composite  damage  model  enforces  a  fracture 
energy  release  rate  which  la  constant  in  the  case  of  a  highly  localized  discrete 
discontinuity.  On  the  other  end  of  the  spectrum  It  reduces  to  the  continuous  dam¬ 
age  theory  for  distributed  fracture.  The  constitutive  formulation  follows  strain- 
softening  plasticity  theories  In  which  the  non-assoclated  flow  rule  Is  recast  Into 
a  non-local  fracture  rule  which  Includes  the  characteristic  fracture  dimension  of 
the  fracture  process  zone.  The  underlying  composite  damage  parameters  for  tension 
and  shear  are  calibrated  from  post-critical  experimental  data  obtained  from  the 
direct  teoslon  and  trlaxial  compression  rests.  These  two  independent  sets  of  ser¬ 
vo-control  Jed  experiments  are  subsequently  utilized  to  study  the  performance  of 
the  computational  post-critical  predictions  with  regard  to  their  ability  to  cap¬ 
ture  distributed  microdefects  as  well  as  single  macrodef ects .  These  two  opposing 
concepts  arc  implemented  within  the  unified  framework  of  the  composite  damage 
theory  which  decomposes  an  elementary  volume  Into  a  series  of  Intact  elastic  zones 
and/or  a  single  daaage  zone  In  which  oriented  tensile  cracking  and  frictional  slip 
take  place  In  a  highly  localized  manner.  The  Implications  on  the  computational 
predictions  of  the  post-critical  behavior  are  then  llluscrated  with  fracture 
Investigations  of  the  direct  tension  test  and  the  trlaxial  compression  test  on 
cylindrical  mortar  specimens  for  which  experimental  data  are  available  for  compar¬ 
ison.  The  uniqueness  of  the  post-cr 1 t leal  response  predictions  Is  examined  with 
the  aid  of  systematic  mesh  sensitivity  studies  of  unnotched  test  specimens  which 
exhlhlc  pronounced  softening  behavior  both  in  tension  and  shear. 

The  disjoint  formats  of  the  composite  damage  model  are  subsequently  Integra¬ 
ted  Into  a  crlaxlal  straln-sof teolng  model  of  non-assoclated  plasticity  for  com¬ 
bined  fracture  In  tension  and  shear.  The  computational  aspects  of  the  strain 
driven  algorithm  are  examined  with  regard  to  the  finite  Increments  which  are  used 
In  computational  Implementation  and  analysis.  To  this  end  an  Implicit  backward 
scracegy  is  adopted  which  assures  that  the  underlying  non-assoclated  fra'cture  rule 
satisfies  the  constitutive  constraints  both  at  the  beginning  and  at  the  end  of  the 
finite  tloe  Interval.  Particular  attention  Is  paid  to  Che  combined  effect  of 
softening  In  tension  and  in  shear. 

Last,  the  composite  daaage  model  is  applied  to  scudy  the  postcrltlcal 
response  behavior  of  prismatic  mortar  specimens  In  the  direct  shear  test.  In  this 
case  contained  fracture  Cakes  place  both  In  the  form  of  tensile  cracking  as  well 
as  frictional  slip  depending  on  the  ratio  between  normal  and  tangential  load 
levels.  In  this  sense  Che  direct  shear  teat  exhibits  all  Ingredients  of  a  highly 
Indeterminate  structure  In  which  progressive  propagation  of  distributed 
olcrudetects  lead*  to  an  overall  degradation  of  stiffness  and  strength  before 
coalescence  Into  a  single  macro-defect  cakes  place  which  Is  ultimately  reeponslble 
for  the  structural  failure. 

2.  DIRECT  TENSION  TEST 

For  simplicity  let  us  first  consider  the  direct  tension  test  In  order  to  Il¬ 
lustrate  the  underlying  fracture  mechanics  Issues  and  streln-sof tenlng  concepts. 
Fig.  1  depicts  the  basic  experiment  layout  and  the  response  curve  for  linear- 
elastic,  perfectly  brittle  behavior.  The  critical  fracture  energy  release  rate 
C;  follows  the  cnecgy  balance  relation  vnlch  reduces  for  prescribed  displacement 
boundary  conditions  to  the  simple  expression 
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where  n  denotec  the  total  potential  energy,  of  the  Ideal  strength  of  a  flawless 
crystalline  material  and  dA  the  change  of  crack  surface  (1).  In  the  spirit  of 
the  Griffith  concept  of  brittle  fracture  Gc  corresponds  to  the  surface  energy 
£  density  y,  while  the  characteristic  crack  dimension  t  defines  the  range  of  the 

Interatomic  bond  forces  In  the  order  of  the  lattice  spacing  for  Ideally  cohesive 
and  crystalline  materials. 
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Frg.l  Direct  tension  test  of  Ideally  brittle  material 


In  the  case  of  cementitious  materials  the  question  arises  if  the  fracture  en¬ 
ergy  release  rate  C.  Is  a  proper  material  property  which  characterizes  the  ten¬ 
sile  cracking  Independently  of  the  particular  geometry  and  boundary  conditions. 

In  fact,  this  Issue  was  raised  p-evlouslv  In  the  contest  of  fragmentation  of  r^cks 
by  Paul  von  Rlttlnger  In  1867  ( 2 1  who  stated  that  the  newly  created  surface  area 
)  la  directly  propoi clonal  to  the  energy  expended  during  fracture,  l.e.  C,  »  cons¬ 

tant.  In  the  recent  past,  chlx  fundamental  fracture  mechanics  concept  was  adopted 
by  Hlllerborg  et  al  (3)  In  formulating  the  fictitious  crack  band  model  for,  con¬ 
crete.  This  discrete  fracture  strategy  was  extended  by  Bazant  et  al  (4|  to  formu¬ 
late  the  crack  band  theory  for  distributed  fracture  within  a  process  zone  of  fin¬ 
ite  wtdch.  Along  similar  lines  the  experimental  and  computational  fracture  stud¬ 
ies  by  Ingraffea  and  Saouma  (5)  were  based  on  the  equivalent  fracture  postulate  of 
a  critical  stress  Intensity  factor  to  constitute  an  Invariant  fracture  property. 
t  In  view  of  this  fundamental  Issue  let  us  examine  whether  the  fracture 

mechanics  hypothesis  applies  to  tensile  cracking  of  cementitious  materials  and  how 
It  can  be  Incorporated  Into  the  strain-softening  formulation  of  an  equivalent  con¬ 
tinuum.  To  this  end  let  us  consider  chc  actual  experimental  data  of  the  direct 
tension  tesc  shown  In  Fig.  2.  We  recall  that  che  cylindrical  mortar  specimen  Is 
tested  In  a  very  stiff,  servo-controlled  MTS-apparatus  which  maintains  stable 
displacement  control  In  the  post-critical  response  regime.  In  fact,  no  noticeable 
change  of  the  geometry  occurs  during  the  entire  cost  except  for  a  fine  cleavage 
^  crack  which  becomes  visible  only  at  the  residual  load  level.  At  this  stage  the 

unnotched  mortar  specimen  separates  Into  two  competent  parts,  while  the  crack 
opening  displacement  localizes  within  the  single  crack  band  between  them. 

In  the  absence  of  noticeable  geometric  distortions  before  separation  of 
the  specimen  Into  two  parts,  che  elementary  uniaxial  Interpretation  of  the  direct 
tension  data  leads  to  chc  following  setffness  and  scrength  properdes  In  terns  of 
nominal  stress  and  scraln  measures: 

i  ^  f  ^  f  '*f 

F.  •  —  -  6.895  MPa ;  a,  •  —  •  1.59  MPa;  c  •  —  -  2.30  x  lO-1*  (2) 

rf  l  a  f  L 
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Fig.  2  Direct  tension  of  HX-glze  specimen 


Assuming  chat  Che  actual  strain  localization  during  tensile  cracking  can  be 
replaced  by  an  equivalent  uniform  distribution  of  olcrodefects  over  the  gauge 
length  L  we  can  extract  the  strain-softening  modulus 

E  -  0.82  £  -  5.653  MPa  (3) 

s 

which  remains  constanc  In  the  case  of  a  bilinear  stress-strain  idealization  of 
the  pre-and  post-critical  response  behavior.  Clearly  If  E,  would  be  a  proper' 
material  property  It  should  have  the  sane  value  for  each  material  particle  within 
the  gauge  length.  Obviously  this  Is  not  the  case  since  Es  is  based  on  a  nominal 
strain  measure  equal  Co  the  overall  axial  displacement  divided  by  the  specimen 
height  which  contains  a  single  tensile  crack.  In  other  words  the  strain-softening 
modulus  F.,  is  a  structural  property  Intimately  tied  to  the  particular  tejt 
configuration  and  fraccure  condition.  Therefore,  It  cannot  be  extended  right  away 
to  different  geometries  and  boundary  conditions  wlcnout  a  rational  concept  benlnd 
this  extrapolation.  In  fact,  the  equivalent  straln-sof cenlng  model  Infers  that 
distributed  cracking  cakes  place  In  a  truly  finite  volume  of  macerlal  In  which 
non-local  theories  for  continuous  media  with  microstructure  simulate  progressive 
fracture  (61.  (7). 

On  the  other  hand  the  displacement-controlled  formation  of  a  single  crack 
band  In  the  direct  tension  specimen  can  be  used  to  determine  the  critical 
fracture  energy  release  rate  In  tension  Cc  c.  According  to  Petersson  (8|  the 
Integral  of  nominal  normal  stress-displacement  curve  shown  In  Fig.  2  yields  a 
first  order  approximation  of  the  critical  fracture  energy  release  race  In  tension 
for  the  particular  mortar  undei  consideration 
u 

r 

C  •  /  0  du  *  44  N/b  (4) 

c .  t  1  1  i 

o 

This  value  agrees  falrlv  well  with  other  fracture  mechanics  data  reported  for 
concrete  and  mortar  (8|.  The  characteristic  fracture  dimension  Is  In  this  case 

0  E 

1  ■  -1^2 —  -  0.12  o  (5) 

‘  °f 


-  ‘4  — 


which  Is  definitely  several  orders  of  magnitude  beyond  the  lattice  spacing  of 
crystalline  solids.  In  fact  Xt  Is  of  the  same  order  of  magnitude  as  the 
entire  NX-slte  specimen  which  gives  rise  to  the  question  if  the  fracture  energy 
release  rate  Is  a  'proper'  fracture  property.  tf  wo  adopt  the  notion  of  a 
fracture  process  rone  with  degrading  material  properties  characterized  by  £4, 
then  an  elementary  series  model  of  the  uniaxial  tension  specimen  [9)  yields  the 
following  value  for  the  width  dc  of  the  characteristic  fracture  process  zone 


This  result  Indicates  that  the  actual  width  of  the  fracture  process  zone  dt  may 
be  considerably  smaller  chan  Xe,  depending  on  the  ratio  between  the  tangential 
material  stiffness  in  fhe  damage  zone  and  the  elastic  properties  of  the  undamaged 
zone.  However,  lc  appears  that  the  width  of  the  fraccure  band  d(  will  be  still 
sufficiently  large  such  chat  the  validity  of  the  underlying  fracture  mechanics 
concept  must  be  questioned. 


2 . 1  Composite  Damage  Model  for  Tensile  Cracking 


In  order  to  compare  Che  compucatlonal  response  predictions  resulting  from 
che  constanc  strain-softening  foraulaelon  as  opposed  to  the  constant  fracture 
energy  release  race  strategy,  we  resort  Co  the  composite  damage  model  Illustrated 
In  Fig.  3  (.See  also  refs.  19  —  11)). 


In  the  case  of  tensile  cracielng  me  elementary  volume  V  is  composed  of  the 
undamaged  elastic  zone  and  Che  localized  crack  band  V,  of  width  1,  and 
longm  h.  in  which  oriented  a  1  c roc racs log  is  confined.  •'^sorting  to  the  maximum 
stro-.s  criterion  fur  cracking,  the  orientation  01  me  tensile  crack  band  is 
defined  by  the  angle  a  between  the  reference  axis  and  the  minor  principal  stress 
Jlrcctlcn.  Figure  1  depicts  a  typical  rectangular  geoaecry  ot  the  elementarv 
volume  of  thickness  c 


V 


where 


V  *  4  b  c 


(7) 


i 

( 

4 
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end 


2  dchcc 


In  this  case  2d:  denotes  the  width  of  the  tensile  crack  band  and  lL  the  Length 
within  the  elementary  volume,  where 


h  -  — 

t  sin  a 


(8) 


whichever  Is  smaller. 


For  straining  up  to  rupture  the  energy  release  of  th'1  .opposite  damage 
element  Is  confined  to  the  localized  crack  band  because  of  rull  unloading  of  the 
elastic  zone.  In  che  simple  case  of  linear  degradation  of  the  tensile  strength 
o^<  the  tensile  energy  release  of  the  composite  element  follows  from 


‘2f 

^  d‘  ^  C 


(9) 


The  critical  fracture  energy  release  rate  for  opening  the  crack  band  with  the 
surface  *At  *  htc  results  In 


where 


(10) 


The  energy  dissipation  characteristics  of  the  equivalent  continuum  element  in 
which  che  crack  band  Is  distributed  uniformly  are  defined  by 


s 


From  the  equivalence  of  the  energy  release  during  cracking  ailf  *  the 

sofcenlng  modulus  E,  In  tension  depends  on  che  fracture  energy  release  rate 
Cct  or  the  alternative  characteristic  fracture  dimension  x t  as  well  as  the 
crack  geometry 


M2) 


This  expression  clearly  demonstrates  that  the  straln-sof  Cenlr.g  modulus  E,  of  the 
equivalent  conclnuun  element  depends  also  on  the  geometry  of  the  elementary  volume 
and  Chus  exhibits  a  pronounced  size  effect.  In  other  words,  E,  varies  with  the 
mesh  size  of  the  particular  finite  element  layout  If  che  fracture  energy  release 
rate  Gcs  and  the  respective  characteristic  fracture  dimension  constitute 
Invariant  fracture  properties.  As  a  result  the  fracture  mechanics  concept  leads 
to  a  non-iccal  foroac  of  the  equivalent  softening  relation  which  Is  fundamentally 
different  from  the  local  c&nscltutlve  formats  of  strain-softening  plasticity  and 
continuous  damage  mechanics  (12,131. 

Before  we  turn  Co  the  computational  study  lec  us  consider  the  material 
Identification  of  che  composite  damage  model.  Contrary  to  previous  proposals  of 
Bazanc  et  il  [ 4 |  the  equivalent  softening  modulus  In  Eq .  12  Is  completely  defined 
by  two  parameters,  the  Initial  tensile  strength  of  and  the  fracture  energy 
release  rate  Cc . ;  or  the  alternative  fracture  dimension  lt.  In  this  case 
there  Is  no  need  to  speculate  on  che  width  of  the  blunt  ersex  band  dt  or  the 
material  benavior  £4  in  the  localized  damage  zone.  In  fact  the  equivalent 
strain-softening  expression  reduces  to  the  following  relatlonsnlp  in  the  case  of 
the  direct  tension  test  of  Section  2.1. 


where  l  *  0.12  m 


and  h  -  a  for  o-0  (13) 

t 


Recall  that  b*-  denotes  the  nelght  of  the  elementary  volume 
transverse  to  the  crack  band.  Clear.;/  the  maximum  size  of 


and  defines  the 
the  mesh  layout 


s  1  ze 
l  s  re- 


L 


i 
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strlctcd  to  b*<  1  •  0.2*  a  which  results  in  perfectly  brittle  sottenlng  of 

the  equivalent  rontinutn  element  with  £<-  •  .  *>n  c he  other  hand  a  decrease  ot 

mem  sire  would  lend  to  a  corresponding  Increase  of  ductility.  In  fact  sirtaally 
Ideal  l  v  plastic  behavior  with  -  S/IO0  is  obtained  if  the  height  of  the 
elementary  volume  reduces  to  *  1.*  era. 

0a  the  other  hand  we  may  asu  oui  selves  under  what  condition  the  constant 
stra l teoing  approach  and  the  constant  fracture  energy  ip  Lease  rate  concept 
coincide.  In  Che  case  ot  our  particular  tension  test  this  corresponds  to  Che 
question  for  vh.t  sl2e  of  b^  do  we  recover  the  softening  modal  us  in  Eq .  3. 


(U) 


The  characteristic  dimension  of  the  fracture  process  zone  Xt  •  0.12  a  and  the 
height  bK  •  0.108  m  of  the  elementary  voluae  will  reproduce  exactly  the  constant 
softening  modulus.  We  recall  that  *  0.108  corresponds  to  the  height  of  the 

actual  test  specimen  in  rig.  2.  This  brings  us  back  to  the  initial  interpretation 
oi  Che  post -c r  i t  leal  response  daca  in  terns  of  the  nominal  sottcr.lr.g  model  and  the 
alternative  fracture  energy  formulation  for  the  localized  cracx  nand.  While  the 
consr  mt  strain-softening  approach  is  strictly  valid  only  for  -.he  particular  test 
configuration  under  cor.slderat  ion ,  the  constant  fracture  energy  re;ease  rat* 
concept  provijes  i  raclonal  me:. nodology  for  translating  this  result  to  other  test 
configurations  with  a  -ilr.gie  crack  band. 


2 . 2  Post -Cr  1 1 1  c 1 1  Re s ponse  Coap' it.irions  o i  Direct  Tension  Test 

During  the  discussion  of  the  direct  tension  test  uata  it  was  observed  that 
both  the  continuous  strain-softening  approach  as  well  as  the  discrete  fracture 
approach  raised  serious  suestiona  with  regard  to  their  validity.  In  ct:e  first 
case  it  was  doubtful  if  the  scraln-so: tcning  of  a  continuum  would  be  ihie  to 
capture  »  single  macrodeiect.  In  the  second  case  it  was  questionably  if  a  nr,- 
slant  fracture  energy  release  rate  would  be  applicable  to  specimen  si. ms  ot  the 
order  it  the  char  n: t  c  i  st  ic  fractur?  dixension. 

For  this  reason  let  us  investigate  the  post-critical  response  predictions  of 
ti..-  !  i  r  c  c  t  tension  t  >'  c  t  using  the  compos- tc  d-anag*  mode)  presence  aoovv.  x.,  this 
•  •nd  thr  upper  halt  ot  the  cylindrical  NX  specimen  in  rig.  2  Is  l  l  rallied  with 
three  different  meshes  -.*•  biquadratic  elements.  Thr  Initial  axi.  :  load- 
displacement  relation  corresponds  to  the  bilinear  Idealization  of  t ne  uniaxial 
tension  test,  whereov  the  specimen  Is  loaded  under  dlsp’ actmenc  cmtrol  with 
radial  restraints  at  the  loading  platens. 

Figure  *  depicts  the  computational  response  prediction  fo  ■  the  case  of 
constant  fracture  energy  release  rate  In  tension  i.e.  Cc  *  ■  -■  N/ra.  The 
relerence  'solution  agrees  fairly  close*./  with  the  results  o‘  *  ue  squar.  -  x  8  mesh 
as  well  as  the  rectangular  2x4  and  '■*  %  "2  idealizations  in  -.ilch  the  softening 
modulus  is  adjusted  according  to  cq .  I?  in  order  to  enforce  i  constant  fracture 
energy  release  rate.  The  constant  fracture  energy  release  rate  formulation 
captures  the  strain  localization  in  ,.he  single  crack  bend  tairl  well*  Figures  b 
and  o  cle»rly  illustrate  the  overall  motion  and  the  redlstr that . on  of  strains  and 
stresses  within  the  tensile  test  specimen  during  the  formation  ot  the  single 
r.ac  roue  I  v  c  t  .  In  fact,  the  comparison  b<  tween  the  resuirs  it  -smu  strength  ,md  the 
solution  it  low  residual  strengci*.  demonstrates  c  he  r  onsecucnr  e-j  of  strain 
i-)ul  izUl  ;n  whim  ini'- lutes  at  the  right  top  corner  of  the  tension  specimen 
N'c-uis**  of  the  stress  concent  i  at  ion  due  to  radial  restraint  at  the  loading  platen. 

Iu  contrast  the  c  onput  a  t  lonal  po*  t-cr  1 1  leal  predictions  tor  'he  •:  ise  or 
constant  so:t'*n|eg  modulus  %  -  0.82  exhibit  the  almost  perfectly  brittle 

r ••spouse  -h*. •  wn  ;n  rig.  7.  Hus  excessive  decrease  of  doctlllfv  Is  a  «il.vct 
cunse-punre  of  :  wn  -.esn  re  f  i  whlcn  Introduces  multi  pic  “racking  into  the 

l  Je.i  1  1  /a  t  i«*n  i:»st"iJ  -t  i  single  m.icroc  rar*  ,  The  si.-, pie  expression  In  h*j  .  I) 

«!•■-)  jr.‘,  t  r  aces  ...arlv  tl-.it  fans  1  r.g  ;  -p  mesh  height,  i.e..  o  ‘  i  /d  n, 

iTpilos  ua  e'Vi*/i:*nt  Increase  of  i*  jnd  the  conj  igif.*  frictufe  «»ntr.:y  release 
lit-*  •»..  -  f  lq  Is  to  remain  constant.  A_s  a  result  the  ••imr.gy  Jlssin.itlon 
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Flu.  5  ?lrecr.  tension  test  -  peak  response  it  load  -top  3 
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capacity  of  t  he  speol-ien  diminishes  to  the  limiting  condition  of  eLastlc  perfectly 
brittle  response  behavior  If  Che  height  ot  the  specimen  Is  subdivided  Into  3 
retried  lesh.  In  fact.  this  conclusion  Is  fairly  obvious  1 1  ve  recall  that  the 
Initial  softening  viope  Ks  ■  0.82  £  was  calibrated  hv  the  clement.irv 
lde.1L1s.ulon  of  the  iircct  tension  test  In  terms  o;  nominal  stress  ind  strain, 
wherebv  the  axial  strain  was  based  on  the  gauge  length  of  the  entire  specimen 
height  containing  a  single  tensile  crack.  Consequently,  it  is  Illusory  to  capture 
the  actual  local l/At Ion  of  the  tensile  crack  band  with  the  constant  strain 
sol  ten  mg  model. 
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'.'icr.cr leal  predictions  'or  c  constant 


1.  TRUXIAL  CCMPKESS ION  TEST 

While  the  fracture  mechanics  concepts  have  been  widely  accepted  for  tensile 
macro-erne <1 ng ,  there  Is  considerable  debate  with  regard  to  shear  fracture  of 
concrete  and  g«»o«t  j  te  r Ills.  In  fact,  most  geotechnical  failure  computations  are 
bavd  on  t  he  traditional  notion  of  ha rden 1 ng/so f ten lng  plasticity  theories 
ll-.lbj,  while  prominent  applied  mechanicians  advocate  fracture  mechanics  models 
for  : he  development  ot  shear  bands  in  geophysical  applications  (16tl7,18|. 

In  cemen t l t Iojs  materials  such  as  concrete  the  situation  Is  aggravated  by 
the  heterogeneous  composition  or  the  mesostructute  which  leads,  for  example.  to 
plittlnr.  tensile  failure  under  uniaxial  compression.  Moreover,  the  traditional 
shear  test”,  tonally  result  In  a  tensile  cracking  mode  ot  fi.actur*  ll  brittle 
materials  -ire  const Jer^d.  In  fact  a  frictional  slip  mode  of  failure  can  only  be 
observed  under  sufficiently  high  confining  stresses.  for  this  ;**..son,  the 
coiivniu-jn.il  tn.ixial  compression  cest  Is  examined  below,  in  wni-h  »  constant 
contlnlng  stress  Is  applied  first  before  the  specimen  I9  failed  under  axial 
1  ■:  1  1  i:ig. 

Klgurr  4  Illustrates  the  NX -size  specimen  and  .1  typical  response  curve  wnlch 
was  iM-iir-  d  l;t  the  -.^rvu-cont roi  ie-J  frlaxiii  doe*,  cell  devel  >pei  at  the  aivcril* 
tv  )t  Coiiralo,  boulder  for  pos t -c r 1 t leal  response  studies  [9].  The  following 
material  ;•  if  ’.meters  were  extracted  from  a  series  of  trlaxlal  compression  tests  it 
•I  1  *.  i<  rent  con: In  lag  pleasures 


—  -Q  — 


Stiffness:  C 


2,873  MPa 

0.20 


Strength:  C0  «  2.868  MPa 

♦o  -  4l-3* 


(IS) 


For  simplicity  let  os  assume  that  the  degradation  of  strength  is  primarily  a 
result  of  decoheston  while  the  internal  friction  remains  constant.  Therefore,  at 
the  residual  strength  level  the  cohesion  reduces  to  zero,  Cj.  ■  0,  and  the  angle 
of  Internal  friction  retains  Its  Initial  value,  *  41.3*. 


7lg.  9  Triaxtal  compression  test  of  hX-slze  specimen 

If  we  replace  the  localization  of  the  shear  band  by  an  equivalent  uniform 
distribution  of  microdetects  over  the  gauge  length,  then  we  can  Interpret  the 
softening  branch  In  terms  of  the  shear  softening  modulus 

C,  -  0.09  C  -  143.6  MPa  (16) 

As  a  result  the  pre-  and  post-critical  response  behavlur  Is  Idealized  by  a 
bilinear  stress-strain  relation  followed  by  perfectly  plastic  behavior  it  the 
residual  strength  levci.  Again  the  definition  of  the  strain-softening  modulus 
C,  Is  a  structural  property  which  Is  Intimately  tied  to  the  particular  test  and 
fracture  configuration.  Therefore,  It  Is  questionable  If  and  how  this  parameter 
tin  be  applied  to  dllferenc  geometries  and  boundary  conditions. 

On  the  other  hend,  the  displacement  controlled  development  of  a  single  shear 
hand  In  che  triaxiai  compression  specimen  can  be  utilized  to  calibrate  the 
critical  fra  fare  energy  release  rate  in  shear.  According  to  Palmer  and  Rice  |ld) 
the  Integral  under  the  nominal  inear  s t rcss-t angen 1 1  a l  displacement  curve  yields  a 
first  order  approximation  of  the  critical  shear  fracture  energy  release  rate  for 
the  particular  mortir  under  consideration.  Let  us  assume  that  a  homogeneous 
stress  slice  prevails  within  the  specimen  up  to  the  formation  of  a  single  shear 
hind  which  is  Inclin'd  lccordlng  to  the  ion r-Cou  1  orab  slip  condition  at.  the  angle  ■) 
*  » / 2  ♦  a/ with  toe  -,j jm  axis  of  principal  stress.  Furthermore,  If  we  neglect 
the  effect  if  elastic  shear  deformations  within  the  two  portions  of  Che  specimen 


separated  by  the  shear  band  shown  In  Fig.  8,  then  the  critical  shear  fracture 
energy  release  rate  can  be  evaluated  by  the  elementary  expression 
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TMs  extremely  large  value  Is  a  direct  result  of  the  tenfold  Increase  of  strength 
and  ductility  lr.  compression  versus  tension.  The  magnitude  agrees  with  first 
order  approximation  of  Gc  s  when  calibrated  froa  alternative  test  daca  of  the 
direct  shear  test.  The  conjugate  value  for  the  characteristic  fracture  dimension 
It  defined  by  the  analogout  expression  for  tensile  cracking  In  Eq.  5 
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If  we  adopt  again  the  notion  of  a  fracture  process  rone  with  degraded  material 
properties  Gj  then  the  width  of  the  shear  band  Is  actually 
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(19) 


However,  we  still  face  a  fracture  process  rone  which  is  extremely  large  and  of 
the  saec  sire  as  che  specimen  height.  Therefore,  In  the  case  or  frictional  slip 
It  Is  highly  questionable  If  the  fracture  mechanics  of  a  single  shear  band 
governs  the  post-crl t leal  response  behavior. 


). 2  Composite  Damage  Model  for  Frictional  Slip 

For  the  sake  of  a  unified  approach  towards  shear  softening,  let  us  resort  to 
che  composite  damage  mcdel  Illustrated  In  Fig.  9,  see  alio  refs.  (9-!l|. 
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Th«  elementary  volume  V  Is  again  composed  of  Che  undamaged  elastic  cone  Ve  and 
Che  localized  shear  band  V,  of  width  d ,  and  length  hs  In  which  ortenced  slip 
It  confined.  Resorting  to  the  Mohr-Coulomb  condition  of  frictional  slip,  the 
orientation  of  the  frictional  shear  band  is  defined  by  the  angle  9  *  ♦  */4 

with  the  major  axis  of  principal  stress.  Elementary  geometry  yields  for  the 
rectangular  configuration  of  the  elementary  volume 

V  •  V  ♦  V  where  V  »  a  b  c  and  V  »  2d  h  c  (20) 

e  S  5  9  3 


In  this  case  2da  denotes  the  width  and  h&  tht  length  of  the  shear  band  within 
the  elcoencary  veluae,  where 
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whichever  Is  smaller. 

for  straining  up  to  rupture  the  energy  release  of  the  composite  element  Is 
confined  to  che  localized  shear  band.  In  the  simple  case  of  linear  degradation 
of  the  cohesion  the  shear  energy  release  of  the  composite  element  follows  from 
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The  critical  fracture  energy  release  rate  Cor  advancing  the  shear  band  by  the 
surface  AA,  '  hs  c  Is  than 
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The  energy  dissipation  characteristics  In  the 
are  defined  by 
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equivalent  continuum  element 


(24) 


From  the  equivalence  AU^  •  AU^  the  softening  modulus  os  in  shear  depends  on 
the  fracture  energy  release  rate  Ccs  or  the  characteristic  fracture  dimension 
1,  and  the  cracK  geometry  according  to 
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This  expression  clearly  demonstrates  that  the  shear  softening  modulus  C,  of  the 
equivalent  continuum  element  depends  on  the  geometry  of  the  elementary  volume  and 
thus  exhibits  a  pronounced  size  effect.  In  other  words,  C5  depends  on  the 
particular  ilntt"  element  mesh  layout  If  the  fracture  energy  release  rate  <;c  , 
and  the  respective  characteristic  fracture  dimension  \s  constitute  invariant 
fracture  properties.  A3  a  result,  the  .under  lying  fracture  mechanics  concept 
loads  to  a  non-local  format  of  the  equivalent  softening  relation  which  la 
fundamentally  different  from  ths  local  constitutive  theories  of  straln-sottenlng 
pfasclcltv  and  continuous  damage  mechanics. 

Wj  conclude  that  the  composite  damage  models  for  tension  and  shear  have 
essential!/  the  same  structure.  In  both  cases  the  fracture  energy  release  concept 
results  In  a  non-local  format  of  the  equivalent  softening  continuum. 

Before  turning  to  the  computational  study  let  us  consider  again  the 
particular  lspects  it  the  crlaxlal  compression  lest  examined  previously  In  Section 
.  1  .  In  tin.-,  case  the  shear  softening  modulus  expression  of  Eq.  25  reduces  to 

C  ^  •  r-  — — -  where  1  •  1,465  m  and  9  *  65.05^ 

s  2l  a 
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Recall  that  a  sin  9  denotes  the  size  of  the  elementary  volume  transverse  to  the 
shear  band.  Clearly  for  reasonable  values  of  with  regard  to  the  NX-slze 
specimen  we  obtain  increasingly  ductile  behavior,  1 .e .  Gs  »  0  as  bl  *  0.  In 
fact,  for  a  single  element  Idealization  of  the  entire  specimen  with  b  •  0.108  m  we 
find  chat  G,  •  0.017  C  which  Is  somewhat  smaller  than  the  Initial  Interpretation 
of  the  post-crltlal  response  data  In  terms  of  the  nominal  shear  sottenlng  modulus, 
G,  •  0.05  C.  However,  from  this  staple  result  it  should  be  apparent  that  a 
convergence  study  with  successive  oesh  refinement  will  result  In  Increasingly 
ductile  response  predictions  if  lf  remains  Invariant.  M ternat  Ively ,  we  should 
ask  ourselves  what  would  happen  if  the  constant  fracture  energy  release  concept 
were  to  be  applied  to  the  solution  of  different  boundary  value  problems  In  view  of 
the  large  value  of  x,  •  l.*b5  a.  Clearly  If  sufficiently  large  geometries  would 
be  considered,  the  opposite  effect  would  take  place.  In  fact  for  Increasing  bl 
values  Increasingly  brittle  post-peak  predictions  would  result,  whereby  the 
maximum  mesh  size  tor  b1  Is  restricted  by  the  condition 

21 

— - -  >  l  (27) 

bl  cos6 

For  ■  t.w65  m  this  condition  results  In  the  Halting  value  tor  h*  <  7.106  a 
and  cue  respective  restriction  for  a*  <  J.JIb  m  for  perfectly  brittle 
conditions,  t.e.  C,.  •  •.  In  other  words  a  constant  fracture  energy  release 
analysis  of  large-scale  geotechnical  problems  would  result  In  a  mesh  size 
restriction  In  the  order  of  aeters. 

J . 2  Pos t-Cr 1 1 1 ca  1  Response  Computations  of  Tr laxlal  Coapres sion  Test 

In  view  of  Che  different  Interpretation  of  the  experimental  data  In  Che 
post-critical  regime,  let  us  examine  the  computational  response  predictions  for 
the  compression  test  using  the  composite  damage  model  above.  The  upper  half  of 
the  cylindrical  NX  specimen  Is  Idealized  with  three  different  mesh  layouts  of 
biquadratic  elements.  The  axial  load-displacement  relation  corresponds  to  the 
bilinear  idealization  of  the  crl.ixul  compression  test,  whereby  the  specimen  Is 
loaded  under  axial  displacement  control  with  radial  restraints  at  the  loading 
platens  and  radial  confinement  o  *  0.69  MPa . 

Figure  10  Illustrate:.  the  computational  response  predictions  for  the  case 
where  che  softening  modulus  reaalns  Invariant,  t.c.  C.3  •  0.05  C.  The  numerical 
result;  Indicate  that  the  post-critical  response  Is  virtually  Independent  of  the 
particular  mesh-slze  and  underestimates  the  ductility  of  the  Idealized  test  data 
In  Fig.  8.  This  conclusion  coincides  with  previous  observations  110,11)  for 
distributed  shear  fracture  which  can  be  auaoarlr  »d  by  stadng  that  the  constant 
softening  formulation  yields  pos t -cr 1 1 lea l  response  predictions  which  are  quite 
Consistent  and  Independent  of  mesh  size.  The  constanc  softening  modulus 
formulation  captures  fracture  In  the  fora  of  a  shear  cone  and  radial  splitting 
within  the  compression  speclaen.  Figures  11  and  12  clearly  illustrate  the  overall 
motion  and  '.he  strain  loc a 1 1  z.a t Ion  along  the  cone,  whereby  the  axlsymmetrlc 
’•deal  l  zat  Ion  imposes  some  restriction  on  the  failure  mode.  However,  the 
comparison  between  the  results  at  peak  strength  and  the  solution  at  the  residual 
strength  level  etoarlv  demonstrates  the  formation  of  the  compressive  shear  cone  at 
the  loading  platen  and  the  subsequent  radial  splitting  at  mid-height. 

In  :on:r.i<J  l  st  :  nc  c ;  in  Che  basic  fracture  siechanlcs  postulate  of  an  Invariant 
••nertfy  release  rate  lead*  to  the  results  depleted  In  rig.  13.  The  excessive  duct¬ 
ility  during  the  forn.it  Ion  of  the  shear  band  leads  to  a  post-critical  response 
predict!  n  which  exhibits  very  little  softening.  Thp  reason  for  the  increased 
ductility  has  be^n  !i  icussed  in  the  context  of  F.q .  2b.  In  fact,  the  large  value 
of  the  characteristic  fracture  dimension  «  l.whi  a  sloply  indicates  chat 
there  Is  no  1  or  ai  \  r.\  c  l  '*n  of  a  single  shear  band  possible  In  the  NX-slre  test  spec¬ 
imen  with  »  height  of  'i.hw'8  *n .  In  other  words  the  entire  specimen  constitutes  the 
fracture  i-ncss  .tone  which  Is  mobilised  In  the  energy  exchange  during  shear  fail¬ 
ure.  \s  i  result  the  o^stant  sotconmg  modulus  for-iulatlon  ror  uniformly  llstri- 
bu t ••  J  ‘ale rode t ecc s  i.  in  this  case  core  appropriate. 

We  on  c  I  u.»»*  t:i-»c  frictionil  slip  rv;>#  fracture  takes  place  In  cementitious 
materials  ;n  i  :ar  -.ore  distributed  tashion  than  tensile  granting.  in  fact,  a 
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single  shear  band  develops  only  under  sufficiently  confined  conditions  when 
prenature  cracking  Is  suppressed  by  compressive  stress  states.  Therefore  fhe 
fracture  oechanlcs  hypothesis  of  a  constant  fracture  energy  release  rate  In  shear 
la  not  of  Immediate  concern  cor  tradlclonal  structures  but  La  of  paramount 
Interest  in  geocechnlcal  and  geophysical  applications. 

4.  DIRECT  SHEAR  TEST 

In  the  final  example  let  us  examine  the  post-critical  behavior  of  the  direct 
shear  test  specimen  In  which  both  tensile  cracxlng  and  frictional  slip  take  place 
s tnul taneously .  In  this  case  distributed  fracture  Is  contained  by  the  two  steel 
box  compartments,  whose  motions  are  monitored  by  two  servo-controlled  actuators 
In  the  stiff  shear  test  frame  built  at  the  L'olverslty  of  Colorado,  Boulder  [20]. 

In  this  sense  the  direct  shear  test  specimen  conatltutea  a  miniature  structure  In 
which  contained  fracture  propagates  In  a  manner  that  la  similar  to  the  failure 
behavior  In  confined  structures  of  high  Indeterminacy.  The  normal  pressure  Is 
applied  to  a  prescribed  load  level  and  held  conscant  during  me  subsequent 
tangential  loading  monitored  under  displacement  control  In  order  to  capture  the 
full  post-critical  response  branch  down  to  Che  residual  load  level.  Depending  on 
che  rado  of  normal  to  tangential  load  level  the  full  range  of  predominantly 
tensile  cracking  to  highly  confined  frictional  slip  Is  mobilized  within  the  'shear 
test"  specimen. 

From  the  different  test  results  reported  in  Ref.  [10]  we  consider  the  prls- 
matlc  0.204  x  0 . 102  x  0.117  o  mottar  specimen  tested  at  the  normal  pressure  an 
■  1.379  MPa.  The  test  specimen  Is  bonded  Co  the  steel  box  compartments  with  the 
aid  of  an  epoxy  compound  fn  order  to  ensure  tull  contact  on  the  tensile  sidewalls 
of  the  shear  box.  Because  of  the  planar  mode  of  fracture  a  two-dimensional  plane 
strain  Idealization  Is  used  to  model  the  specimen  and  contiguous  shear  box.  The 
finite  elemenc  aesn  layout  Includes  the  actual  test  specimen  as  well  as  che  upper 
and  lower  shear  box  compartments  In  order  to  ensure  that  the  tangential  displace¬ 
ment  control  lies  In  tha  prevalent  shear  plane  st  mid-height  of  the  specimen. 

The  Idealization  of  the  direct  shear  test  Is  shown  In  Pig.  14  for  the  3x6  mesh 
of  the  specimen.  Altogether  683  degrees  of  freedom  were  used  to  Interconnect  73 
biquadratic  Isoparametric  elements  which  Include  tha  Idealization  of  the  steel  box 
ana  the  loading  frame  In  order  to  avoid  overturning  effects  due  to  eccentric 
tangential  loading,  see  also  ref.  (21). 

The  composite  damage  model  for  tensile  cracking  and  frictional  slip  are 
combined  for  describing  constant  softening  or  alternatively  constant  fracture 
energy  release  In  tension  and  shear.  The  actual  properties  of  ihc  mortar 
specimen  were  calibrated  before  In  the  context  of  the  direct  tension  and”  trlaxial 
compression  ccscs.  Therefore  che  direct  shear  test  constitutes  an  entirely  Inde¬ 
pendent  example  of  different  size  which  will  asserc  the  previous  findings  on  dis¬ 
tributed  and  localized  fracture  In  tension  and  shear.  The  elastic  response  Is 
defined  as  before  by  E  *  6,895  MPa  and  v  ■  0.20.  The  Initial  fracture  condition 
Is  governed  hy  che  modified  Kohr-Coulomb  condition  Cn  •  2.87  MPa  and  *  41.3* 
augmented  by  the  tension  cut  o'f  o{  *  1.59  MPa.  The  post-critical  behavior 
Is  in  one  case  described  by  the  constant  softening  model  where  £,  •  0.82  E  and 
G,  -  0.05  C  and  In  the  other  case  by  the  constant  fracture  energy  release  rates 
Cr.  ,  t  *  N/m  and  Cc  _  3  *  2,452  N/m.  In  addition,  a  third  post-critical 
strategy  Is  pursued  in  which  the  fraccure  energy  release  rate  In  tension  is  kept 
r-oiscnnt,  t  *  G  u/a. while  the  shear  softening  modulus  remains  Invariant, 

G,  •  o.O 5  G.  rot  frictional  slip  a  non-assoc la  ted  flow  rule  Is  used  throughout 
In  order  to  capture  che  splitting  tension  mode  of  fracture  during  shearing  with  a 
flow  potential  In  the  form  of  the  maximum  stress  criterion.  In  the  case  of 
tensile  cracxing  this  flow  rule  reduces  to  the  traditional  associated  flow  rule 
used  to  monitor  the  growth  of  tensile  cracking  strains. 

Figure  15  depicts  the  overall  load-deformation  behavior  In  terms  of  the  noml - 
nal  shear  stress  In  the  prevalent  shear  plane  versus  the  relative  tangential  dis¬ 
placement  hetween  the  upper  and  lower  shear  box  compartments.  The  three  compul  i- 
clon.il  predictions  exhibit  considerable  variation  as  coopered  to  the  experlfli  ntal 
results.  The  ,-onsCanc  shear  soften!"*;  modulus  approach  shows  the  best  agreement 
with  regard  to  both  strength  and  overall  dual  1 1  ty,  whl  le  the  constant  fracture 
energy  formulation  Indicates  far  too  ductile  behavior.  Clearly,  In  view  of  the 
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Fig.  15  Direct  shear  test  -  numerical  predictions  for  various 

softening  strategies 


confined  fracture  propagation  within  Che  shear  specimen  one  could  argue  that  a 
single  slip  band  only  develops  In  the  final  stage  at  the  residual  strength  level. 

In  addition,  this  locallied  shear  band  Involves  extensive  tensile  cracking  In  a 
direction  which  differs  substantially  from  the  prevalent  shear  plane.  Theiefore, 

one  could  conclude  that  contained  fracture  Is  beyond  the  realm  of  fracture 

•echanlcs  of  single  aacrodefects.  However,  In  engineering  applications  this  (Inal 
node  of  fracture  In  cementitious  materials  Is  always  accompanied  by  extensive 
tones  of  distributed  damage  In  which  the  constant  aoftenlng  model  prevails. 

Figure  16  Illustrates  the  shearing  motion  In  terms  of  the  total  displacement 
vectors  for  load  steps  corresponding  to  the  peak  and  post-critical  regimes. 

Observe  the  localization  Into  a  wide  shear  band  which  Is  accompanied  by 
considerable  vertical  dilatation  due  to  tensile  splitting.  The  total  strain 
distribution  Is  depicted  In  Fig.  17  at  the  same  loading  steps.  The  plots 
emphasize  the  localization  of  shear  strains  In  the  prevalent  shear  plane  while 
the  upper  and  lower  portions  of  the  specimen  unload  progressively.  The 
associated  stress  redistribution  during  progressive  fracturing  is  shown  in  Fig. 

18  with  plots  of  principal  stresses  at  the  two  load  steps.  Recall  that  the  unor 

principal  stress  Is  Initially  aligned  with  the  vertical  pressure  on  *  1.379 

MPa.  It  rotates  upon  continuous  shearing  by  almost  80*  In  order  to  accommodate 
the  cccpresslve  strut  action  between  the  compressive  slde-valls  of  the  shear  bov. 
The  Initial  cracking  Is  confined  to  the  tensile  corners  of  the  specimen  and 
exnlblts  little  Influence  on  the  overall  response  behavior.  The  splitting  renslon 
that  Is  parallel  to  the  compressive  st'uts,  leads  ultimately  to  a  single  shear 
band  at  mid-height  In  the  form  of  an  S-shaped  fracture  surface. 

The  results  clearly  confirm  r '  he  direct  shear  test  is  definitely  not  a 
slcple  material  test  when  cement  1 1 io-s  materials  are  considered.  In  fact,  t  -e 
contained  mode  of  fracture  results  In  a  complex  combination  of  tensile  cracking 
and  frictional  slip  which  la  highly  distributed  before  final  localization  of  a 
macrodt'fect  taxes  place. 


Pit*  16  Motion  of  the  direct  sheer  test  specimen 


PRINCIPAL  STRESSES 


FIr.  17  Stress  redistribution  within  direct  shear  test  specimen 


PRINCIPAL  STRAINS 


Fis.  13  Strain  redistribution  within  direct  shear  test  soeclnen 
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5.  CONCLUDING  REMARKS 


The  principal  result  of  the  proposed  composite  damage  model  Is  a  unified 
approach  cowards  localized  and  distributed  fracture  In  Che  fora  of  tensile 
cracxltag  and  frictional  slip.  In  the  case  of  localized  fracture  the  composite 
damage  model  provides  a  non-local  softening  law  of  Che  equivalent  continue  which 
depends  on  the  characteristic  fracture  dimension  as  well  as  the  geometric 
configuration  of  the  fracture  band  within  the  elementary  volume  of  finite 
dimension.  In  fact,  the  composite  damage  model  recovers  In  this  case  the  constant 
fracture  energy  release  concept  without  resorting  to  the  fracture  mechanics  of 
discrete  oacrodefects.  On  the  other  end  of  the  spectrtm  the  composite  damage 
model  degenerates  to  che  continuous  damage  theory  of  distributed  fracture  when  the 
volume  fracdon  of  localized  damage  approaches  unity.  In  this  case  the  softening 
modulus  Es  of  the  equivalent  continuum  constitutes  a  local  material  property 
which  holds  polntvtse  and  Is  Independent  of  the  particular  fracture  configuration. 

The  cocposlte  damage  aodel  was  applied  to  the  post-critical  response  behavior 
of  the  direct  tension  and  crlaxlal  compression  tests  on  cylindrical  mortar 
specimens  for  which  test  data  were  available  for  comparison.  While  the  constant 
fracture  energy  release  concept  captured  the  localization  of  a  single  crack  band 
In  tension,  che  constant  softening  Qodel  reproduced  distributed  fracture  In  shear. 

In  the  final  example  of  the  direct  shear  test  both  tensile  cracking  as  well 
as  frictional  slip  modes  of  contained  fracture  were  mobilized,  similar  to  the 
fjllure  behavior  In  highly  confined  structures.  In  this  case  the  combination  of 
the  constant  fracture  release  concept  for  tensile  cracking  and  the  constant 
softening  formulation  for  snear  slip  furnished  the  best  agreement  with  the 
experimental  post-critical  response  behavior. 

Clearly,  the  previous  exposition  was  based  on  gross  slmplleatlon  In  order  to 
arrive  at  elementary  algebraic  relations  wnlch  clearly  exhibit  the  fundamental 
differences  of  dlscrlbuced  and  highly  localized  fracture.  In  reality  fracture  Is 
an  evolutionary  process  In  which  the  localized  damage  zone  Increases  gradually  and 
where  the  critical  fracture  energy  release  rate  varies  In  accordance  with  the 
experimental  degradation  of  strength  in  the  post-critical  response  regime.  In 
addition  Inelastic  dissipation  mechanisms  take  place  which  are  mobilized 
particularly  In  shear,  As  such  the  composite  damage  model  above  Is  a  first  order 
approximation  which  captures  the  limiting  conditions  of  the  local  si ze- Independent 
softening  relation  as  well  as  the  non-local  fracture  energy  release  rate  concepc. 
For  the  continuous  transition  from  distributed  to  localized  fracture  an 
evolutionary  damage  process  should  provide  an  Immediate  solution  to  combine  these 
two  extreme  cases. 
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EXPERIMENTAL  AND  CONSTITUTIVE  ASPECTS  OF  CONCRETE  FAILURE 


By  Kaspar  Wiliam1,  M.ASCE,  Bryan  Hurlbut2,  and  Stein  Sture3,  M.ASCE 


ABSTRACT 

The  failure  analyses  of  concrete  structures  are  plagued  (1)  by  the  con¬ 
stitutive  difficulties  in  describing  the  transition  from  brittle  to  duc^ 
tile  fracture  of  concrete  under  increasing  confinement,  and  (ii)  by  the 
nunerlcal  intricacies  related  to  structural  response  predictions  of  con¬ 
tained  material  instabilities  due  to  strain-softening  and/or  due  to  lack 

#  of  normality  in  conjunction  with  non-associated  plastic  flow  concepts. 

•  This  paper  presents  novel  experimental  observations  of  post-critical 
response  phenomena  when  cylindrical  concrete  specimens  are  loaded  under 
displacement  control  in  the  form  of  direct  tension  and  triaxial 
compression  experiments.  The  investigation  focuses  on  the  problem 

9  whether  or  not  "proper"  material  properties  can  be  extracted  from  such 

experiments  which  characterize  not  only  stiffness  and  strength  behavior 
j  but  also  the  ductility  in  the  post-peak  response  regime.  In  this 

i  context  a  fracture  energy  based  strain-softeni  ig  plasticity  formulation 

j  will  be  developed  in  terms  of  the  fundamental  Mode  I  type  cracking  which 

•  accounts  for  the  transition  from  brittle  to  ductile  fracture  under 

%  increasing  triaxial  confinement. 

INTRODUCTION 

The  failure  behavior  of  cementitious  materials  such  as  mortar  and 
concrete  can  be  broadly  classified  into  (i)  tensile  cleavage  in  the  form 

#  of  Mode  I  type  cracking,  and  (li)  decohesive  faulting  in  the  form  of 
Mode  II  type  frictional  slip.  The  tensile  fracture  mode  typically  leads 
to  a  sudden  discontinuity  at  peak  strength  and  a  very  brittle  appearance 

j  of  the  ioad-deforraation  response  when  the  direct  tension  experiment  is 

considered.  Recent  experiments  on  notched  concrete  specimens  by 
Reinhardt  [15],  Gopalaratnan  and  Shah  [6]  as  well  as  the  direct  tension 

#  experiments  by  Hurlbut  [9]  on  unnotched  specimens  all  indicate  that 
little  hardening  precedes  the  onset  of  cleavage  failure  which  results  in 
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the  separation  of  the  test  specimen  into  two  competent  parts.  In 
contrast,  shear  fracture  exhibits  far  more  ductile  behavior,  whereby 
hardening  in  the  pre-peak  regime  is  followed  by  a  gradual  softening 
regime  during  the  process  of  progressive  decohesion  down  to  the  residual 
strength  level  which  is  characterized  by  granular  friction  only.  The 
triaxial  compression  experiments  on  concrete  by  van  Mier  (13)  and 
Hurlbut  (9)  illustrate  clearly  the  increase  of  ductility  with  increasing 
confining  pressure.  In  fact  the  splitting  tension  mode  of  failure  in 
uniaxial  compression  transforms  gradually  into  a  faulting  mode  of 
failure  with  increasing  confinement  and  ultimately  diffuses  into 
distributed  damage  throughout  the  entire  specimen.  As  a  result  under 
sufficiently  high  confinement  there  is  a  third  response  regime  which 
does  not  exhibit  any  softening  whatsoever  and  which  is  governed  by 
corapactive  hardening  due  to  pore  collapse  and  related  Inelastic 
processes  of  the  composite  microstructure  of  concrete. 

Since  the  failure  of  plain  and  reinforced  concrete  structures  depends  to 
a  large  extent  on  the  ductility  and  thus  the  governing  post-peak 
response  characteristics  of  the  concrete  the  main  thrust  of  this 
investigation  is  directed  towards  the  softening  behavior.  We  will  focus 
on  recent  experimental  evidence  obtained  at  the  University  of  Colorado 
within  a  series  of  displacement  controlled  post-critical  experiments  on 
low  strength  concrete  specimens  featuring  shear  dilatancy,  hardening  as 
well  as  softening  in  tension  and  shear.  There  have  been  various 
attempts  to  correlate  these  global  phenomena  to  different  processes  in 
the  raicrostructure  such  as  debonding  between  the  cement  matrix  and  the 
aggregate,  progressive  microcracking  within  the  cement  matrix  and 
ultimately  fracture  of  the  aggregate.  Quite  noteworthy  is  the  recent 
work  by  Ortiz  (14]  who  used  the  mixture  theory  of  interacting  continua 
to  model  these  evolutionary  processes  within  each  phase  of  the 
composite.  This  development  Is  closely  related  to  the  continuous  damage 
formulation  of  Mazars  (12]  who  expressed  Mode  1  type  cracking  and  Mode 
II  type  shear  fracture  In  terras  of  maximum  tensile  strain  conditions. 

In  the  following  we  will  concentrate  on  the  fundamental  issues  of 
strain-softening,  the  transition  from  brittle  to  ductile  fracture,  and 
the  decrease  of  shear  dilatancy  with  increasing  confinement.  Special 
attention  will  be  paid  to  the  extraction  of  fracture  properties  from  the 
post-critical  response  regime  which  are  Independent  of  a  particular  test 
geometry  and  the  given  set  of  boundary  conditions.  This  task  is  compli¬ 
cated  by  the  fact  that  the  post-critical  response  behavior  is  accompan¬ 
ied  by  severe  strain  localization  and  actual  discontinuities  within  the 
test  specimen.  For  this  purpose  we  will  expand  the  composite  fracture 
formulation  which  has  previously  been  developed  by  the  authors  (19]  in 
order  to  incorporate  the  fracture  energy  release  rate  concept  within  the 
strain-3of tening  formulation  of  an  equivalent  elastic  plastic  continuum. 

POST-PEAK  EXPERIMENTS 

In  conjunction  with  a  substantial  research  effort  on  finite  elements  and 
localized  failure  in  cementitious  materials  two  series  of  experiments 
were  recently  completed  which  Involved  stable  post-crit J.caL  loading  of 
cylindrical  concrete  specimens  in  direct  tension  and  triaxial  coropres- 


slon,  (NX  size  samples  with  D  “  2.125  in.  (54  mm)  diameter  and  H  -  4.25 
in.  (108  mm)  height).  The  concrete  mix  was  based  on  a  water-cement 
ratio  of  W/C  »  0.833,  a  mix  grading  of  C:S:G  -  1:3.16:3.19  and  a  maximum 
aggregate  size  of  3/8  in.  (9.5  mm)  which  in  turn  gave  a  uniaxial 
compressive  strength  of  f'c  •  3.2  ksi  (22.1  MPa)  at  28  days.  The 
concrete  was  cast  into  1x3x3  ft.  (0.3  x  1.0  x  1.0  ra)  slabs  which 
were  kept  in  a  fog  room  until  testing  except  for  the  period  for  prepar¬ 
ation  of  the  specimens.  The  cylindrical  concrete  samples  were  cored 
from  the  slab  and  milled  with  a  diamond  grinding  wheel  to  the  prescribed 
height  with  a  tolerance  of  AH  ■  ±0.02  in.  (±0.5  mm).  The  exposed  voids 
at  the  cylindrical  surface  of  each  specimen  were  filled  with  putty  in 
order  to  prevent  intrusion  of  the  elastomeric  membrane  into  the  concrete 
voids  in  the  case  of  high  confining  pressures. 

Two  series  of  post-peak  experiments  were  carried  out  with  the  aid  of  a 
MTS  loading  system  consisting  of  a  standard  110  kip  (489  kN)  loading 
frame  with  servo  control  and  function  generator  units.  These  were 
comprised  of  a  series  of  direct  tension  and  a  series  of  triaxlal 
compression  experiments  which  will  be  described  below. 

Direct  Tension  Experiment 

In  order  to  clarify  our  understanding  of  the  tensile  fracture  mechanisms 
the  following  objectives  were  pursued  with  the  aid  of  direct  tension 
experiments: 

(i)  Uniqueness  of  the  load-displacement  response  in  the 
post-peak  regime, 

(il)  Progressive  damage  and  stiffness  degradation  during  loading 
and  unloading  cycles,  and 

(ill)  Size  effects  related  to  different  height  specimens. 

A  special  technique  was  developed  by  Hurlbut  [9)  to  stabilize  the 
post-critical  response  behavior  of  the  cylindrical  tension  specimens. 

The  basic  difficulty  of  maintaining  stable  post-peak  control  was 
resolved  by  Incorporating  two  axial  Linear  Virtual  Displacement 
Transducers  (LVDT)  which  were  attached  to  the  loading  platen  as 
illustrated  in  Fig.  1.  In  contrast  to  the  recent  tension  tests  by 
Reinhardt  (15]  and  Shah  [6]  unnotched  specimens  were  utilized  in  order 
to  avoid  any  bias  towards  fracture  mechanics  by  introducing  a  well 
defined  macrofect  in  the  sample  to  start  with.  Therefore  the  actual 
location  of  tensile  cracking  was  undetermined  and  the  stroke  of  the 
loading  actuator  had  to  be  controlled  by  the  two  axial  LVDT's  in  order 
to  monitor  the  post-peak  response  behavior.  In  return  a  premature 
fracture  mechanics  interpretation  of  the  tensile  cracking  process  was 
avoided  and  the  question  of  a  volume  versus  surface  dominated  failure 
process  could  be  explored. 

It  1 3  important  to  keep  in  mind  that  the  tensile  test  specimen  was 
placed  Into  the  loading  frame  12  hrs  before  actual  testing.  The 
cylindrical  sample  was  attached  to  the  specially  prepared  solid  loading 
platens  by  means  of  a  thin  layer  of  a  high  strength  structural  epoxy 
cur'd  during  the  12-hr.  period  under  a  slight  compressive  load  of  Pz  - 
-101)  lbs  (-445N).  Clearly  the  epoxy  layer  provided  some  lateral 
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restraint  and  was  therefore  responsible  for  a  non-uniforra  stress  state 
near  the  end  section.  In  fact,  the  stress  concentration  at  the  corner 
of  the  specimen-platen  interface  led  In  some  cases  to  premature  bond 
failure  or  delamination. 

Fig.  2  depicts  the  axial  load-extension  response  of  a  typical  direct 
tension  experiment  conducted  at  a  displacement  rate  of  1.67  x  10-6 
in/sec.  (4.24  x  10~8  m/sec).  We  observe  a  nearly  linear  ascending 
branch  up  to  the  peak  load.  At  this  point  the  axial  tension  suddenly 
decreased  without  noticeable  formation  of  a  crack.  Upon  close 
Inspection  a  single  transverse  hairline  crack  appeared  around  the 
periphery  of  the  sample  at  an  overall  axial  displacement  of  1.5  x  10-3 
In.  (0.04  mm)  when  the  residual  tensile  strength  diminished  to  near  zero, 
load  level.  At  that  time  the  specimen  had  separated  into  two  competent 
parts  which  remained  entirely  intact.  Upon  continuous  loading  the 
increasing  displacements  reflected  simply  further  crack  opening 
accompanied  by  rigid  body  movement  of  the  two  separated  portions  of  the 
specimen. 

Considering  the  reproducibility  of  the  post-critical  response  branch 
under  raonotonic  loading  the  question  arises  If  the  residual  tensile 
strength  can  be  recovered  during  loading-unloading  cycles  and  if  the 
associated  stiffness  is  affected  by  the  degrading  strength. 
Loading-unloading  experiments  were  performed  in  order  to  explore  the 
residual  tensile  strength  and  the  associated  tangent  stiffness.  Fig.  3 
depicts  the  resulting  load-displacement  response  when  the  stroke  is 
reversed  to  zero  displacement  and  subsequently  Increased  again.  The 
plot  Illustrates  clearly  that  tensile  bonds  still  remain  intact  In  the 
softening  branch  and  that  the  residual  tensile  strength  is  fully 
recovered  during  reloading.  In  fact  the  post-peak  load-displacement 
response  curve  approaches  that  of  the  raonotonic  experiment  very  much 
along  the  line  of  strain-softening  elastic-plastic  behavior  with 
substantial  permanent  deformations.  We  observe  a  slight  decrease  of 
elastic  stiffness  In  the  first  cycle  and  a  substantial  reduction  when 
the  tensile  crack  has  fully  developed  very  much  along  the  results  on 
notched  tension  specimens  reported  recently  by  Reinhardt  [15].  There  Is 
some  elastic-inelastic  coupling,  however,  it  Is  not  nearly  so  pronounced 
as  to  fully  support  the  progressive  fracturing  theory  of  Dougill  (5]  or 
the  continuous  damage  theory  of  Mazars  (12),  Janson  and  Hult  (10], 
Krajlnovlc  and  Fonseka  (11]  which  resort  to  secant  formulations  and  thus 
fully  reversible  deformations  upon  unloading. 

In  order  to  probe  the  fracture  mechanics  hypothesis  of  a  single 
macrodefect  an  additional  raonotonic  tension  test  was  carried  out  with  a 
specimen  of  half  the  height  and  the  same  cross-section.  In  the  case  of 
a  constant  fracture  energy  release  rate  the  area  under  the  P-u  diagram 
(actually  the  nominal  axial  stress-extension  plot)  should  remain 
constant,  assuming  that  the  tensile  strength  is  not  affected  by  the  two 
different  volumes  and  their  different  probability  of  failure.  On  the 
other  hand  if  the  concept  of  a  continuous  distribution  of  microdefects 
were  to  be  adopted  along  the  line  of  the  smeared  cracking  approach,  then 
the  area  under  the  P-u  diagram,  i.e.,  the  strain  energy  dissipation 
during  the  transverse  crack  process  should  be  reduced  by  half.  Fig.  4 
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illustrates  the  results  of  the  half-size  specimen  test.  The  plot 
clearly  indicates  that  the  area  under  the  P-u  diagram  is  hardly  affected 
by  halving  the  height  of  the  tension  specimen.  Therefore  the  results 
support  the  postulate  of  a  constant  fracture  energy  which  ultimately 
controls  the  localization  of  distributed  microdefects  into  a  single 
macrodefect. 

Trlaxlal  Compression  Experiments 

In  order  to  clarify  our  understanding  of  shear  fracture  the  following 
objectives  were  pursued  with  the  aid  of  axisymmetric  trlaxlal 
compression  experiments: 

(i)  Transition  from  brittle  to  ductile  fracture  at  various 
levels  of  confinement, 

(ii)  Shear  dilatancy  during  progressive  shear  fracture,  and 
(ill)  Stiffness  degradation  during  loading  -  unloading  cycles. 

Conventional  trlaxlal  compression  experiments  were  conducted  in  a 
Hoek  cell  device  (8)  which  had  been  modified  at  the  University  of 
Colorado,  Boulder,  in  order  to  measure  axial  and  lateral  deformations  of 
cylindrical  NX-size  specimens.  The  Hoek  cell  which  is  illustrated  in 
Fig.  5  has  a  capacity  of  up  to  ar  *  10  ksi  (68.9  MPa)  confining 
pressure.  It  is  instrumented  for  measuring  lateral  deformations  at  the 
mid-height  with  the  aid  of  six  strain  gaged  cantilevers  spaced  at  60° 
Intervals  around  the  circumference.  The  cell  pressure  is  provided  by 
hydraulic  fluid  and  is  transmitted  to  the  specimen  by  a  flexible 
polyurethane  membrane  of  0.05  in.  (1.3  ram)  thickness.  The  confining 
pressure  has  no  component  in  the  axial  direction  of  the  cylindrical 
specimen.  The  axial  load  was  applied  with  the  same  110  kip  (489  kN)  MTS 
loading  system  with  servo  control  and  function  generator  units  capable 
of  stabilizing  the  post-critical  softening  branch. 


!•' guru  0  Tri.ixi.il  Compression  Appoi  nt  u.-; 


* 


8 


All  the  trlaxial  compression  tests  were  conducted  by  pressurizing  first 
the  cell  to  the  prescribed  confining  pressure  which  was  held  constant 
while  the  axial  Load  was  applied  under  stroke  control  at  a  displacement 
rate  of  1.67  x  10-4  ln/sec  (4.24  x  10~6  m/sec). 

Fig.  6  depicts  the  axial  load-displacement  behavior  at  four  different 
confining  levels  o  -  100,  500,  1000  and  2000  psi  (0.69,  3.45,  6.89 
and  13.79  MPa).  Again  experimental  results  are  plotted  in  the  form  of 
global  load-displacement  results  rather  than  In  the  form  of  stress  and 
strain  measures  since  localization  of  fracture  takes  place  in  the 
softening  regime  which  requires  a  proper  homogenization  concept.  The 
four  plots  clearly  demonstrate  the  increase  of  ductility  with  increasing 
confinement  which  suppresses  the  softening  tendencies.  The  partial 
unloading  cycles  indicate  very  little  hysteretic  energy  dissipation. 

The  linear  elastic  unloading-reloading  behavior  hardly  exhibits  any 
degradation  cf  stiffness.  This  suggests  that  the  elastic-plastic 
coupling  remains  negligible  and  conjectures  little  damage  accumulation 
of  the  elastic  material  texture  not  only  in  the  pre-  but  also  in  the 
post-peak  regime.  The  resulting  failure  modes  showed  a  distinct 
transition  from  splitting  tension  under  low  confinement,  over  shear 
faulting  under  intermediate  confinement,  to  diffuse  damage  at  high 
confinement  when  no  softening  was  observed  altogether. 

The  dilatation  characteristics  are  shown  in  Fig.  7  in  the  form  of 
lateral  versus  axial  displacement  response  curves.  Again  the  lateral 
expansion  diminishes  rapidly  with  increasing  confinement  pressure.  This 
indicates  that  volumetric  dilatation  accompanies  strain-softening  and 
that  the  transition  from  brittle  to  ductile  fracture  is  related  to  the 
transition  from  volume  dilation  to  volume  compaction. 

If  the  axial  and  radial  displacements  in  Fig.  7  were  converted  into  nom¬ 
inal  strains  their  ratio  would  remain  fairly  constant  in  the  softening 
regime  and  would  decrease  rapidly  with  increasing  confinement.  It  was 
therefore  felt  that  the  inelastic  dilatancy  ratio  a  ■  eV  could 

provide  a  proper  fracture  characteristics  which  would,  however,  depend 
on  the  confining  stress.  The  post-peak  response  at  low  confinement  is 
accompanied  by  large  inelastic  volume  dilation  which  decreases  rapidly 
under  increasing  confining  stress  and  results  in  inelastic  compaction 
rather  than  dilation. 

The  trlaxial  strength  results  are  summarized  in  Fig.  8  with  a  plot  of 
the  test  data  in  the  Rendulic  plane  where  the  data  are  shown  on  the 
compressive  meridian.  The  diagram  depicts  the  maximum  strength  envelope 
in  conjunction  with  the  residual  strength  data.  The  residual  strength 
diminishes  to  zero  at  zero  confinement  and  approaches  rapidly  the 
maximum  strength  envelope  with  increasing  confinement  at  the  transition 
point  from  brittle  to  ductile  fracture  when  softening  i3  suppressed 
altogether.  Since  the  entire  test  program  was  focused  on  the  softening 
behavior  under  low  confinement  no  attempt  was  made  to  extract  the 
hardening  behavior  which  had  been  explored  in  a  separate  research 
effort  by  Stankowski  and  Gersile  (17). 
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COMPOSITE  FRACTURE  MODEL 


Up  to  now  the  test  results  were  Intentionally  described  In  terras  of 
overall  load-displacement  diagrams  because  of  the  strain  localization 
during  progressive  fracturing  in  the  softening  regime.  Therefore,  we 
have  to  ask  ourselves  whether  or  not  useful  information  can  be  extracted 
from  the  softening  data  which  is  net  only  generic  to  the  particular  test 
specimen  but  which  can  also  be  used  in  the  analysis  of  actual 
structures.  In  order  to  clarify  the  underlying  issues  let  us  consider 
the  "simple"  configuration  of  the  direct  tension  experiment.  In  this 
case  the  elementary  study  of  an  axial  tension  specimen  yields  some 
surprising  results  with  far-reaching  conclusions. 


Cons  ider 
with  the 
relation 


a  uniaxial  tension  specimen  of  height  H  and  cross-section  A 
uniaxial  tensile  strength  ot ,  The  incremental  stiffness 
follows  the  elementary  axial  force-deformation  relation 
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where  Ep  denotes  the  tangent  modulus,  az  the  nominal  stress  rate  and 
u z  the  axial  rate  of  extension.  If  we  assume  a  simple  bilinear 
material  law  with  Ep  ■  E  for  hardening  and  E-p  ■  -E  for  softening, 
uniform  degradation  of  strength  without  localization  will  result  In  the 
bilinear  response  depicted  in  the  left  plot  of  Fig.  9  for  monotonically 
increasing  extension.  Moreover,  the  total  strain  energy  dissipation  for 
tensile  fracture  is 

U  -  —  ah  (2) 

f  E 


which  is  thus  proportional  to  the  height  H  of  the  uniaxial  tension 
specimen.  This  implies  that  in  the  limit  the  energy  dissipation  would 
diminish  to  zero  as  the  height  of  the  specimen  were  reduced  to  zero,  a 
fact  which  has  been  rightly  criticized  by  Bazant  et  al.  at  several 
occasions  [1,2,3].  On  the  other  hand,  if  we  would  adopt  the  notion  of  a 
constant  fracture  energy  release  rate  for  the  formation  of  a  single 
crack,  then  the  energy  expended  for  generating  a  crack  of  unit  area  is 


(3) 


which  should  remain  constant  irrespective  of  the  size  of  the  specimen. 
This  could  be  simply  accomplished  by  adjusting  the  softening  slope 
according  to  the  size  of  the  softening  zone,  i.e.  E-p  ■  -E/3  for  half 
the  specimen  height,  in  order  to  maintain  a  constant  fracture  energy 
release  rate  as  indicated  in  the  right  plot  of  Fig.  9. 


In  order  to  probe  this  notion  of  constant  versus  variable  fracture 
energies  Involved,  let  us  compare  the  two  direct  tension  experiments 
illustrated  in  Figs.  2  and  A.  Converting  the  axial  load  into  nominal 
axial  stress,  the  area  under  the  at~Uf  diagram  can  be  readily 
evaluated,  and  yields  the  following  average  values  for  the  full  and  the 
half  size  NX-  specimens: 


I 

G  *  0.322  lbf/in  ■  56  N/ra  for  H  »  4.25  in  »  108  mm 

f  1 

(4) 

G1  -  0.354  lbf/in  -  62  N/m  for  H  »  2.125  in  *  54  mm 

f 

1/2 


The  fracture  energy  release  rate  slightly  increases  rather  than 
decreases  with  decreasing  specimen  height.  Therefore,  the  direct 
tension  experiments  fully  support  the  fracture  energy  release  rate 
concept.  Thi9  Is  not  surprising  since  the  specimen  ultimately  fails  in 
the  form  of  a  single  localized  raacrodefect  and  not  in  the  form  of 
continuous  debonding  of  raicrodcfects  which  are  uniformly  distributed 
over  the  entire  length  of  the  specimen. 
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In  the  following  discussion  we  will  adopt  this  rationale  to  develop  a 
strain-softening  plasticity  model  which  accounts  for  the  localization  of 
fracture  and  which  monitors  constant  fracture  energy  release  rates  in 
tension  and  shear.  Conceptually,  this  approach  corresponds  to  the 
modification  of  the  softening  modulus  according  to  the  geometry  very 
much  along  the  line  of  the  elementary  example  above. 

Fracture  Based  Strain-Softening  Model  for  Tension 

The  composite  fracture  formulation  is  in  essence  an  attempt  to 
incorporate  constant  fracture  energy  release  rate  in  the 
strain-softening  formulation  of  an  equivalent  elastic-plastic 
continuum. 

For  simplicity  let  us  consider  the  elementary  volume  V  =»  abt  which  is 
subjected  to  the  tensile  state  of  stress  depicted  in  Fig.  10.  The 
maximum  3tress  criterion  results  in  tne  localized  crack  band  At  * 
btt  which  is  oriented  normal  to  the  d.rectio*.  of  major  principal 
stress  "h".  Incipient  cracking  is  controlled  by  the  initially  isotropic 
fracture  condition 

F(c,k)  -  on  -  oc  =  0  where  ■  Oj  with  Oj>o?>03  (5) 

and  ot  -  oc(k) 

For  the  first  active  crack  the  normal  tensile  stress  coincides  with  the 
major  principal  stress  on  *  oj.  The  tensile  strength  ot  describes 
the  degradation  of  bond  stresses  across  the  crack  band  as  a  function  of 
the  tensile  crack  opening  Uf  In  analogy  to  the  fictitious  crack  band 
model  of  Hillerborg,  Modeer  and  Peterson  [7]. 

Restricting  our  attention  to  Mode  I  type  cracking  the  fracture  potential 
Q(e)  coincides  with  the  fracture  condition  F( e,k)  along  the  line  of  an 
associated  flow  rule  and  degenerates  in  the  case  of  primary  crack 
initiation  to  the  unit  vector  in  three  space.  The  associated  fracture 
r»'le 
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defines  ’ireotlcn  and  magnitude  of  the  tensile  strain-sof coning  response 
In  accordance  to  the  flow  rule  of  traditional  plasticity  theory. 

In  view  of  the  strain-space  setting  of  the  computational  algorithm,  the 
consistency  condition 
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leads  ro  the  tensile  fracture  parameter  X, 


where 


and 


«  E  t 


Ef  +  in  E  n 


for  in  Ec  >  0  Progressive  Cracking 


(8) 


for  a\CEc  <  0  Elastic  Behavior 


•  •  • 

where  the  additive  decomposition  of  strain  rate  c  =  ce  +  cf  ^a8  been 
inferred  and  where  E  denotes  the  elastic  material  law  in  three  space. 

Let  us  now  consider  the  hardening  or  rather  the  softening  modulus  Ef 
which  provides  the  key  to  incorporate  the  fracture  energy  concept  within 
the  framework  of  strain-  softening  plasticity.  Formally,  the  softening 
modulus  Ef  in  Eq.  8  is  defined  by 
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In  order  to  control  the  fracture  energy  release  rate  via  the  strain¬ 
softening  modulus  Ef  we  will  pursue  the  following  approach:  Consider 
the  invariance  of  the  softening  branch  for  different  specimen  heights  in 
the  oc-uf  space  and  the  inherent  definition  of  fracture  energy  in 
terras  of 
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Then  it  is  only  natural  to  monitor  softening  directly  in  terras  of  the 
fracture  modulus  E^  which  denotes  the  slope  of  Of-Uf  in  Figure  10. 

The  mapping  between  the  crack  opening  displacement  duf  and  the 
equivalent  tensile  fracture  strain  dk  •  def  leads  to  the  following 
definition  of  the  fracture  energy  based  strain-softening  modulus 

3  a  3uf  3ef  t  ,t  1/2  - - . - 

E,  -  ( - )  <n  where  c,m< c,  t,>  ■  /<  e {>  l+<  e2> ?K  ti>  l  (11) 

*  3uf  3ef  3cf  f  f  t  +  t  *■  f  *  f 


In  Eq .  11  the  MacCauley  brackets  extract  the  tensile  components  of  the 

fracture  strain  rate  tf.  This  provides  the  key  for  the  homogenization 
of  the  localized  crack  discontinuity  over  the  elementary  volume  of  size 
V  -  abt.  In  fact,  if  we  denote  by  ht  the  height  of  the  elementary 
volume  normal  to  the  crack  band  At  *  btt,  then  the  projection  along 
the  normal  vector  "m”  yields  the  geometric  "gage  length"  h  for  the 
evaluation  of  the  softening  modulus  for  cracking  in  direction  "1" 
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As  a  consequence  of  this  definition  of  the  strain-softening  modulus  Ef 
the  resulting  tangential  material  law 


where  (E^  *  E 


E  r»  <n*E 
E^+  n 


(13) 


will  depend  on  the  mesh  size  through  the  geometry  of  the  elementary 
volume  similar  to  non-local  continuum  theories. 


The  composite  fracture  theory  la  based  on  the  equivalence  of  the 
fracture  energy  in  the  crackband  Af  in  the  composite  element  with  that 
of  the  equivalent  continuum  element  V.  This  energy  equivalence  was 
previously  utilized  by  the  authors  (19]  in  order  to  express  the  surface 
dominated  fracture  process  in  terras  of  an  equivalent  volume  dominated 
failure  phenomenon 

J  c^du^dA  -  /  cCdtjdV  (14) 

At  V 

where  the  individual  terras  are  defined  in  Fig.  10.  Assuming  linear 
softening  to  simplify  the  algebra  the  tensile  energy  release  for  the 
formation  of  a  single  crack  band  AA^  within  the  composite  element  is 
s imply 
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The  energy  release  characteristic  of  the  equivalent  continuum  element  is 
defined  by  the  energy  dissipated  during  softening  up  to  rupture 


2  E, 


abt 


where 


(16) 


From  the  equivalence  of  the  energy  release  during  the  formation  of  the 
CE  EQ 

crack  band,  AUt  *  AUt  ,  the  governing  strain-softening  modulus  Ef 
in  the  homogenized  continuum  can  be  expressed  in  terms  of  the  fracture 
modulus  Ej  or  the  equivalent  fracture  energy  release  rate 
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(17) 


In  the  case  of  uniaxial  tension  in  the  axial  direction  the  geometric 
properties  of  the  crack  band  are  readily  expressed  in  terms  of  those  of 
the  eLeraentary  volume,  i.e.  Ac  -  btt  and  ht  ■  a,  thus 
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We  observe  that  this  expression  is  identical  with  the  definition  of  thf 
strain-softening  modulus  Ef  developed  previously  in  Eq .  12.  This 
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result  clearly  demonstrates  that  strain-softening  in  the  equivalent 
continuum  element  depends  on  the  geometry  of  the  elementary  volume  and 
thus  exhibits  a  pronounced  size  effect.  In  other  words,  the  fracture 
modulus  Ef  has  to  be  adjusted  to  the  mesh  size  of  the  particular 

finite  element  Idealization  if  the  fracture  energy  release  rate  G^ 
should  remain  constant.  As  a  result  the  fracture  mechanics  concept 
leads  to  a  non-local  format  of  the  governing  softening  properties  which 
Is  fundamentally  different  from  the  local  (geometry  independent)  format 
of  strain-eof toning  plasticity  and  progressive  fracture  theories. 


In  conclusion  let  us  consider  the  particular  concrete  tested  in  the 
direct  tension  experiment.  For  linear  softening  the  fracture  modulus  is 
readily  evaluated  from  the  nominal  uniaxial  tensile  strength  <j  -  400 
psi  (2.76  MPa)  and  the  fracture  energy  release  rate  for  Mode  I  type 
cracking,  G*  •  0.322  lb/in.  (56  N/m). 
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—  -  -  248.45  biP.  -  -  67.44  — 
in3  mm3 


(19) 


The  characteristic  crack  length  which  defines  the  range  of  debonding 
within  an  elastic  perfectly  brittle  comparison  solid  is  in  this  case 


C*  E 


5.63  in  -  143  am  where  E  *  2800  ksi  =■  19305  MPa.  (20) 


Since  the  characteristic  length  is  larger  than  the  specimen  height, 
linear  elastic  fracture  mechanics  does  hardly  apply.  Along  the  arae 
vein,  previous  suggestions  by  Bazant  et  al  [1]  of  three  times  me 
aggregate  size  for  the  width  of  the  cracl  band  are  also  open  to 
question. 

The  numerical  stability  of  strain-softening  computations  Imposes 
definite  restrictions  on  the  localization  of  fracture  within  a  given 
finite  element  mesh  which  were  recently  discussed  by  the  authors  [20). 

We  recall  that  the  softening  formulation  infers  unstable  material 
behavior  with  a  corresponding  energy  release.  As  long  as  these  material 
instabilities  are  localized  and  contained  within  the  overall  structure, 
stability  In  the  large  can  be  monitored  by  displacement  control  since 
the  structural  stiffness  remains  positive  definite.  On  the  other  hind, 
since  stability  in  the  small  Is  violated  by  the  loss  of  positive 
definiteness  of  the  tangent  material  lav  Ej  in  Eq  .  13,  we  can  examine 
the  limiting  condition  when  (E  j-  *  ->»,  which  corresponds  to  perfectly 
brittle  behavior.  In  the  i_ase  of  tensile  softening  thi  takes  place 
when  the  denominator  in  the  fracture  cjrrectii  i  of  the  elastic  mater. al 
stiffness  vanishes,  i.e.  when  the  following  condition  holds 

Ej  *  rC  E  r  -  0 


(21) 
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Let  U8  consider  for  definiteness  the  case  of  plane  strain  and  the 
normality  vector  n  in  Eq.  6.  Then  the  limiting  condition  for  the 
softening  modulus  reduces  to 


g  <  £  _ ^  ~v _ 

f  ”  ( 1  + v)  ( 1-2  v) 


where  E, 


Ed  ht 


In  other  words  the  mesh  size  ”h"  must  be  such  that  the  following 
constitutive  constraint  is  satisfied 


1  —  v 

(lVy)(l-2y) 


(23) 


For  the  concrete  properties  above  and  v  ■  0.2  this  constraint  results  in 
the  maximum  height  of  the  elementary  volume 

h  <  - - — -  -  12.52  in  -  318.06  ram  (24) 

“  248.45  1.2  *  0.6 

for  which  stable  post-critical  strain  control  may  be  maintained.  We 
observe  that  this  mesh  size  constraint  is  quite  restrictive  if  large 
scale  massive  concrete  structures  such  as  darns  or  pressure  vessels  are 
considered.  It  is  a  direct  consequence  of  the  fracture  energy  basis  of 
the  equivalent  strain-softening  model  which  increases  the  ductility  with 
decreasing  mesh  size  and,  vice  versa,  increases  the  brittleness  with 
increasing  mesh  sizes. 


Strain-Softening  Regime  and  Transition  of  Brittle-Ductile  Failure 


In  the  following  let  us  probe  whether  the  concept  of  a  composite 
fracture  model  for  Mode  I  type  tensile  cracking  could  be  extended  to 
the  Mode  II  type  shear  fracture  in  order  to  capture  the  pronounced 
softening  response  in  low  confined  compression.  In  the  previous 
approach  Day  (4)  and  che  authors  (19)  utilized  the  Mohr-Coulomb 
condition  as  point  of  departure  augmented  by  an  independent  tension 
cut-off.  In  that  case  tensile  debonding  wa3  controlled  by  the 

fracture  energy  release  rate  G^  for  Mode  I  type  cracking  similar  Co 
che  strain-softening  formulation  for  crack  bands  delineated  above. 

On  the  other  hand,  decohesion  due  to  frictional  slip  was  monitored  by 

the  fracture  energy  release  rate  for  Mode  II  type  fracture  governing 

the  formation  of  shear  bands.  Because  of  the  tenfold  increase  of 

uniaxial  compressive  versus  tensile  strength,  the  value  of  G^  was 

approximately  50  times  larger  than  G^  when  faulting  along  a  single 

shear  plane  was  considered  along  the  line  of  the  Mohr-Coulomb  failure 
hypothesis.  At  this  stage  it  was  noted  that  mixed  mode  shear  fracture 
is  actually  initiated  by  tensile  splitting  and  the  formation  of  a  single 
shear  hand  is  thus  the  result  of  coalescing  crack  feathers  along  an 
inclined  slip  plane.  Therefore,  the  postulate  of  an  independent 
faulting  mechanism  is  really  the  consequence  of  the  limited  resolution 
of  currenc  toois  to  trace  each  Individual  tensile  crack  preceding  the 
formation  of  a  single  shear  band. 
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Since  tensile  cracking  and  shear  faulting  are  clearly  interrelated 
we  resort  here  to  a  continuous  failure  condition  which  comprises  both 
tensile  as  well  as  shear  fracture  within  a  single  criterion  according  to 
Hoek  and  Brown  [8].  This  empirical  criterion  was  originally  developed 
to  characterize  the  triaxlal  strength  of  rock  specimens  which  were 
failed  in  the  conventional  triaxial  compression  experiment.  The  threfe 
parameter  failure  model  neglects  rhe  influence  of  the  intermediate 
principal  stress  and  can  be  viewed  as  a  combination  of  the  Mohr-Coulorab 
criterion  with  the  tension  cut-off  condition  very  much  along  the  line  of 
the  parabolic  failure  envelope  originally  proposed  in  1933  by  Leon  for 
concrete  [16]. 

1/2 

F(«,  k)  »  oi  -  o3  -  (cc2  -  mo  a,)  *  0  (25) 

oc  denotes  the  uniaxial  compressive  strength  and  ra  and  c  define  the 
frictional  and  cohesive  strength  material  parameters  to  be  calibrated 
from  the  uniaxial  tension  experiment  and  a  third  condition  when 
"cohesive"  damage  is  considered  with  c<l.  Neglecting  hardening  behavior 
before  the  peak  strength  is  reached  the  uniaxial  compression  experiment 
yields  with  oi  “  0  and  c  ■  1 

2  1/2 

-  Os  *  (c  oz)  where  a  ■  f  *  3.2  ksi  *  22.06  MPa  (26) 

J  c  c  c 

Thus  under  uniaxial  compression  o3  -  -3.2  ksi  as  long  as  intact  cohesive 
behavior  prevails  with  c  ■  1.  The  uniaxial  tension  experiment  on  the 
other  hand  leads  to  calibration  of  m  via  solution  of  the  quadratic 
equation  Eq .  25  with  oj  -  ot  and  o3  =»  0 

2  2 
C°c  '  °t 

m  »  -  where  *  0.9  ksi  -  2.76  MPa  (27) 

a  a 
c  t 

This  condition  yields  ra  ■  7.9  fur  the  frictional  strength  parameter  at 
peak  strength.  Fig.  11  illustrates  the  fit  of  the  failure  envelope  with 
the  experimental  strength  value?  in  terras  of  the  Mohr  diagram  depicting 
the  normalized  shear  strength  vjrsus  effective  normal  stress, 

T  n“ °n ’ 

Let  us  now  examine  whether  the  triaxial  strength  envelope  could  be  used 
to  describe  softening  behavior  in  tension  and  shear,  whereby  both  the 
uniaxial  tensile  strength  and  the  uniaxial  compressive  strength  diminish 
to  zero  in  the  post-peak  regime.  Fig.  12  delimits  the  strain-softening 
regime  in  terras  of  the  maximum  strength  and  the  residual  strength 
envelopes  which  coalesce  at  the  transition  point  P  of  brittle-ductile 
failure.  The  direction  of  inelastic  fracture  strain  rates  has  been 
superimposed  on  a  separate  diagram  in  order  to  illustrate  normality  and 
feasible  formats  of  the  corresponding  plastic  potential. 

If  we  rewrite  the  failure  condition  Eq.  25  then  strain-hardening/ 
softening  can  be  accommodated  by  the  variation  of  the  two  parameters  m 
and  c . 
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F(«,k)  •  x  -  (c  -  mo)  -  0 


where  t 


o1-a3|/ac  and  o“  o^/  <?, 


(28) 


with  ra  =7.9  and  c  “  1 
o  o 


We  recognize  that  m  ■  ra(k)  describes  the  variation  of  the  frictional  and 
c  “  c(k)  that  of  the  cohesive  resistance.  A  possible  interpretation 
considered  previously  in  (4]  assumes  that  softening  is  a  purely 
decohesive  mechanism  where  c(k)  +  0.  This  corresponds  to  a  vertical 
shift  of  the  failure  surface  to  pass  ultimately  through  the  origin. 

This  would  however  imply  that  the  triaxlal  strength  under  high 
confinement  would  also  degrade  while  Fig.  6  clearly  demonstrates  that 
the  opposite  is  true.  In  fact  Fig.  8  illustrates  that  there  is  a  clear 
transition  from  brittle  to  ductile  fracture  between  -0.5  -1.0 

ksi  (-3.45  to  -6.89  MPa)  confining  stress  beyond  which  continuous 
plastic  flow  and  increasing  hardening  tendencies  prevail.  Therefore, 
decohesion  should  be  accompanied  by  simultaneous  frictional  hardening  as 
proposed  by  Vermeer  and  de  Borst  (18)  such  that  the  two  mechanisms 
compensate  each  other  at  the  brittle-ductile  transition  point  P. 


It  has  been  noted  previously  by  Mazars  (12)  and  others  that 
fracture  is  really  induced  by  tensile  strain  rather  than  stress.  Thus 
the  location  of  the  transition  point  P  should  roughly  correspond  to  a 
vanishing  elastic  tensile  st.ain  in  the  radial  direction  if  the 
cylindrical  specimen  is  considered 

er  “  I  (or  "  v  (oe  +  °z})  (29) 

g 

This  yields  the  following  condition  for  c  *  0  and  o  ■  o  : 

r  r  0 

a  “  — —  a  and  o*  ■  8  o,  for  v  ■  0.111  (30) 

r  1-v  2  3  1 


The  tensile  strain  condition  plots  in  Fig.  12  as  a  straight  line  through 
the  origin.  For  a  Poisson's  ratio  of  v  *  0.2  the  tensile  strain  cut-off 
corresponding  to  o3  *  4a^  does  hardly  Intersect  the  maximum  strength 
envelope.  However,  the  value  v  *  0.111  corresponds  to  the  maximum 
strain  criterion  with  o3“-oc+8oi  where  the  ratio  of  the  uniaxial 
compressive  to  uniaxial  tensile  strength  oc/  at*8  coincides  with  the 
biaxial  Griffith  criterion,  see  e.g.  (16).  In  this  case  the  zero 
tensile  strain  condition  in  Eq.  30  does  intersect  the  maximum  strength 
envelope  very  near  the  experimental  transition  point  P  for  brittle  and 
ductile  fracture  which  separates  softening  from  hardening  response 
regimes.  Therefore  we  could  regard  the  tensile  cap  beyond  the  residual 
strength  envelope  as  a  transient  strength  regime  which  ultimately 
collapses  to  the  residual  strength  level  upon  continuous  tensile 
straining.  In  fact  the  trace  of  the  residual  strength  envelope  is 
captured  fairly  well  by  the  following  modification  of  Eq .  28  for  purely 
frictional  resistance 


1/2 

Fr(«,k)  ■  t  —  (-m^a)  “  0 


where  ra  “11.96  and  c  =0 
r  r 


(31) 
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Prom  the  directions  of  the  Inelastic  strain  rates  in  Fig.  12  we  observe 
that  the  trace  of  the  maximum  strength  envelope  furnishes  an  acceptable 
format  for  the  inelastic  potential  Q(  a)  to  be  used  as  fracture  rule  for 
strain-softening  under  mixed  mode  fracture  processes. 

Fracture  Based  Strain-Softening  Model  for  Mixed  Mode  Fracture 

Considering  the  strain  softening  formulation  for  Mode  I  type  fracture  in 
the  previous  section, let  us  examine  the  pertinent  generalization  to 
mixed  mode  fracture  along  the  line  of  the  combined  softening  model  for 
tension  and  shear. 

Starting  from  the  maximum  strength  condition  Eq .  28  we  recognize  that 
decohesion  is  accompanied  by  frictional  hardening.  During  failure  the 
fracture  rule  defines  the  direction  and  magnitude  of  inelastic  fracture 
strains  in  the  homogenized  continuum. 

c.  »  A,  n  where  n  -  —  (32) 

f  f  3c 


From  Fig.  12  we  recognize  that  the  direction  of  inelastic  strain  rate  is 
fairly  well  defined  by  normality  to  maximum  strength  envelope. 

Therefore,  the  fracture  potential  is  in  this  case  described  by  Eq.  28 
with  m  ■  m0  and  c  ■  cQ 

1/2 

Q(c)  “  t  -  (c  -mo)  *  0  where  m-7.9  (33) 

o  o  o 

and  c  ■  1 
o 

The  gradient  of  the  fracture  potential  defines  the  direction  and  thus 
the  inelastic  dilatation  in  principal  stress  space 


3Q 

.  ro 

i  *  1  °  _ 

3oj 

/  .1/2 

2  (c*-m  o) 
o 

3Q 

n  -  —  - 

3Q 

1 

m  ,  . 

0 

3c 

302 

°C 

3Q 

-1 

303 

Similar  to  the  Mohr-Coulomb  condition  there  is  no  inelastic  component  in 
the  direction  of  intermediate  principal  stress.  We  also  note  that  the 
fracture  component  in  the  minor  direction  does  not  depend  on  the 
frictional  resistance  m  because  the  normal  stress  o  in  Eq .  28  is  defined 
only  in  terms  of  the  major  principal  stress  Oj  and  not  in  terms  of  the 
average  value  of  t ne  major  and  minor  normal  stress  used  in  the 
Mohr-Coulomb  condition.  For  computational  reasons  the  fictitious 
cohesion  c*"T2+m0o  has  been  introduced  in  order  to  accommodate 
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"elastic"  trial  stress  states  beyond  the  failure  condition  which  are 
restricted  to  major  principal  stress  values  Oj  < 

The  fracture  parameter  if  determines  the  magnitude  of  inelastic 
straining  during  mixed  mode  fracture.  With  the  aid  of  the  consistency 
condition  F  »  0  the  fracture  parameter  if  reads  in  the  case  of  a 
non-associated  flow  rule  where  F  *  Q 


•  in  it  e  t  • 

-  -  for  m  E  t  >  0 

E^  +  mt  E  ii 


(35) 


Xf  * 


0 


for  tn  Ec  <  0 


In  this  case  the  vectors  ra  and  m  denote  the  gradients  of  the  current 
loading  surface  F(«,k)  and  the  fracture  potential  Q(  c) .  the  latter  of 
which  was  previously  developed  in  terms  of  principal  stresses.  In 
analogy  to  Eq.  34  we  find  with  c*  ■  x2  +  mo 


- 

■ 

3F 

i  + 1 

3  o  ^ 

3F 

1 

»  - 

0 

302 

a 

c 

3F 

3o3 

-1 

(36) 


The  principal  issue  of  the  non-associated  softening  formulation 
revolves  around  the  proper.def inition  of  Ef .  Formally  we  obtain  from 

the  consistency  condition  F  -  0  the  following  expression  for  the 
instantaneous  strain-softening  modulus 


E 


f 


,  3F  3k  *  3Q  .  .  ,  ,  s 

( - )  — -  where  k  -  k(ef) 

3k  3ij  3e  f 


(37) 


We  recall  that  decohesion  is  accompanied  by  frictional  hardening 
both  of  which  are  mobilized  by  tensile  straining  during  compressive 
splitting  in  the  brittle  failure  regime.  Therefore  it  is  nacural  to 

adopt  again  the  definition  of  equivalent  tensile  fracture  strain  e, 
in  Eq .  11  which  played  the  role  of  the  kinematic  damage  metric  of  Ele 
homogenized  continuum  subjected  to  tensile  cracking. 


dk  -  de.  where  e 
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*  Nl/2 


<el>f  +  <e2  >f  +  <C3> 


(38) 


We  note  that  the  MacCauley  brackets  only  extract  the  positive,  i.e.  the 
tensile  fracture  strain  rates  from  the  inelastic  strain  rate  tensor 
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In  case  of  compressive  splitting  the  chain  rule  of  differentiation 
leads  to  the  following  definition  of  the  strain-softening  modulus  in 
terras  of  the  tensile  fracture  behavior 
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3u, 
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— 

(39) 

3c 

3ot 

3uf 

3ef 

3,f 

da 

Recollecting  that  tensile  splitting  in  compression  as  well  as  shear 
faulting  is  preceded  by  distributed  microcracking  and  that  the  actual 
crack  surface  is  in  this  case  a  multiple  of  that  in  the  direct  tension 
experiment,  we  resort  to  the  fundamental  definition  of  Mode  I  type 
cracking.  For  the  current  concrete  this  fracture  property  is 
characterized  by  G*  -  0.322  lbf/in.  (56  N/m)  which  corresponds  to  the 
area  under  the  tensile  softening  law  in  Fig.  13. 

The  individual  terras  in  Eq .  39  are  defined  as  follows 

3F  ”1 

—  - - with  dm  =  -4.06  dc  (40) 

.  ,,  v  1/2 

3c  2(c  -  rac) 

The  interdependence  between  the  cohesive  strength  parameter  and  the 
tensile  strength  results  in  our  case  to 

*  2.5  where  a  -  0.4  ksi  -  2.76  MPa  (41) 

Zr,  t 


The  relationship  between  the  degrading  tensile  strength  ot  and  the 
crack  opening  displacement  uf  was  previously  characterized  by  the 
fracture  modulus  of  the  uniaxial  tension  experiment,  see  Fig.  10. 
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d 


dot  d2G* 
du^  du2 


(42) 


An  exponential  fit  of  the  normalized  test  data  is  shown  in  Fig.  13  which 
leads  to  the  following  decaying  softening  function 
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(4  3) 


where  ct  denotes  the  uniaxial  tensile  strength  and  ur  the  crack 
opening  displacement  at  rupture  when  the  interfacial  bond  stresses  have 
diminisned  to  zero  (in  our  case  ur  =  0.002  in.  -  0.05  mm,  from  Fig.  2). 

T.t  remains  to  relate  the  crack  opening  displacement  uf  to  the 
tensile  fracture  strain  ef  of  the  equivalent  continuum.  In  Fig.  7  we 
have  seen  that  the  lateral  extension  is  over  ten  times  larger  than  the 
axial  extension  in  the  direct  tension  experiment  shown  in  Fig.  2.  This 


should  be  of  no  surprise  since  the  splitting  compression  mode  of  failure 
Is  far  more  distributed  which  Is  In  turn  responsible  for  the  eight  to 
tenfold  increase  of  the  uniaxial  compressive  over  the  tensile  strength 
values . 


A  fracture  energy  argument  similar  to  that  developed  for  tensile 
cracking  leads  the  following  relationship  between  the  crack  opening 
displacement  and  the  equivalent  tensile  fracture  strain  in  compression 


duf  V 

\  =  "C 


(44) 


where  the  volume  fraction  represents  the  microcrack  density  In 

compression.  We  observe  that  Eq .  44  plays  the  fundamental  role  of  the 
geometric  “gage  length"  and  has  the  simple  physical  interpretation  of 
the  crack  spacing  hc  in  the  case  of  splitting  compression.  From  the 
global  observations  it  may  be  evaluated  approximately  in  terms  of  the 

ratio  between  the  fracture  energy  release  rate  in  tension  G,  and  that 

II  r 

in  shear  .  Thereby  it  is  understood  that  the  second  quantity  by 

Itself  is  already  the  result  of  homogenizing  distributed  tensile  cracks 
into  a  single  shear  band  with  the  surface  area  Ag  which  has  the  same 
order  of  magnitude  as  the  crack  surface  in  direction  tension,  l.e.  At 
=  Ag .  From  the  equivalence  of  fracture  energy 

g]  AA  -  G^1  AA„  (45) 

f  c  f  t 


we  can  find  the  unknown  fracture  surface  under  splitting  compression 
f  rora 
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,11 
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(46) 


This  yields  in  the  case  of 
27.2  lbf/in  (4763  N/ra)  the 
spacing 
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II 

the  uniaxial  compression  test  with  Gf  * 
following  measure  of  longitudinal  crack 


h  = 


.1 

jf_ 

,11 


0.05  in  *>  1.28  mm 


(47) 


which  should  be  compared  with  the  characteristic  length  -  5.63  in. 

Ill 

(143  mm)  in  the  case  of  Mode  I  type  cracking.  The  ratio  G^.  / G ^ 
represents  the  magnification  of  tensile  crack  surface  when  compared  with 
that  of  direct  tension  test.  This  ratio  increases  with  increasing 
confining  pressure  in  order  to  simulate  the  diffuse  damage  mechanism 
around  and  beyond  the  brittle-ductile  transition  point.  Using  the 
experimental  data  described  earlier,  a  first  order  approximation  yields 
the  following  expression  for  the  fracture  energy  ratio  as  a  function  of 
the  minor  principal  stress  a 3 
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(48) 


rhis  ratio  reduces  to  one  In  the  case  of  trlaxial  tension  when  <o3>_  = 

0.  It  Increases  rapidly  with  increasing  confinement  o3  <  0  and  reduces 
the  strain-softening  modulus  to  Ef  -*■  0  as  hc  -*-0  which  corresponds 
to  perfectly  plastic  behavior  when  o3  <  0. 


In  conclusion  let  us  summarize  the  tangential 
frictional  hardening  and  non-associated  decohesion 


where 
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material  law  for 
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The  strain-softening  formulation  above  is  based  on  the  fracture 

energy  release  cate  G^  and  the  concomitant  exponential  decay 
function  for  E<j  using  energy  equivalence  concepts  to  homogenize  the 
underlying  failure  phenomena.  Fig.  14  illustrates  the  performance  of 
the  fracture  energy  based  strain-softening  plasticity  model  when  the 
post-peak  response  of  the  uniaxial  compression  test  is  considered.  In 
this  case  only  the  exponential  softening  law  for  tensile  cracking  in 
Fig.  13  was  utilized,  in  order  to  predict  the  softening  response  under 
compression. 


CONCLUDING  REMARKS 


A  strain-softening  elastic-plastic  formulation  has  been  presented 
for  the  description  of  rhe  fracture  mechanisms  in  concrete  which  range 
from  brittle  Mode  I  type  cracking  in  tension  to  ductile  mixed  mode 
fracture  in  compression.  The  composite  fracture  theory  is  based  on  the 
fracture  energy  content  of  the  homogenized  continuum  which  is  subjected 
to  equivalent  strain  softening  mechanisms  in  tension  and  shear. 

We  conclude  that  the  combined  mechanism  of  cohesive  softening  and 
frictional  hardening  covers  the  entire  failure  spectrum  of  concrete 
which  is  monitored  solely  by  the  fracture  modulus  Ej  of  the  underlying 
Mode  I  type  cracking.  We  note  that  an  overall  constitutive  formulation 
for  trlaxial  concrete  behavior  should  include  a  hardening  cap  in 
compression  in  addition  to  the  tensile  strain-softening  mechanism  for 
brittle  failure  in  tension  and  shear. 
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Summary 

The  paper  addresses  computational  issues  of  strain-softening  calculations 
In  so  1 ids •  The  main  focus  is  directed  towards  the  nunerical  simulation  of 
the  localization  of  softening  within  structural  components  with  simultan¬ 
eous  unloading  of  the  elastic  neighborhood.  The  spatial  aspect  of  the 
localization  process  requires  special  provisions  of  the  incremental/ 
iterative  solution  strategy  in  order  to  survive  on  one  hand  local  material 
Instabilities  and  capture  on  the  other  hand  bifurcation  of  the  solution 
path  within  load  Increments  which  are  truly  finite.  The  performance  of 
several  solution  algorithms  will  be  assessed  with  two  simple  model 
problems  which  illustrate  the  strong  interaction  between  the  numerical 
integration  of  strain-softening  constitutive  relations  and  the  spatial 
localization  of  contained  softening  within  elastic  unloading  zones. 

Introduction 

Strain-softening  computations  encounter  puzzling  phenomena  during  the 
numerical  simulation  of  distributed  damage  and  localized  fracture  which 
transcend  previous  experience  with  the  inelastic  analysis  of  "stable" 
materials.  In  fact  these  numerical  difficulties  are  exacerbated  by  the 
loss  of  'normality"  and  the  dominance  of  shear  dilatancy  under  mixed  mode 
fracture  which  results  in  predominant  straining  orthogonal  to  the  loading 
direction.  On  one  side  the  failure  analysis  of  brittle  solids  and  struc¬ 
tures  Is  very  susceptible  to  the  formulation  of  strain-softening  constitu¬ 
tive  relations  which  are  equivalent  in  the  sense  of  strength,  stiffness  as 
well  as  fracture  energy  to  the  actual  discontinuous  failure  process.  On 
the  other  side,  the  numerical  Implementation  of  strain-softening  models  is 
far  from  successful  if  accepted  and  conventional  numerical  tools  for 
strain-hardening  plasticity  are  simply  extended  into  the  softening 
regime.  In  fact  it  is  the  localization  phenomenon  which  leads  to  simul¬ 
taneous  loading  and  unloading  in  adjacent  elements  which  has  to  be  moni¬ 
tored  properly.  Therefore  many  numerical  algorithms  which  perform  very 
well  In  the  case  of  weak  continuous  nonl  inearitles  are  nor.  successful  at 
all  when  a  hard  discontinuity  takes  place  within  a  contiguous  material  In 
a  finite  load  Interval  which  draws  the  demarcation  between  the  localized 
strain-softening  regime  and  the  elastically  unloading  surrounding. 


The  intent  of  this  paper  is  to  clarify  the  computational  issues  of  strain¬ 
softening  calculations.  To  this  end  the  composite  fracture  model  for 
strain-softening  will  be  reviewed  together  with  local  and  global  unique¬ 
ness  conditions  for  non-associated  hardening/softening  solids  in  order  to 
delineate  guidelines  for  successful  strain-softening  computations. 
Subsequently  three  incremental/i terative  solution  methods  will  be 
scrutinized  with  regard  to  their  ability  to  capture  strain  localization. 
Two  simple  model  problems  will  be  adopted  in  order  to  illustrate  the 
respective  performance  of  initial  load  and  tangential  stiffness  methods 
when  applied  to  strain-softening  elastic-plastic  computations.  In 
conclusion  the  constant  arc  length  method  will  be  examined  for  stabilizing 
the  post-critical  response  calculations  beyond  the  range  of  direct 
displacement  control. 

Composite  Fracture  Model  for  Strain-Softening 

Strain-softening  or  rather  the  degradation  of  strength  is  normally  attri¬ 
buted  to  the  change  of  internal  microstructure  of  the  material.  Two 
philosophies  have  recently  emerged  for  the  interpretation  of  the  post- 
critical  response  regimes  of  stroke-controlled  experiments. 

(i)  The  progressive  fracture  (1]  and  continuous  damage  theories  (2) 
start  from  the  notion  of  distributed  raicrodefects  which  coalesce 
continuously  during  progressive  failure.  The  simultaneous  decrease 
of  strength  and  stiffness  essentially  results  in  a  secant-type 
material  formulation  where  no  permanent  deformations  accompany  the 
failure  process  during  loading-unloading  cycles.  This  concept 
corresponds  to  a  parallel  arrangement  of  elastic-per fectly  brittle 
elements  which  fail  progressively  under  Increased  deformation. 

(11)  The  blunt  crack  approach  [3]  and  the  composite  fracture  model  (4)  on 
the  other  hand  start  from  the  notion  of  a  constant  energy  release 
rate  which  accompanies  the  progressive  failure  process  in  the  form 
of  a  contained  crack  band.  The  simultaneous  decrease  of  strength 
and  the  Increase  in  irreversible  fracture  deformations  essentially 
result  in  a  tangent-type  material  formulation  where  no  degradation 
of  stiffness  takes  place  as  long  as  the  elastic  recovery  process 
remains  uncoupled  from  the  progressive  failure  process.  This 
concept  corresponds  to  a  series  arrangement  of  intact  elastic  and 
strain-softening  failure  zones  very  much  along  the  line  of  strain- 
softening  plasticity  which  is  based  on  the  additive  decomposition  of 
deformation  rates  into  elastic  and  plastic  components. 

In  view  of  recent  post-critical  experiments  on  concrete  (5]  the  composite 
fracture  model  will  be  adopted  which  collapses  to  previous  strain- 
softening  formulations  of  e las 1 1 c- pi  as 1 1 c  solids  [6,7]  when  the  fracture 
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energy  release  concept  Is  abandoned.  For  simplicity  we  will  focus  on 
tensile  softening,  i.e.  the  composite  fracture  formulation  of  an  equiva¬ 
lent  strain-softening  continuum  which  maintains  constant  fracture  energy 
release  rate  during  Mode  I  type  cracking.  For  the  extension  to  general 
three-dimensional  conditions  with  Mode  II  and  mixed  mode  failure  the 
reader  is  referred  to  Ref.  (51- 

in  the  case  of  tensile  cracking,  fracture  initiation  is  defined  by  the 
maximum  tensile  9tress  criterion 

F(#,q)  ■  a  -  a  •  0;  where  o  *  Oi  and  oi  >  o5  >  o> 
nt  n  1  1  —  4  —  3 

and  oc  -  oc(q)  (1) 

For  initially  isotropic  conditions  the  normal  tensile  stress  o  in  an 

n 

active  crack  coincides  with  the  major  principal  stress,  i.e.  o^  " 
a j .  The  tensile  strength  ot  describes  the  degradation  of  bond 
stresses  in  terras  of  the  internal  variables  q  which  characterize  the 
change  of  the  Internal  microstructure  during  softening.  In  the  case  of 
single  1  ocallzed  macrodefects  the  fracture  energy  release  rate  G*  will 
provide  the  missing  link  between  the  macroscopic  measure  of  strength  and 
the  Internal  fracture  state.  The  consistency  condition 


•  3FC  .  3F  .  t  .  t  . 

F  *  —  «  ♦  —  q-n  £(t  -  ef)  +  p  q  -  0  (2) 

dC  3q  1 

leads  to  the  following  strain-based  loading  conditions  for  progressive 
cracking 

for  nCE  c  >  0  then  tf  *  0  (3) 

and  q  *  0 

Fig.  1  illustrates  the  progressive  degradation  of  tensile  strength  of 
when  the  elastic  trial  stress  rate  has  a  positive  component  in 

the  direction  of  the  gradient  n  ■  3F/3®.  The  constitutive  behavior  is 
characterized  by  the  following  evolution  laws: 


El astic: 
Fracture: 

Internal  State: 


i  -  E(e  -  ef) 

’  f  .  9Q 

i,  *  X,  m  where  rt  -  — 


(4) 


When  the  elastic  response  behavior  is  not  affected  bv  the  progressive 
fracturing  process  the  elasticity  tensor  remains  Invariant,  F,=c.onst ,  and 
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does  not  depend  on  the  internal  variable  q.  In  this  simple  case  the 
elastic  stiffness  and  the  degradation  of  strength  during  the  progressive 
fracture  process  are  decoupled  entirely  which  is  partly  supported  by  the 
experimental  evidence  reported  in  Ref.  ( 5 1 . 


The  evolution  of  fracture  strains  is  characterized  by  the  fracture 
rule  which  defines  direction  and  magnitude  of  the  macrosco¬ 

pic  strain-softening  response  via  the  gradient  of  the  fracture  potential 
Q*Q(«)  and  the  fracture  parameter  A^.  In  the  case  of  Mode  I  type 
cracking  in  the  direction  of  the  major  principal  stress  the  associated 
fracture  rule  with  tn  *  r  degenerates  to  the  unit  vector  in  three  space 
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The  general  setting  for  the  internal  state  variable  q  follows  the  early 
proposal  by  Kroner  (8]  who  introduced  the  material  tensor  €  in  order  to 
map  changes  of  the  microstructure,  i.e.  the  moments  of  dislocation  pile- 
ups  in  metals  onto  the  macrostructure.  In  our  case  of  tensile  cracking  we 
will  resort  to  this  generalization  of  the  traditional  definitions  of 
strain/and  work  hardening  In  order  to  incorporate  the  concept  of  a 
constant  fracture  energy  release  rate. 


To  this  end  let  us  consider  the  elementary  volume  V-hbt  which  is  composed 
of  an  intact  elastic  portion  and  a  contained  crack  band  of  the  area  A»bt 
where  t  denotes  the  thickness.  Fig.  2  illustrates  the  composite  fracture 
element  and  the  equivalent  continuum  element  which  dissipate  the  same 
fracture  energy  when  the  tensile  strength  diminishes  to  zero  residual 
level  during  increased  straining.  In  the  post-critical  response  regime 
the  tensile  softening  localizes  into  a  single  tensile  crack  band,  the 
formation  of  which  requires  the  following  supply  of  separation  energy 

u 

C*  -  /  rot  duf  (6) 

0 

Note  that  this  expression  is  a  first  approximation  for  the  extraction  of 
the  critical  fracture  energy  release  rate  from  direct  tension  experiments, 
whereby  the  Inelastic  dissipation  mechanisms  which  accompany  the  coales¬ 
cence  of  distributed  microdefects  into  a  single  macrodefect  have  been 
neglected.  If  we  adopt  the  basic  postulate  of  fragmentation  theory  where 
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the  surface  area  generated  during  fracture  is  proportional  to  the  energy 
expended  (9) 

G .  »  -  —  *  const.  (7) 

f  9  A 

then  strain-softening  may  be  monitored  directly  in  terms  of  the  raw 
experimental  data  Involving  the  nominal  tensile  strength  o  and  the 
overall  extension  uf.  Because  of  the  localized  character  of  the  tensile 
crack  and  the  elastic  unloading  process  of  the  remaining  region,  the  over¬ 
all  extension  uj  is  really  the  crack  opening  displacement  in  the  final 
stage  of  separation. 


The  question  remains  how  can  this  relative  displacement  measure  be  conver¬ 
ted  into  an  equivalent  strain  in  order  to  be  amenable  to  strain  rather 
than  displacement-softening  computations?  To  achieve  this  we  resort  to 
the  notion  of  an  equivalent  continuum  element  which  dissipates  the  same 
amount  of  fracture  energy  in  the  o  -e^  space  as  does  the  composite 
fracture  model  in  the  o^-u^  space.  An  elementary  energy  consider¬ 
ation  leads  to  the  following  relationship  between  the  softening  modulus 
Ef  and  the  fracture  modulus  Ej 


E 


f 


dot  d  2cj 

where  E.  *  -  =  — 

d  duf  du‘ 


(8) 


Clearly  in  the  case  of  the  direct  tension  test  the  gage  length  V/A-=h 
corresponds  to  the  height  of  the  specimen  and  thus  to  the  crack  spacing 
which  converts  the  overall  extension  into  a  nominal  strain,  du^.  * 
hdef.  Therefore,  spatial  discretization  with  a  finite  element  mesh 
size  of  that  particular  height  would  yield  the  correct  fracture  energy, 
while  a  finer  mesh  resolution  would  yield  far  more  brittle  response 
predictions,  if  the  fracture  modulus  Ef  were  not  adjusted  to  facilitate 
the  corresponding  element  height.  On  the  other  hand  if  an  elemental 
philosophy  were  adopced  In  which  h  corresponds  to  the  height  or  rather  the 
governing  element  size  in  the  direction  of  the  major  principal  stress  then 
we  have  to  ascertain  that  localized  cracking  is  contained  to  a  single 
element  at  a  time.  In  this  sense  the  fictitious  crack  approach  by  Hlller- 
borg  et  al.  1 1 0 ]  and  the  blunt  crack  approach  by  Bazant  ec  al.  13]  solve 
the  spatial  problem  by  following  a  particular  fracture  topology  similar  to 
the  discrete  fracture  model  of  Ingraffea  and  Saouma  (11).  In  contradic¬ 
tion  the  current  composite  fracture  model  does  not  involve  remeshing  but 
requires  a  solution  strategy  for  the  numerical  integration  of  the  strain- 
softening  relations  which  captures  the  strain  localization  within  a  single 
element  in  order  to  monitor  the  proper  fracture  energy  release  rate. 
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For  the  construction  of  a  general  fracture-based  strain-softening  theorv 
we  have  to  nap  the  crack,  opening  displacement  du^  and  the  equivalent 
tensile  fracture  strain  de^  into  three  dimensions.  To  this  end  let  us 
return  to  the  general  formulation  outlined  in  Eqs.  1-5.  Introducing  the 
evolution  laws  in  Eq .  4  into  the  consistency  condition  in  Eq .  2,  we  obtain 
the  missing  expression  for  the  fracture  parameter 

t  . 

"  — 1 - — -  where  E,  »  -pCG  n  (9) 

r  t  1 

Ef  +  n  E  m 


If  we  resort  to  a  strain  definition  of  the  internal  state  variable 
d<fdij  with  €-1  then  the  metric  of  tensile  deformation  is  defined  by 
the  equivalent  tensile  fracture  strain 


t  1/2  r  , 

€f  -  <  tf  cf  >  -  /  <  +  <  e 


2  >f  +  <  *3 


(10) 


which  extracts  from  the  strain  tensor  only  the  positive  tensile  compo¬ 
nents.  The  fracture  modulus  Ef  in  Eq .  9  characterizes  Che  strain¬ 
softening  behavior  of  the  equivalent  continuum  as  follows 
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where 
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For  uniaxial  tension  this  expression  reduces  to  the  relationship  in  Eq .  8, 
i.e.  Ef»E<jh,  where  h  defines  the  characteristic  dimension  of  the 
elementary  volume  in  the  direction  of  major  tensile  straining.  Note  that 
this  geometric  quantity  lends  itself  to  a  more  fundamental  micromechanical 
Interpretation  since  h  »  V/A  is  nothing  but  the  reciprocal  value  of  the 
mlcrocrack  density,  i.e.  the  crack  surface  per  unit  volume  within  the 
sol  id . 


Combining  the  evolution  laws  for  the  elastic  and  inelastic  fracture 
behavior  we  arrive  at  the  tangential  constitutive  relation  to  be  compared 
with  non-assoclated  elastic-strain  softening  formulations 


ET* 


where  E„ 


E  - 


E  m  r  E 


Ej  +  «  Ecn 


(12) 


The  orinclpal  difference  Is  the  definition  of  the  strain-softening  modulus 
Ef  in  Fq  .  ll  which  monitors  directly  the  fracture  energy  release  rate 

via  F..h 
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On  the  other  hand  the  fracture-based  formulation  reproduces  the  tradition¬ 
al  strain-softening  models  which  describe  the  degradation  of  strength 
directly  In  strain  rather  than  displacement  space. 


Uniqueness 

Let  us  consider  for  a  moment  local  and  global  uniqueness  conditions  of  the 
|  elastic  plastic  strain-softening  solution.  In  the  past  considerable 

'  discussion  has  revolved  around  the  well-posedness  of  strain-softening 

computations  with  repeated  reference  to  Radamard  [12].  The  argument 
relates  to  wave  propagation  in  elastic  solids  whereby  the  hyperbolic 
nature  of  the  characteristic  equation  and  Che  concomitant  speed  of  wave 
propagation  Is  limited  by  the  positive  definiteness  of  the  tangential 
material  law,  l.e. 

>  c  *Cir  with  c  >  C  (14) 


The  local  stability  postulate  of  frrucker  [13]  leads  to  a  similar  conclu- 

i 

I  sion,  since  the  condition 
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iC  c  >  0  (15) 

is  synonymous  with  the  positive  definiteness  of  Ef  according  to  Eq . 

(14).  Clearly  this  stringer;  requirement  leads  to  9evere  restrictions  of 
the  material  model  to  positive  behavior  wlch  unlimited  ductility.  There¬ 
fore,  if  we  intend  to  capture  the  degradation  of  strength  by  strain- 
softening  material  models  we  nave  to  resort  to  contiguous;  uniqueness 
arguments  which  are  considerably  less  restrictive.  In  this  context  we 
recall  the  original  work  of  Melan  [14]  who  studied  uniqueness  of  strain¬ 
hardening  elastic-plastic  solutions.  The  classical  approach  revolves 
around  the  solution  of  an  auxiliary  problem  in  which  we  examine  the 
condition  under  which  two  different  solutions  may  exist  at  the  time  tn+i 
when  a  unique  solution  exists  at  the  time  tn.  For  quasistatlc  condi¬ 
tions  the  proof  of  uniqueness  of  the  solution  of  the  initial  value  problem 
results  in  the  question  whether  the  following  relation  can  vanish  for 
o. ‘i2  and/or  * i 2 

/  ( C;  ~  C])t  (i2  ~  ij)dV  =*  j  Act  it  dV  *  0  (lb) 
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I/)cal  uniqueness  is  therefore  assured  whenever  the  governing  tangential 
material  law  is  positive  definite  i .e.  when 

T 

x  ET  *  >  0  whe re  i  *  4t  ’  («2  '  t|)  (17) 

In  the  case  of  the  non-assoclated  strain-softening  law  in  Eq.  12  this 
constraint  leads  to  the  following  limiting  Inequality  examined  originally 
by  Mroz  (13] 


_  .  tc  .  (xC E  m)  (nCE  x) 

E.  +  n  t  o  >  - 

i  t 

x  E  x 


(18) 


This  leads  in  the  worst  case  to  the  "hard"  restriction  on  the  strain¬ 
softening  modulus  (16]  that 


2E,  > 


(®atE 


,1/2,  t_  ,1/2 

ra)  (n  E  r) 


-  rCE 


rn 


(19) 


Clearly,  for  associated  behavior  with  n  ■  m  we  recover  the  classic  state¬ 
ment  of  He lan  (14]  that  uniqueness  is  guaranteed  for  strain-hardening 
behavior  with  Ef>0. 


On  the  other  hand,  when  the  strain  rate  satisfies  either  orthogonality 
condition 

xC  fc  d  •  0  (20) 

or  xC  E  «  *  0 

then  we  recover  the  far  less  restrictive  uniqueness  condition  for  the 
strain-softening  modulus 

-  E  <  nCE  m  (21) 


Therefore  in  the  most  fortuitous  case  of  orthogonality  between  the  elastic 
stress  rate  and  the  gradients  n  ,m  uniqueness  is  assured  as  long  as  the 

softening  modulus  -Ef  is  smaller  than  r1  Em,  see  also  Ref.  (7],  In 
the  case  of  uniaxial  tension  this  condition  implies  that  the  softening 
slope  must  be  less  than  the  elastic  modulus,  i.e. 

-  E  <  E  (22) 

As  a  result  the  choice  of  softening  modulus  is  quite  restrictive.  Recal¬ 
ling  that  the  fracture  based  softening  formulation  introduces  with  h  a 
characteristic  dimension  into  the  definition  of  Ef  *  E^h  such  that  the 
restriction  in  Eq .  22  results  in  a  following  limitation  of  the  mesh  size  h 
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(23) 


h  <  h 


E 

-E 


This  last  expression  clearly  illustrates  the  interrelationship  between  the 
fracture-based  strain-softening  formulation  and  the  element  size  in  which 
fracture  energy  is  released.  In  fact  the  ratio  between  the  elastic  and 
the  fracture  moduli  governs  the  maximum  mesh  size  hc  for  which  weak 
uniqueness  is  ascertained. 

In  this  context  Valanis  (17]  demonstrated  very  recently  that  rate  depen¬ 
dent  constitutive  relations  for  strain-softening  viscoelastic  and  visco- 
plastic  solids  satisfy  the  local  uniqueness  condition  as  long  as  the 
elastic  component  remains  positive. 

Since  the  pivotal  uniqueness  condition  in  Eq.  16  involves  the  volume 
Integral  over  the  entire  domain,  let  us  examine  the  implication  related  to 
global  uniqueness,  see  also  Ref.  (6].  To  achieve  this  we  subdivide  the 
domain  into  a  zone  V-Vf  in  which  the  positive  material  operator  assures 
local  uniqueness  in  the  strong  sense  and  into  a  zone  Vf  in  which  the 
non-positive  material  operator  may  lead  to  the  loss  of  uniqueness. 


/  it  dv  +  J  A«C  At  dv 


0 


(24) 


V-V 


f 


V. 


Clearly,  the  concomitant  finite  element  discretization  results  In  the 
question,  if  the  tangent  stiffness  remains  positive  definite,  i .e . 


Ar  Kp  at  >  0  where  K^,  *  K  - 

with  K  •  /  BC  E  B  dv 


(25) 


and  Kp  •  /  B  E^  B  dv 


In  the  case  of  contained  strain-softening  computations  the  splitting  of 
the  tangential  stiffness  operator  into  a  positive  elastic  component  K  and 
into  an  non-positive  contribution  Kf  leads  to  the  question:  Under  which 
circumstances  will  the  softening  material  operator  contaminate  the  overall 
stiffness  and  lead  to  a  loss  of  positlvenass  of  Clearly,  theie  is 

no  general  answer  to  this  protlem,  because  the  topology  and  boundary 
conditions  determine  the  extent  of  the  localized  strain-softening  region 
Vf  for  a  given  load  history.  As  a  guideline  one  might  speculate  that 
sharp  localization  of  the  softening  zone  will  enhance  the  confinement  and 
will  thus  retard  the  loss  of  post t iveness .  Therefore,  the  numerical 
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integration  of  the  differential  constitutive  relations  has  to  capture  the 
spatial  localization  of  strain-softening  as  well  as  contiguous  elastic 
unloading  of  the  neighborhood. 
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In  conclusion  let  us  consider  global  stability  which  corresponds  to  the 
local  statement  in  Eq.  15.  In  this  case  the  change  of  stress  leads  to  the 
following  exchange  of  internal  and  external  energy 


i  ;  tV* '-i  i' 


v-v. 


t  dv  -  i  (  Pg  w  dS  ♦  I  /  pCus  ds 

2  s  1  s 

o  u 


(26) 


Adsorpt  ion 


Re  lease 


Positive 


Positive/Negative 


where  body  forces  are  neglected  and  a  small  displacement  gradient  theory 
is  presianed.  The  softening  in  region  Vf  results  in  a  release  of  energy 
which  is  monitored  by  the  fracture-based  softening  law  in  Eq .  12.  On  the 
other  hand  the  rate  of  external  work  by  the  prescribed  surface  tractions 
p8  is  always  positive,  while  that  of  prescribed  surface  displacements 
us  can  be  positive  or  negative.  This  clearly  demonstrates  the  stabiliz¬ 
ing  Influence  of  uisplacement  control  which  provides  the  computational 
tool  for  tracing  the  post-critical  response  behavior.  On  the  finite 
element  level  Eq .  26  reduces  to  the  incremental  equilibrium  relation 


(K  -  Kf )  r 


R 

P 


+  R 

u 


(27) 


If  we  transform  the  non-positive  softening  stiffness  Kf  into  an  equiva¬ 
lent  force  contribution  on  the  right  hand  side  we  obtain 


K  r  -  fi,  f  R  +  R  (28) 

f  p  u 

where  Ff  denotes  the  crack  driving  force  defined  as 

Rf  -  Kfr  •  /  BeE  t{  dV  (29) 


Note  that  the  fracture  strain  rate  i ^  monitors  crack  extension  via  the 
fracture  energy  release  rate  embedded  in  Eqs.  9  and  11.  We  recognize  that 
if  ■  E  cf  in  expression  Eq .  29  corresponds  to  the  crack  extension  stress 
In  fracture  mechanics  which  monitors  the  separation  of  nodes  between  con¬ 
tiguous  elements  when  a  discrete  fracture  approach  Is  adopted  (10,11). 

The  principal  difference  between  this  scheme  and  the  discrete  fracture 
approach  Is  chat  equivalent  strain-softening  takes  place  within  an  element 
rather  than  at  the  Interface  between  adjacent  elements.  This  fracture 
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energy  equivalence  is  assured  as  long  as  strain-softening  is  localized  in 
a  single  element  and  does  not  spill  over  into  contiguous  elements. 

Model  Problems 

In  conclusion  we  examine  the  stability  and  uniqueness  of  the  numerical 
solution  with  the  two  elementary  model  problems  shown  in  Fig.  3.  The 
series  arrangement  of  two  axial  force  members  of  different  cross-section 
was  selected  in  order  to  investigate  the  performance  of  three  different 
incremental/iteratlve  solution  methods  for  the  nuoerical  integration  of 
the  strain-softening  relations.  Of  particular  Interest  was  the  question 
whether  the  computational  strategies  were  able  to  capture  the  localization 
of  fracture  in  the  weaker  element,  and  if  a  unique  solution  could  be 
obtained  independently  of  the  size  of  the  load  increment.  Of  further 
interest  was  the  overall  computational  efficiency  and  the  limiting  condi¬ 
tions  for  convergence  in  the  light  of  the  theoretical  considerations 
delineated  in  the  previous  two  sections.  Particular  care  was  exercised  to 
implement  the  three  algorithms  in  their  full  generality  and  not  to  bias 
che  solution  effort  by  taking  advantage  of  the  simple  structural  configur¬ 
ation  for  which  hand  calculations  suffice  to  determine  the  reference 
solution . 

The  odd  geometry  and  material  parameters  were  chosen  to  fit  the  direct 
tension  experiments  on  cylindrical  concrete  specimens  reported  in  Ref. 

(3).  A  linear  soften.ng  model  was  adopted  to  monitor  the  fracture  energy 
release  rate  Gf  ■  0.322  lb/in.  which  was  required  to  separate  the 
tension  specimen  into  two  competent  parts.  For  the  height  h  •  2.125  in. 
of  each  bar  element  the  equivalent  strain-softening  modulus  Ef  is 
obtained  as  follows: 

, 

E,  ■  E.h  •  i  _  h  «  0.530  x  10b  psi  (30) 

f  d  2  rl 
Gf 

with  *  40C  psi 

and  E  ■  2.8  x  10b  psi 

Figure  4  illustrates  the  element  deformations  at  three  stages  under  the 
prescribed  end  displacement  rp>0.  The  results  clearly  illustrate  that 
beyond  peak  softening  leads  to  a  rapidly  increasing  extension  of  the  weax 
element  while  the  stronger  elements  contract  due  to  elastic  unloading. 
Therefore  severe  strain  localization  takes  place  during  monotoni cal ly 
Increasing  loading  of  the  structure  under  displacement  control  which 
results  In  a  discontinuous  strain  distribution  along  the  axis,  whereby  the 
deformation  rates  in  the  softening  and  the  elastically  unloading  element 
nave  opposite  directions. 
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Fig.  4  Localization  of  Strain  -  Softening  Regime 


Figure  5  shows  the  representative  behavior  of  the  two  bar  structure  when 
the  continuation  process  is  resolved  into  six  relatively  fine  increments. 
The  iterative  refinement  within  each  increment  demonstrates  rapid  conver¬ 
gence  of  the  strain-softening  and  elastic  unloading  deformations  in  the 
two  elements.  The  accompanying  relaxation  of  "rear'  stresses  is  plotted 
when  an  elastic  predictor  is  used  along  the  line  of  the  initial  load 
iteration  method  which  dates  back  to  the  original  work  in  Refs.  (18,19). 

The  algorithm  is  summarized  in  the  Appendix  which  corresponds  to  the 
elastic  plastic  solution  strategy  developed  in  Ref.  (20).  We  observe  that 
the  initial  elastic  predictor  step  at  the  beginning  of  each  increment  is 
iteratively  refined  until  the  constitutive  constraint  and  equilibrium  are 
satisfied  simultaneously.  It  is  intriguing  chat  the  associated  inelastic 
flow  corrections  lead  to  a  stress  relaxation  which  is  particularly 
pronounced  in  the  stronger  element  B  subjected  to  elastic  unloading. 

In  the  previous  study  the  incrementation  was  chosen  in  such  a  way  that  the 
initial  predictor  step  was  causing  strain-softening  in  the  weaker  element 
only.  In  this  case  all  three  algorithms  summarized  in  the  Appendix 
captured  the  localization  of  strain-softening  in  element  A,  whereby  the 
tangent  modulus-return  scheme  was  most  efficient  because  of  its  rapid 
convergence  within  two  Iterations.  However,  this  conclusion  only  holds 
for  the  refined  staging  strategy  in  which  incipient  softening  is  traced 
for  one  element  at  a  time.  Therefore  in  a  complex  structure  the  continua¬ 
tion  process  would  have  to  be  subdivided  into  extremely  fine  increments  in 
order  to  capture  separate  localization  events  in  each  softening  element. 
This  procedure  corresponds  to  the  original  proposals  of  elastic-plastic 
computations  at  the  end  of  the  1960's  in  which  the  load  incrementation  was 
adapted  to  Incipient  yielding  of  the  most  critical  element  within  the 
mesh.  In  computational  practice  we  have  to  resort  to  much  bolder  staging 
strategies  which  are  independent  of  the  different  localization  points  in 
the  continuation  process.  In  fact  the  determination  of  the  next  localiza¬ 
tion  point  is  in  general  quite  costly  particularly  when  arc-control  is 
used.  Therefore  a  more  pragmatic  approach  resorts  to  the  following 
strategy  for  a  prescribed  load  increment  : 

(1)  determine  from  the  elastic  predictor  step  Che  most  critical 
trial  stress  state  within  the  structure 

(il)  monitor  the  structural  response  due  to  softening  within  this 

critical  element,  while  preventing  all  other  elements  from  strain 
1 ocal 1 zat  ton. 

This  strategy  was  certainly  successful  when  the  present  model  problems 
were  analyzed  with  large  load  increments,  which  might  otherwise  have 
resulted  in  the  simultaneous  softening  of  all  elements.  Because  of  the 
extensive  search  and  sorting  effort  this  approach  Is  equally  costly  since 
it  involves  continuous  exchange  between  structure  and  element  oriented 
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oriented  data  within  large  scale  finite  element  programs.  Therefore  it  is 
highly  desirable  that  the  solution  algorithm  capture  localization 
automatically  wichout  resorting  to  the  artifacts  above. 

To  this  end  we  return  to  our  two  model  problems  and  examine  the  perform¬ 
ance  of  the  three  solution  algorithms  summarized  in  the  Appendix.  The 
load  is  applied  in  two  large  increments;  the  first  is  purely  elastic  while 
the  second  leads  to  simultaneous  softening  of  all  elements.  Fig.  6  de¬ 
picts  the  localization  of  deformations  within  the  structure  as  predicted 
by  the  initial  load  method.  Observe  the  snap  back  and  the  elastic  recov¬ 
ery  of  the  stronger  elements  within  the  second  iteration  of  the  second 
Increment.  The  initial  load  method  captures  the  strain  localization  in 
spite  of  the  fact  that  the  initial  elastic  predictor  ste,  mobilizes 
strain-softening  in  all  elements.  In  fact  this  snap  back  in  the  stronger 
elements  is  captured  by  two  essential  ingredients  of  the  algorithm 

(1)  the  initial  elastic  stiffness  is  utilized  throughout  the  analysis 
along  he  line  of  Sq .  28. 

(li)  during  elastic* Inelastic  transition  the  penetration  point  is 
located  through  the  contact  parameter,  £. 

In  fact  it  is  t l.o  latter  argument  which  results  in  elastic  snap  back  of 
the  stronger  elemencs  because  of  predominant  softening  in  the  weak 
element.  In  the  case  of  the  two  bar  model  problem  the  initial  load  method 
determines  the  increment  of  the  free  displacement  degree  of  freedom  Ar  j 
according  to  the  following  Iteration  scheme 

<k  +  k_)  Ar}  +  1  -  arJ  ♦  AR  (31) 

A  B  *  f  u 

where  kA  and  kg  denote  the  elastic  stiffnesses  of  elements  A  and  B. 

The  crack  driving  force  due  to  strain-softening  in  the  two  elements  is 
given  by 

AR*  -  (AX*  AE).  ♦  (AX*  AE)_  (32) 

r  r  A  id 

where  the  fracture  parameter  for  uniaxial  tension  is 

Aa1.  -  (1  -  £l)  — -  AC  -  1.232  (1  -  C1)  Ac  (33) 

f  Ef  ♦  E 

Tnc  contact  parameter  f  locates  che  penetration  point  for 
elastic-inelastic  transition  according  to  linear  interpolation 
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a  -  a 

y  n 


o  .  -  o 
n+1  n 


where  0  <  5  £  1 

for  o  .  >  o 
n+1  y 


o  denotes  the  initial  value  of  the  tensile  strength  o  -  o  (e,)  and 
y  y  t  f 

o  and  o  .  are  the  initial  stress  and  the  current  value  of  the 
n  n+ 1 

trial  stress  state.  We  observe  that  £  -  0  for  continuous  loading  and 

reduces  to  £«1  if  o*  .  <  ov.  The  strain  increment  Ac  In  Eo.  33  can 

n+1  —  y 

be  expressed  in  terms  of  the  displacement  increment  which  yields  for  the 
crack  extension  force 

ArJ  -  1.232  kA(l  -  -  1.232  -  cJ)(Arp  "  Arl>  (35: 

where  Ar2  *  Arp  denotes  the  prescribed  displacement  increment.  In  the 

case  when  the  weak  element  has  arrived  at  the  elastic  limit  £  ■  0  and 

A  0 

when  the  sufficiently  small  load  increment  is  applied  such  chat  "  1 
then  the  iterative  corrector  equation  (Eq  .  31)  reduces  to 


(k  +  k  )  Ar1/1 


1.232  k  Art  +  k  Ar 


Convergence  of  the  initial  load  iteration  is  assured  since 


1.232  k, 


k  +  k 
A  B 


0. 58  <  1 


We  observe  that  strain-softening  results  in  overrelaxation  of  the  fracture 
parameter  AX^  in  Eq.  33  because  of  the  magnification  factor  E/(Ef  +  E)  * 

1.232  >  1.  Therefore  we  could  ask  the  question  what  is  the  limiting 

£ 

softening  modulus  Ef  for  which  convergence  breaks  down  when 


- -  1 


Ef  +  E  kA  +  kB 


l 

The  convergence  limit  is  reached  when  E-  -  E,  ■  -1423  ksi.  For  the 

I  1  r 

fracture  energy  release  rate  of  G^-0.322  lb/in.  this  limitation  results 
in  the  maximum  element  height  of 

h  -  5.47  in  (39) 

c 

below  which  convergence  can  be  assured.  Note  that  this  global  convergence 
limitation  is  considerably  more  restrictive  than  the  weak  uniqueness  con- 
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ditlon  in  Eq  .  22  which  was  based  only  on  the  local  material  consideration 

-Ef  <  E  -  2800  ksl  (40) 

Clearly  these  elementary  convergence  observations  are  complicated  by  the 
general  transition  problem  which  has  to  be  considered  for  large  load 
Increments,  when  both  elements  extend  into  the  softening  regime.  However, 
from  the  higher  stress  state  in  the  weak  element  A  we  conclude  that 
<  £  which  signals  predominant  softening  in  clement  A  and  results  in 
the  snap-back  of  the  stronger  element  B  into  the  elastic  recovery  regime 
according  to  Eq.  35. 

Let  us  now  examine  the  performance  of  the  modified  tangent  modulus  return 
scheme  summarized  in  the  Appendix.  This  particular  strategy  was  developed 
by  Vermeer  (16]  for  the  solution  of  non-associated  shear  dilatancy  prob¬ 
lems  and  it  perforaed  best  among  the  numerous  variations  of  tangential 
methods  which  were  recently  assessed  in  Ref.  (21).  This  strategy  essen¬ 
tially  consists  of  three  parts  an  elastic  predictor  step  and  for  plastic 
loading  a  subsequent  plastic  corrector  step  which  is  combined  with  the 
out-of-balance  equilibrium  step  depicted  in  the  flow  chart  of  the  Appen¬ 
dix.  Unfortunately,  many  recent  publications  focus  on  a  particular  aspect 
of  the  algorithm  such  as  the  plastic  corrector  step  and  leave  the  Import¬ 
ant  details  of  initial  staging  and  synthesis  with  the  out-of-balance 
equilibrium  correction  to  the  reader.  In  any  case  the  plastic  correction 
Involves  in  general  the  solution  of  elastic-plastic  transition  problem  for 
the  contact  point  denoted  by  the  subscript  "c"  and  the  tangential  computa¬ 
tion  of  the  intermediate  state  "I”  which  is  subsequently  returned  to  the 
current  yield  surface  by  the  plastic  flow  corrector.  The  essential  step 
is  in  this  case  the  calculation  of  from  the  incremental  format  of 

the  consistency  condition  using  a  Taylor  expansion  about  the  intermediate 
configuration  and  assuming  that  the  total  strain  Increment  during  flow 
correction  remains  constant.  This  yields  for  the  plastic  parameter  the 
following  expression 


<  F  («j^  o  )  > 

AX  .  - - - 

Ej  ♦  E  m 


(41) 


when  the  intermediate  trial  stress  state  exceeds  the  current  yield  condi¬ 
tion.  In  contrast  to  the  initial  load  method  which  returns  after  each 
Iteration  to  the  Initial  stage  "n"  the  state  variables  are  now  continuous¬ 
ly  updated  In  the  spirit  of  the  Newton-Raphson  iteration  strategy. 
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In  the  case  of  our  single  degree  of  freedom  model  problem  the  tanecnt 
modulus  return  scheme  results  In  tha  following  Iterative  corrector  of  the 
incremental  equilibrium  relation 

(kA  +  kB}T  =*  R  -  si  +  +  (42) 

For  uniaxial  conditions  and  linear  hardenir.g/sof tening  laws  the 

intermediate  predictor  state  will  remain  on  the  yield  surface,  thus  ARf 

■  0.  Simillarly  if  we  monitor  the  load  history  such  that  only  the  weaker 

c 

element  is  subjected  to  Incipient  softening  ARf  *  0  and  the  tangent 
stiffness  is  confined  to  element  A.  Since  there  are  no  external  loads 
applied  at  the  free  degree  of  freedom,  Eq.  42  reduces  in  the  case  of  the 
fine  staging  scrategy  to 

( k. a  T  +  kg)  Ar[+l  -  -(Ao1)A  +  (AaSg  (43) 

where  o*  ■  o.  +  E_  Ac*  -  AX*E 
A  A  T  A  f 

and  cjg  -  og  ♦  E  A  e* 

where  the  righthand  side  represents  the  out-of-balance  force  of  the  free 
degree  of  freedom.  For  the  fine  staging  strategy  the  tangent  modulus 
return  scheme  captures  immediately  softening  of  the  weak  element  A  and 
elastic  unloading  of  the  stronger  element  B  resulting  in  a  convergent 
solution  after  one  iteration. 

However,  in  the  case  of  the  coarse  staging  strategy,  where  both  elements 

mobilize  softening  the  tangent  modulus  return  scheme  results  in  divergent 

oscillations  since  (kA  +  kg>x  is  now  negative  definite.  As  a  result 

the  tangent  modulus  return  scheme  breaks  down  prematurely  because  of  its 

inability  to  snap  back.  Only  in  the  case  of  extremely  careful  staging  of 

the  load  increments  the  tangent  approach  captures  the  strain  localization 

within  the  structure.  Ultimately  failure  will  occur  when  the  tangent 

stiffness  looses  positive  definiteness  which  leads  to  the  same  restriction 

£ 

of  the  softening  modulus  £f  ■  -1.423  ksl  as  the  initial  load  method. 

The  third  method,  the  initial  modulus  return  scheme  in  the  Appendix  is 
very  similar  to  the  tangent  modulus  return  method  except  for  the  initial 
stiffness  which  is  retained  throughout  the  iteration  process.  The  plastic 
flow  corrector  step  has  been  discussed  extensively  in  Ref.  (22)  and 
corresponds  essentially  to  the  flow  correction  from  the  intermediate 
state,  Eq .  41.  The  performance  of  this  algorithm  is  very  similar  to  the 
initial  load  method.  For  fine  staging  convergence  is  obtained  in 
approximately  the  same  number  of  iterations.  For  coarse  staging  snap  back 
is  oscillatory  and  not  as  prompt  as  in  the  initial  load  method.  However, 
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the  same  strain  localization  is  attained  in  the  weak  element  when  cwo 
large  load  increments  are  used. 
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We  conclude  that  the  localization  of  strain-softening  computations  imposes 
severe  restrictions  on  the  staging  strategy  when  tangential  approaches  are 
used.  The  initial  stiffness  methods  are  In  this  regard  less  restrictive 
and  capture  strain  localization  irrespectively  of  the  magnitude  of  the 
load  Increments.  All  methods,  however,  break  down  when  the  tangent 
stiffness  looses  positive  definiteness.  Therefore  an  in-depth  effort  Is 
presently  on  the  way  to  explore  the  capabilities  of  the  arc-length  methods 
(23, 24]  to  capture  strain  localization.  Initial  results  for  elastic 
softening  material  laws  indicate  very  promising  results  when  the  modified 
constant  arc  length  normal-projection  method  (23]  is  extended  by 
orthogonality  constraints  between  the  initial  displacement  predictor  and 
the  Iterative  displacement  corrections  using  the  Initial  stiffness 
throughout  the  iteration  history.  In  this  case  stable  strain  localization 
was  obtained  for  our  model  problems  using  strain-softening  moduli  Ef  < 
-1.423  ksl  which  fully  contaminated  the  positive  definiteness  of  the 
tangent  stiffness  operator.  These  results  will  be  detailed  in  a 
forthcoming  publication. 
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ABSTRACT 

The  composite  fracture  model  will  be  delineated  which  monitors  constant 
fracture  energy  release  rates  1'n  tension  within  an  equivalent  strain 
softening  plasticity  formulation.  Its  predictive  value  will  be  Illustrated 
with  the  one-dimensional  model  problem  of  a  uniaxial  tension  member  which 
has  been  previously  used  to  query  the  well-posedness  of  quasistatic  and 
dynamic  strain  softening  computations. 

INTRODUCTION 

It  Is  widely  accepted  that  post-critical  experiments  under  stroke  control 
are  not  only  accompanied  by  changes  of  the  Internal  mlcrostructure  but  also 
by  changes  of  the  macrostructure  In  the  form  of  localized  fracture  zones 
which  penetrate  the  test  specimen  of  finite  dimension.  Therefore  strain 
softening  requires  homogenization  of  the  non-uniform  deformation  process  In 
order  to  extract  ductility  properties  from  servo-controlled  experiments 
which  do  not  only  hold  for  a  specific  test  configuration.  In  the  case  of 
plain  concrete  the  underlying  Issue  of  a  single  macrodefect  versus  distribu¬ 
ted  microdefects  Is  further  complicated  by  the  conposlte  nature  of  the 
mortar-aggregate  mixture.  Clearly,  the  tortuous  character  of  the  macro¬ 
defects  and  the  progressive  localization  of  distributed  microdefects  cast 
considerable  doubt  on  the  validity  of  fracture  mechanics  on  one  hand  and  on 
the  relevance  of  continuous  damage  mechanics  on  the  other  hand.  The  princi¬ 
pal  question  centers  around  the  issue  If  progressive  failure  Is  a  surface  or 
a  volume  dominated  process.  In  direct  tension  there  Is  little  argjment  that 
the  formation  of  a  discrete  macrocrack  poses  a  fracture  mechanics  problem  In 
which  the  energy  expended  for  the  crack  formation  Is  approximately  propor¬ 
tional  to  the  surface  area  generated,  l.e.  Gf  ■  constant.  However,  In  the 
case  of  splitting  conpresslon  and  even  more  In  the  case  of  faulting  modes  of 
failure  this  problem  Is  far  from  resolved  and  has  been  recently  forwarded  In 
the  form  of  the  provocative  question,  If  the  formation  of  "shear  bands  Is 
sheer  nonsense"  in  the  case  of  brittle  materials. 
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In  the  subsequent  development  the  conposlte  fracture  model  will  be  brief¬ 
ly  outlined  which  Incorporates  the  constant  fracture  energy  release  rate 
concept  Into  the  strain  softening  formulation  of  an  equivalent  elastic 
plastic  solid.  The  general  formulation  will  be  subsequently  applied  to  the 
solution  of  a  one-dimensional  model  problem  In  which  a  uniaxial  specimen  Is 
subjected  to  increasing  tension.  Because  of  tensile  softening,  strain 
localization  takes  place  In  a  transverse  crack  band  wnerever  the  stress 
amplitude  first  exceeds  the  governing  tensile  strength.  The  trapping  of 
tensile  deformations  appears  Irrespective  of  the  underlying  strain  softening 
mechanics  which  may  resort  to  singular  fracture  or  distributed  damage  modes 
of  failure.  The  strain  localization  slnply  reflects  the  continuum  approxi¬ 
mation  of  the  actual  crack  discontinuity  and  the  concomitant  separation  of 
the  specimen  Into  two  unloading  components. 

In  the  recent  past  considerable  discussion  has  centered  aoout  the  well- 
posedness  of  quasistatic  and  dynamic  strain  softening  computations  [1],  the 
uniqueness  of  rate-sensitive  material  formulations  [2],  the  advantages  of 
viscoelastic  and/or  viscoplastic  strain  softening  models  [3],  and  alterna¬ 
tively  the  need  for  non-local  continuum  theories  along  the  line  of  the 
Imbricate  fornulatlon  of^long  range  deformation  measures  [*.5].  Some  of 
these  arguments  will  be  addressed  In  the  context  of  the  elementary  model 
problem  for  which  the  composite  fracture  formulation  provides  both  a  unique 
and  mesh-insensltl ve  solution  of  the  fracture  energy  [6],  In  fact,  we  will 
observe  that  the  wave  propagation  problem  Is  conputatlonal ly  less  demanding 
than  the  corresponding  quasistatic  strain  softening  analysis  of  contlnua. 

COMPOSITE  FRACTURE  MODEL 

The  composite  fracture  fornulatlon  was  recently  developed  by  the  authors 
In  order  to  Introduce  the  constant  fracture  energy  release  rate  concept  Into 
strain  softening  confutations  of  an  energy-equivalent  continuum.  The  ori¬ 
ginal  approach  In  ref.  [7]  started  from  two  Independent  softening  mechanisms 
for  Mode  I  type  debonding  in  tension  and  for  Mode  II  type  decohesion  In 
shear  considering  some  limited  Interaction  [18],  8ased  on  new  experimental 
evidence  on  low  strength  concrete  which  was  tested  under  stroke  control  In 
direct  tension  and  low  confined  compression,  a  fracture  energy  based  strain 
softening  plasticity  model  was  subsequently  developed  In  ref.  [8].  In  this 
case  decohesion  under  low  confined  compression  was  expressed  In  terms  of 
Mode  I  type  fracture  which  was  monitored  In  terms  of  an  equivalent  tensile 
rather  than  a  shear  strain  measure.  The  point  of  departure  for  the  fracture 
based  strain  softening  plasticity  model  was  the  Leon  criterion  of  trlaxlal 
strength  [9]  which  corresponds  to  a  parabolic  general  1 zatlon  of  the  Mchr- 
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Coulomb  criterion  augmented  by  the  appropriate  tension  cut-off  conditions. 
This  formulation  resorts  to  a  single  tensile  strain  controlled  softening 
mechanism  which  describes  splitting  conpression  and  shear  faulting  In  terms 
of  the  fundamental  Mode  I  type  fracture  process  [19]. 

In  view  of  the  nature  of  the  elementary  model  problem  below  we  will  focus 
on  the  composite  fracture  model  for  Mode  I  type  cracking  In  tension.  In 
this  case  the  formulation  corresponds  to  the  strain  softening  plasticity 
model  of  Nilsson  and  Oldenourg  [10]  when  the  fracture  condition  degenerates 
to  the  maximum  stress  criterion  of  Ranklne: 

F(#,q)  •  on  -  0t(q)  •  0  where  op  •  ox  (1) 

and  jj  >  9j  ^  o3 

For  Isotropic  conditions  the  normal  tensile  stress  on  coincides  with  the 
major  principal  stress.  I.e.  on  ■  ox.  For  permanent  crack  memory  sub¬ 
sequent  cracking  Is  monitored  separately  along  the  three  directions  defined 
oy  the  first  crack  occurrence.  The  tensile  strength  parameter  describes 
the  degradation  of  bond  stresses  In  terms  of  the  Internal  variables  qx 
which  characterize  the  tensile  fracture  deformations  e|  along  the  ‘*ree 
crack  directions.  This  "orthotroplc"  softening  strategy  deviates  sightly 
from  the  Isotropic  formulation  used  previously  [8]  In  which  an  equlvilent 
tensile  strain  cf  was  extracted  from  the  Inelastic  strain  tensor  which 
monitors  the  degradation  of  tensile  strength.  For  rotating  stress  direc¬ 
tions  the  orthotroplc  approach  requires  however  additional  provisions  along 
the  line  of  the  rotating  crack  model  of  Schnobrlch  et  al.  [11]. 

In  the  case  of  localized  macrodefects  the  fracture  energy  release  rate 
Gf  describes  the  degradation  of  tensile  strength  ot  In  ten*  of  the 
crack  opening  displacement  u*.  Therefore  the  fracture  energy  formulation 
of  strain  softening  has  to  establish  a  relationship  between  the  crack 
opening  displacement  uy  and  the  nominal  tensile  strain  measure  cf  In  the 
equivalent  continuum.  To  this  end  let  us  consider  the  elementary  volume 
ve  »  hebt  depicted  In  Fig.  1  which  Is  composed  of  an  Intact  elastic 
region  and  a  contained  crack  band  of  the  area  Af«  bt,  where  t  denotes  the 
thickness.  Fracture  of  the  conposlte  fracture  element  should  release  the 
same  amount  of  fracture  energy  as  that  of  the  equivalent  continuum  element 
when  the  crack  band  propagates  through  the  elementary  volume.  A  simple 
energy  consideration  leads  to  the  following  definition  of  the  strain  soft¬ 
ening  modulus  If  of  the  equivalent  continuum  In  terms  of  the  fracture 
modulus  £(j  and  the  underlying  fracture  energy  release  rate 
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Clearly  In  the  case  of  a  uniaxial  tension  specimen  the  characteristic  length  j 

Ve/Af  corresponds  to  the  height  of  the  elementary  volume  hfi  which 
converts  the  crack  opening  displacement  Into  the  overall  extension  of  the 
composite  fracture  element  after  elastic  unloading.  Therefore  the  charac¬ 
teristic  length  plays  the  role  of  the  gage  length  which  calibrates  the 
"nominal "  fracture  strain  In  the  form 
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A  more  fundamental  viewpoint  regards  the  volume  fraction  Af/Ve  as 
microcrack  density  which  maps  the  fracture  surface  area  Into  an  equivalent 
volume  quantity.  Note  that  Eq.  3  Is  the  key  to  the  homogenization  of 
localized  macrodefects  Into  a  fracture  energy  equivalent  continuum. 

In  a  conputatlonal  environment  the  characteristic  length  Ve/Af 
corresponds  to  a  representative  dimension  of  the  mesh  size,  l.e.,  the  cubic 
root  of  the  Jacobian  If  an  Isoparametric  finite  element  formulation  Is 
utilized.  This  elemental  approach  ensures  that  the  proper  amount  of 
fracture  energy  Is  released  when  a  single  element  Is  subjected  to  complete 
softening.  Clearly  the  localization  of  strain  softening  within  the 
structural  Idealization  Is  a  separate  Issue  which  requires  computational 
provisions  In  order  to  capture  snap-back  phenomena  during  the  Iterative 
correction  phase  of  the  Incremental  staging  strategy  [6]. 

The  remaining  fornulation  of  fracture  energy  based  strain  softening 
follows  the  traditional  derivation  of  strain  controlled  elastic  plastic 
solids.  Resorting  to  an  associated  flow  rate  and  the  usual  additive 
decomposition  of  the  deformation  rate  we  find  In  the  absence  of 
elastic-plastic  coupling  the  following  evolution  laws  for  the 


Elastic  Response:  •  ■  E(*  -  «^)  (4) 

Fracturing  Response:  t,  ■  1  r  where  r  ■  —  (5) 
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The  consistency  condition  f  •  0  leads  to  the  strain-based  loading  condition 
for  progressive  cracking 
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when  the  trial  stress  rate  ij  ■  E  c  has  a  conponent  In  the  positive 
direction  of  the  normal  n.  The  fracture  parameter  follows  from 


i  .-tin 

f  Ef  ♦  n 
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where  the  fracture-based  strain  softening  modulus  was  defined  in  Eq.  2 
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This  formulation  for  a  single  active  crack  readily  extends  to  cracking  In 
two  and  three  orthogonal  crack  planes,  whereby  trial  stress  states  In  the 
vertex  region  require  the  simultaneous  solution  of  two  contiguous  fracture 
parameters  and  x£  along  the  line  of  simultaneous  plastic  flow  In  the 
case  of  nultlple  loading  surfaces. 

As  a  result  the  tangential  material  law  of  the  equivalent  elastic  strain 
softening  continuum  Is  defined  by 
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where  C-j.  -  £ 
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The  fracture  energy  basis  of  the  strain  softening  modulus  In  Eq.  8 
ensures  that  the  fracture  energy  released  within  each  elementary  volume 
remains  constant  Irrespective  of  the  size  of  the  finite  element.  He  note 
that  the  weak  requirement  for  local  uniqueness  [6]  infers  that  the 
denominator  in  Eq.  9  remains  positive,  l.e.. 


♦  n*  E  n  >  0  (10)  j 

j 

! 

This  restriction  infers  for  uniaxial  conditions  that  1 

i 

-Ef<  E  (ID 

i 

which  results  in  the  following  limitation  of  the  characteristic  dimension  j 

of  the  elementary  volume  } 
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This  condition  Imposes  severe  restrictions  on  the  choice  of  the  maximum  mesh 
size  hfflax  <  h#. 

MODEL  PROBLEM 

Let  us  examine  the  computational  aspects  of  the  composite  fracture  model 
when  applied  to  the  analysis  of  a  uniaxial  specimen  subjected  to  Increasing 
tension. 

In  the  quasistatic  case  of  a  uniform  tension  specimen  the  localization  of 
the  softening  zone  Is  a  problem  of  Its  own.  The  homogeneous  deformation 
problem  requires  first  deterministic  or  probabilistic  Identification  of  the 
crack  location  via  Initial  Imperfections  in  the  form  of  variations  of  the 
tensile  strength  or  the  cross-sectional  properties.  In  the  case  of  the 
discrete  fracture  approach  of  Ingraffea  and  Saouma  [12]  and  the  fictitious 
crack  model  of  Hllleborg  et  al.  [13],  a  real  or  fictitious  crack  element 
would  have  to  be  Incorporated  Into  the  mesh  lay-out  In  the  form  of  a 
discrete  interface  element  which  monitors  the  fracture  energy  release  rate 
directly  In  terms  of  the  ensuing  crack  opening  displacement.  In  the  case  of 
the  blunt  crack  approach  of  Bazant  et  al.  [14]  the  Initial  mesh  lay-out 
would  also  have  to  be  changed  In  order  to  accommodate  the  blunt  crack 
element  at  the  crack  location  which  now  has  finite  rather  than  zero  width  In 
the  order  of  three  times  the  aggregate  size. 

In  contrast  the  composite  fracture  model  does  not  require  remeshing  as 
long  as  strain  softening  Is  localized  within  a  single  element.  In  this  case 
the  proper  amount  of  fracture  energy  Is  released  Irrespective  of  the  loca¬ 
tion  of  the  “weak"  element  along  the  axis  of  the  tension  specimen.  In  view 
of  the  Incremental  format  of  strain-softening  computations  It  Is  however 
crucial  that  strain  softening  and  the  concomitant  deformation  trapping  will 
be  localized  within  a  single  element.  This  Infers  that  the  Iterative 
corrector  phase  of  the  solution  strategy  must  capture  elastic  snap  back 
within  elements  for  which  the  initial  predictor  step  might  erroneously 
signal  softening.  As  a  result  the  tangent  modulus  return  scheme  might  not 
only  lead  to  premature  difficulties  with  the  loss  of  positiveness  of  the 
tangent  stiffness  matrix  but  also  to  wrong  fracture  energy  release  predic¬ 
tions  when  no  special  provisions  are  made  for  "iterative  reversibility" 
within  the  load  increment.  In  fact  it  was  shown  recently  by  the  authors  [6] 


that  the  classic  initial  load  method  and  the  initial  modulus  return  scheme 
are  better  suited  for  capturing  strain  localization  in  a  single  element  when 
the  contiguous  elements  are  slnultaneously  unloading  elastically. 

Quasistatic  Study 

Fig.  2  illustrates  the  geometry  and  the  material  properties  of  the  model 
problem.  The  tension  specimen  is  subjected  to  monotonlcal 1y  Increasing  end 
displacements  which  result  In  a  uniform  tension  In  the  absence  of  Initial 
Imperfections.  The  finite  element  Idealization  of  half  the  bar  with  three 
elements  of  different  length  results  in  a  uniform  stress  distribution  in  the 
entire  bar,  and  only  numerical  round-off  eventually  leads  to  localized  soft¬ 
ening  In  the  short  element  A.  Fig.  3  depicts  the  overall  solution  in  terms 
of  the  nominal  stress  versus  axial  extension  diagram  of  half  the  bar,  the 
area  under  which  provides  a  measure  of  the  fracture  energy  release  rate. 
Reduction  of  the  cross-sectional  area  from  A  ■  1.0  to  Aw  »  0.99  leads  to 
predetermined  strain  localization  and  deformation  trapping  In  the  weak  ele¬ 
ment,  a  concept  which  has  previously  proposed  by  Prevost  and  Hughes  [15]  in 
order  to  Initiate  localization.  In  the  case  of  the  conposite  fracture 
model.  In  which  the  strain  softening  modulus  depends  on  the  size  of 
the  weak  element  h^  the  fracture  energy  release  rate  remains  constant 
when  either  of  the  three  elements  is  weakened.  In  the  case  of  the  tradi¬ 
tional  strain  softening  formulation  In  which  Ef  does  not  depend  on  mesh 
size,  the  corresponding  results  Indicate  that  the  ductility  and  thus  the 
underlying  fracture  energy  release  rate  strongly  depends  on  the  size  of  the 
weak  element  h^,  see  also  Crlsfleld  [16], 

The  axial  distribution  of  tensile  strains  Is  shown  In  Fig.  4  for  the 
constant  fracture  energy  and  the  constant  strain  softening  formulations. 

The  conposite  fracture  model  results  In  all  three  cases  In  Identical 
fracture  energy  release  values  and  thus  In  a  corresponding  decrease  of 
tensile  strain  level  when  the  size  of  the  weak  elements  Increases  from  h^ 

•  5  to  hg  »  IS  and  hc  •  30  (see  also  Fig.  3a).  The  constant  strain 
softening  model  in  contrast  leads  to  the  same  level  of  tensile  straining  In 
each  of  the  weakened  elements  and  thus  to  a  drastic  decrease  of  ductility 
and  corresponding  fracture  energy  as  the  size  of  the  weak  element  decreases 
(see  also  Fig.  3b).  In  the  limit  this  would  result  In  zero  fracture  energy 
expended  when  softening  localizes  In  an  element  of  infinitesimal  size.  This 
disturbing  aspect  of  the  traditional  strain-softening  approach  has  been 
rightly  criticized  In  the  past  by  Bazant  et  al.  [4,14].  We  might  add  that 
the  computational  difficulties  which  arise  with  maintaining  stable  control 
of  the  post-critical  response  predictions  are  exacerbated  by  Increasing 


brittleness  and  could  become  prohibitive  In  the  case  of  actual  snap-back 
conditions.  The  numerical  results  above  were  obtained  with  the  Incremental 
Initial  load  Iteration  method  In  which  the  Initial  elastic  stiffness  matrix 
remains  constant  during  the  entire  solution  phase.  There  was  no  need  for 
arc  length  stabilization  since  the  underlying  tangent  stiffness  maintained 
positive  definiteness  In  all  three  cases,  as  detailed  In  ref.  [6]. 

Dynamic  Study 

In  the  second  part  of  this  study  the  same  model  problem  was  analyzed  with 
the  aid  of  the  explicit  finite  element  code  DEGDYN  which  was  developed  by 
one  of  the  authors  [17]  for  the  dynamic  Investigation  of  concrete  dams.  In 
this  case  the  uniaxial  bar  was  subdivided  Into  9,  21  and  97  elements  of 
equal  length  and  subjected  to  the  constant  velocity  c  ■  0.50  at  both  ends, 
see  Fig.  5.  In  contrast  to  the  quasistatic  problem,  no  special  provision 
was  necessary  to  capture  localization  of  the  strain  softening  zone  in  the 
stress  wave  problem.  The  analytic  solution  was  recently  reviewed  by  Bazant 
and  Belytschko  [5]  for  which  a  well  posed  solution  was  found  Irrespective  of 
the  loss  of  hyperboliclty  of  the  underlying  characteristic  wave  equation  in 
the  softening  regime. 

Figs.  6  and  7  compare  the  numerical  results  of  the  composite  fracture 
model  with  those  of  the  traditional  strain  :oftening  approach  In  which  the 
softening  modulus  Is  kept  constant  Irrespective  of  the  mesh  size.  The 
resulting  stress  histories  and  the  associated  values  of  tne  fracture  energy 
release  rates  clearly  demonstrate  the  extreme  mesh  sensitivity  of  those 
quantities  when  E^  ■  const.  On  the  other  hand,  the  composite  fracture 
model  maintains  constant  fracture  energy  release  rates  without  resorting  to 
non-local  continuum  theories  along  the  line  of  the  recent  Imbricate 
continuum  proposal  of  ref.  [4]  for  resolving  this  basic  deficiency  of  the 
traditional  strain  softening  approach. 

Fig.  8  Illustrates  the  various  degrees  of  localized  deformation  trapping 
at  the  center  of  the  tension  specimen.  Clearly  It  Is  a  question  of  spatial 
resolution  to  capture  tensile  cracking  and  actual  separation  of  the  bar 
within  the  strain-softening  formulation  of  an  equivalent  continuum.  We 
observe  that  both  f emulations  reproduce  the  actual  separation  phenomenon 
quite  well.  In  fact,  the  conventional  strain  softening  approach  leads  to  a 
sharper  approximation  of  the  localized  strain  discontinuity  at  the  center  of 
the  bar  where  e  *  •  as  he  -*■  0.  On  the  other  hand,  the  composite  fracture 
model  enforces  constant  fracture  energy  release  rates  Irrespective  of  the 


particular  element  size  In  view  of  the  actual  separation  process, 
it  appears  questionable  if  non-local  continuum  theories  are  really  appro¬ 
priate  because  they  prevent  sharp  strain  localization  in  the  neighborhood  of 
the  actual  discontinuity  and  distribute  softening  and  thus  failure  over  a 
long-range  regime. 

CONCLUDING  REMARKS 

A  strain  softening  plasticity  model  was  developed  based  on  the  conposite 
fracture  concept  of  homogenizing  fracture  discontinuities  into  an  energy 
equivalent  continuum.  In  view  of  the  one-dimensional  nature  of  the  static 
and  dynamic  model  problems  the  main  thrust  of  the  exposition  was  focused  on 
mode  I  type  cracking.  The  composite  fracture  forrrulation  has  been  extended 
in  [8,19]  to  mixed  mode  fracture  in  the  form  of  splitting  compression  and 
shear  faulting  modes  of  failure  in  order  to  develop  a  unified  approach 
towards  softening  in  tension  and  compression. 
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Composite  Fracture 


Equivalent  Continuum 
Fig.  1  Composite  Fracture  Model 


Fig.  2  Quasistatic  Model  Problem 
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Fig.  3  Global  Response  Behaviour 
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Fig.  7  Mesh  Sensitivity  of  Fracture  Energy 
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